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20.  Abstract  (Cont.) 


The  effects  of  fluctuations  in  the  Coulomb  potential  due  to 
charged  impurities  in  high-purity  n-type  III-V  semiconductors  are 
examined  at  low  temperatures.  Assuming  that  charged  donors  and  accep¬ 
tors  are  randomly  distributed  at  high  temperatures,  we  conclude  that 
the  donors  are  selectively  filled  at  low  temperatures  leaving  non-random 
distributions  of  charged  and  filled  donors.  The  potential  fluctuations 
from  these  distributions  can  approximately  account  for  a  number  of 
experimental  observations  on  low-temperature  high-purity  GaAs  including 
the  apparent  decrease  of  donor  binding  energy  with  increasing  impurity 
concentration  observed  in  Hall  measurements,  the  sharpness  and  tempera¬ 
ture  dependence  of  the  ls-2p  transition  along  with  the  diffuse  ls-con- 
duction  band  edge  observed  in  photoconductivity  experiments,  and  the 
fact  that  experimentally  observed  photoconductivity  lineshapes  are 
narrower  than  those  previously  predicted. 

Static  strains  in  piezoelectric  semiconductors  give  rise  to  an 
electric  field  or  potential  which  can  have  an  effect  on  the  electrical 
properties  of  the  material.  We  have  calculated  the  electric  potential 
due  to  the  strain  field  arising  from  a  random  distribution  of  point 
defects.  This  potential  contributes  a  term  to  the  mobility  that  is 
proportional  to  T1/2  and  a  Hall  factor  of  1.10.  Crude  estimates  of 
strain  strengths  indicate  that  this  scattering  mechanism  may  contribute 
significantly  to  the  mobility  of  electrically  rather  pure  III-V 
semiconductors  below  room  temperature  when  neutral  impurity  concentra¬ 
tions  are  greater  than  lO^-Scm-3.  The  mechanism  may  also  constitute 
a  dominant  one  in  the  mobility  of  some  III-V  alloys  at  fairly  low  tem¬ 
peratures.  The  existence  of  strain  induced  electric  potentials  also 
provides  at  least  a  possible  mechanism  whereby  different  donors  can 
have  different  lineshapes  as  measured  in  photoconductivity  experiments. 

A  model  of  Cr  in  GaAs  which  is  consistent  with  a  large  body  of  ex¬ 
perimental  data  has  been  developed.  It  relies  on  recent  spectroscopic 
models  and  on  our  interpretation  of  redistribution  and  electrical  data, 
all  of  which  indicate  the  existence  of  Cr  complexes.  The  existence  of 
rapidly  diffusing  interstitial  Cr  donors  is  assumed  and  justified. 

The  model  offers  a  unified  picture  of  the  effects  of  implantation  on 
the  Cr  profile.  It  contains  mechanisms  for  compensation  and  redistribu¬ 
tion,  which  offer  an  explanation  of  the  semi-insulating  properties  of 
Cr  doped  GaAs  and  of  the  two  apparently  incompatible  classes  of  diffu¬ 
sion  and  anneal  data.  The  redistribution  depends  on  how  the  Cr  was 
incorporated  and  on  the  vacancy  concentration  profiles.  A  study  of 
representatives  of  the  two  classes  of  redistribution  data  allows  us  to 
estimate  a  lower  limit  of  interstitial  Cr  diffusion  constant  and  of 
the  vacancy  diffusion  lengths  in  GaAs. 

During  the  annealing  of  ion-implanted  Cr-doped  GaAs,  Cr  often  re¬ 
distributes  and  accumulates  at  the  surface.  Although  this  behavior  has 
been  attributed  to  strain  fields  and  other  mechanisms,  the  widths  of 
these  accumulation  regions  suggest  that  electric  fields  due  to  surface 
states  are  a  limiting  factor  in  Cr  redistribution.  For  this  reason  we 
have  developed  a  thermodynamic  model  for  Cr  redistribution  which  takes 
into  account  the  electric  field  due  to  surface  states.  A  qualitative 
fit  to  SIMS  data  on  annealed  unimplanted  GaAs  samples  can  be  obtained 
with  this  model.  We  have  also  used  applied  voltages  during  annealing 
to  modify  the  amount  of  band  bending  and  Cr  buildup  at  the  surface. 

This  experiment  indicates  that  the  accumulated  ions  are  positively 
charged.  We  conclude  from  these  experiments  that  electric  fields  play 
a  significant  role  in  the  redistribution  of  Cr  at  GaAs  surfaces. 
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20.  Abstract  (Cont.) 

In  this  report  we  present  a  model  of  impurity  redistribution 
during  epitaxial  growth  of  semiconductors  due  to  diffusion,  drift 
in  the  built-in  electric  field,  and  the  growth  itself.  Examples 
show  unexpected  impurity  profiles  for  the  minority  species  and  indi¬ 
cate  how  undesirable  conductivity  regions  can  occur. 

The  incorporation  of  Group  IV  impurities  as  donors  and  as  accep¬ 
tors  in  high  purity  epitaxial  GaAs  has  been  investigated  using  photo- 
thermal  ionization  spectroscopy  and  variable  temperature  photolumin¬ 
escence  to  detect  donors  and  acceptors  respectively.  Samples  from 
several  sources  of  high  purity  LPE ,  AsClo-VPE,  ASH3-VPE,  MOCVD,  and 
MBE  grown  GaAs  were  measured  to  establish  the  typical  residual  impuri¬ 
ties  present  and  their  relative  concetnrations .  For  AsH3-VPE,  MOCVD 
and  MBE  GaAs,  impurity  incorporation  data  are  presented  as  a  function 
of  III/V  ratio.  The  relative  incorporation  of  amphoteric  impurities 
as  donors  and  acceptors  is  compared  with  the  model  of  Teramoto  (1972) 
for  LPE  and  its  extension  to  ASCI3-VPE  by  Ashen  et  al.  (1975) . 

Photothermal  ionization  measurements  on  Si  doped  GaAs  grown  by 
molecular  beam  epitaxy  (MBE)  have  indicated  that  the  impurity  peak 
previously  assigned  to  Si  was  incorrect.  Data  leading  to  the  new 
identification  are  presented,  and  the  results  leading  to  the  earlier 
identification  are  reexamined.  Implications  of  the  new  identification 
on  the  importance  of  Si  as  a  residual  donor  in  GaAs  grown  by  various 
techniques  are  discussed. 

Photothermal  ionization  spectroscopy  has  also  been  used  to  deter¬ 
mine  the  residual  donor  species  present  and  their  relative  concentra¬ 
tions  in  the  highest  purity  MBE  n-GaAs  yet  reported.  Data  are  presented 
for  samples  grown  in  two  different  MBE  growth  reactors;  one  using 
elemental  As  and  the  other  using  cracked  AsH3  as  the  arsenic  source. 

Jn  spite  of  the  substantial  differences  between  growth  systems,  the 
donor  backgrounds  are  quite  similar. 

To  directly  examine  possible  impurity  interactions  or  complexing 
in  GaAs,  multiply-doped  epitaxial  layers  were  prepared.  Samples 
doped  both  homogeneously  and  inhomogeneous ly  with  H2O  and  H2S  exhibit 
effects  which  indicate  that:  (1)  H^O  produces  free  carrier  compensa¬ 
tion  and  deep  donor  behavior;  and  (2)  H2O  affects  the  incorporation 
and/or  diffusion  of  sulfur.  Experiments  are  being  performed  to 
obtain  more  insight  into  these  interactions . 
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1.  RESEARCH  OBJECTIVES 


The  purpose  of  this  work  is  to  investigate  the  inter¬ 
actions  among  impurities,  defects,  and  complexes  in  GaAs 
which  adversely  affect  the  yield  and  performance  of 
high-speed  GaAs  integrated  circuits.  To  achieve  this 
objective  the  following  experimental  approach  is  being 
employed:  Impurities  are  introduced  into  GaAs  by  gas-phase 
doping  during  epitaxial  growth  and  by  ion- implantation 
into  bulk  samples.  Annealing  is  performed  under  controlled 
atmospheres  with  applied  fields.  The  resulting  samples 
are  characterized  by  a  combination  of  low  and  high  temper¬ 
ature  resistivity  and  Hall  measurements,  differential 
capacitance  measurements,  far  and  near  infrared  photo¬ 
conductivity  measurements,  photoluminescence,  electron 
spin  resonance,  and  secondary  ion  mass  spectroscopy. 

These  data  are  then  analyzed  with  impurity  incorporation, 
redistribution,  complex  formation,  and  other  models. 

In  the  theoretical  part  of  the  program  various  aspects 
of  the  experimental  results  are  investigated  to  obtain 
new  insights  into  impurity  and  defect  interactions. 
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2.  PHOTOCONDUCTIVITY  SPECTRA  FOR  HYDROGENIC  DONOR  COMPLEXES 

Although  the  study  of  the  spectra  of  shallow  donors  in  semiconductors 
is  a  rather  old  subject,  many  problems  of  understanding  the  linewidths, 
lineshapes,  and  even  the  line  positions  remain  unsolved.  This  is  somewhat 
unfortunate  since  photoconductivity  experiments  cannot  only  detect  shallow 
donor  Impurities  that  have  densities  of  10  /cm  or  less,  but  can  also 
yield  resolvable  lineshapes.  With  such  a  fantastically  sensitive  spectro¬ 
scopic  tool,  one  might  hope  to  obtain  a  wealth  of  information  about  the 
structure  and  environment  of  shallow  donor  impurities. 

One  major  difficulty  is  that  different  species  of  shallow  donors 
ought  to  yield  identical  lineshapes  whose  positions  are  controlled  by 
central  cell  corrections  and  whose  intensities  are  proportional  to  the 
densities  of  the  species.1  As  has  been  discussed  in  the  literature,  this 
is  not  the  case  quantitatively  and  is  sometimes  not  even  the  case  qualita¬ 
tively.  For  example,  linewidths  for  different  species  in  GaAs  are  usually 

different1”^  and  often  even  lineshapes  are  grossly  different.4"6 
2.1  EFFECTIVE  MASS  APPROXIMATION 

Towards  a  partial  understanding  of  these  discrepancies  we  investi¬ 
gate  the  electrostatic  effects  on  the  energy  levels  of  a  hydrogenic 
shallow  donor  complex  in  the  presence  of  an  applied  magnetic  field. 

For  our  purposes  here,  a  simple  shallow  donor  is  described  by  an 
effective  mass  Hamiltonian  where  the  potential  is  due  to  a  point 
charge  of  magnitude  e.  A  complex  (or  complex  shallow  donor)  is  also 
described  by  an  effective  mass  Hamiltonian  but  the  potential  is  due 
to  a  charge  density  p(r)  whose  Integrated  weight  is  e  but  which  in 
general  is  not  the  charge  density  of  a  single  point  charge. 
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The  Hamiltonian  for  an  electron  in  the  field  of  a  donor  and  a  magnetic 
field  I  which  defines  the  z-axis  is 

H  *  (p2/2m*)  “>c(xpy-ypx)  +|-m*^(x2+y2)  -  (e/e^JdV  p(r*  )/|r-r'  j , 
Jd3r*  p(r')  =  e  (1) 

where  wc  »  eB/m*c  is  the  cyclotron  frequency.  We  are  assuming  that  the 
effective  mass  approximation  is  valid  and  thus  m*  and  e0  are  the  effective 
mass  and  dielectric  constant  of  the  medium.  We  shall  further  assume  that 
p(r)  Is  zero  outside  of  a  volume  whose  dimensions  are  small  compared  to  a 
Bohr  radius.  The  effects  of  charged  impurities  many  Bohr  radii  from  the 
center  of  the  charge  distribution  can  be  treated  as  perturbations. 

The  effective  mass  approximation  yields  excellent  results  In  a  wide 
variety  of  semiconductors7  and  Its  validity  will  not  be  discussed  here. 

Of  course,  in  a  real  material,  p(r)  Is  not  the  actual  charge  distribution 
but  Is  the  actual  charge  distribution  minus  the  charge  distribution  of  the 
perfect  lattice.  Further,  the  use  of  a  dielectric  constant  In  and  very 
near  the  charge  distribution  is  clearly  not  valid.  The  details  of  the 
interactions  in  this  volume  are  contained  in  the  central  cell  correction 
to  the  energy  of  the  Is  state  which  we  shall  assume  to  be  given.  These 
corrections  are  quite  small  In  GaAs  because  the  effective  Bohr  radius  is 
99  A. 
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2.2  ANALYSIS 

In  this  section  we  outline  the  calculator  of  the  effects  of  a  complex 
on  the  energy  levels  of  a  shallow  donor.  Initially  we  ignore  the  effects 
due  to  any  remote  charges.  In  reduced  units,  the  zero  order  Hamiltonian, 

Hq,  Is  that  for  a  simple  donor  at  the  origin, 

H0  =  -V2  +  (y/1)0/3<J>)  +  \yZp2  -  2/r  ,  (2) 

where  lengths  and  energies  are  taken  in  units  of  the  effective  Bohr  radius, 

O  O  A  O  O 

aQ  *  eQ  ir/m*e  ,  and  the  effective  Rydberg,  R  =  m*e  72^15.  The  dimensionless 

2  2  2 

magnetic  field  strength  is  y  where  y  -  -fko _/2R,  p  =  r  sin  0  ,  and  the 
magnetic  field  $  defines  the  z  direction.  In  what  follows,  we  shall 
always  assume  y  Is  much  larger  than  any  other  (dimensionless)  perturbation. 

We  define  the  complex  by  a  set  of  charges  Z.e  at  the  points  r^  where 

£z,  -1  .  (3) 

1  1 

(If  the  charge  distribution  Is  a  continuum,  the  appropriate  sums  can 
trivially  be  converted  to  Integrals.)  Thus  H',  the  perturbation  due 
to  the  complex,  can  be  written  as  a  sum  of  terms,  one  from  each  member 
of  the  complex. 

H1  - Z »4 
1  1 

V,  -  -EZ^Ir-r,!"1  -  r'1)  .  (4) 

At  this  point  it  Is  important  to  note  that  the  origin  of  the  charge 
distribution  is  Irrelevant.  Thus,  if  we  translate  all  of  the  charges 
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where  the  Y^fi)  are  the  usual  spherical  harmonics  and  fi  and  fi'  refer  to 

the  solid  angles  for  r  and  r-  respectively.  As  noted  earlier,  we  have 

assumed  that  the  charge  distribution  has  dimensions  d  and  that  d  «  1 

in  reduced  units.  Assuming  that  the  origin  in  the  problem  is  chosen  to 

lie  within  a  distance  d  of  the  charge  distribution,  the  contribution  to 

2 

the  energy  of  an  s  state  from  r<  r^  will  be  of  order  d  and  will  be 

independent  of  the  angles  8^.  The  corresponding  energies  of  other  states 

will  be  proportional  to  the  fourth  or  higher  powers  of  d.  Such  contributions 

can  be  lumped  into  a  central  cell  correction  and  will  temporarily  be  Ignored 

while  we  concentrate  on  the  contributions  from  Eqs.  (4)  and  (6a). 

Although  the  final  answers  must  be  independent  of  origin,  some 

choices  are  more  convenient  than  others.  The  perturbations  described  by 

Eq.  (6a)  contains  a  sum  over  i  and  the  tth  term  is  proportional  to  d  . 

If  we  chose  the  origin  to  eliminate  the  l  -  1  term,  then  first  order 

2 

perturbation  theory  will  yield  the  correct  answer  to  order  d  .  Any  other 

choice  will  require  second  order  perturbation  theory  to  obtain  answers  to 
2 

order  d  .  This  choice  is  accomplished  by  choosing  the  origin  so  that 


r 
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EziV»  •  <z> 

that  is,  so  that  the  dipole  moment  of  the  charge  distribution  vanishes. 

2 

Now,  to  order  d  ,  H'  can  be  written  with  only  the  £» 2  terms  as 

H'  ’-(8"/5)3^Z1r?Y2ra(i!1)jWn)/r3  •  (8) 

Only  the  ma0  term  in  Eq.  (8)  will  contribute  in  lowest  order 
perturbation  theory.  The  energy  shift  of  the  state  x  can  be  written  as 


aex  a  -^(xM^) 

(9a) 

ctg(x)  *  <x|(3cos2B  -  l)/r3|x> 

(9b) 

QC«i )  3  j  ^  2^r2(3cos2  6f  -  1)  . 

(9c) 

The  r^  in  Eq.  (9c)  are  determined  by  Eq.  (7)  and  6i 

refers  to  the  ith  charge 

in  the  coordinate  system  defined  by  the  magnetic  field.  Equation  (9c)  will 
be  presented  in  a  more  convenient  form  in  the  next  section.  Before  presenting 
our  calculations  of  <Xq(x)  for  the  Is  and  2p_  states,  we  wish  to  briefly 
comnent  on  the  coupling  of  terms  arising  from  the  complex  with  terms 
describing  the  electric  field  from  remote  charges. 

The  Hamiltonian  describing  an  electron  interacting  with  an  electric 
field  t  generated  by  remote  charges  is 

Hs  a  -et«r/eor3  .  (10) 

Since  this  term  has  an  odd  parity,  bound  states  (in  the  presence  of  a 
magnetic  field)  are  unaffected  by  H$  to  first  order.  They  are,  of  course, 
affected  to  second  order  and  this  Hamiltonian  gives  rise  to  the  second 


order  Stark  shift.  As  has  been  extensively  discussed  in  the  literature, 
different  donors  are  under  the  influence  of  different  electric  fields 
generated  by  remote  defects.  This  leads  to  a  characteristic  lineshape 
that  is  usually  observed  in  photoconductivity  experiments. 

Since  in  general  there  will  be  a  "length"  vector  3  associated  with 
the  charge  distribution  of  a  complex,  one  might  have  expected  a  term  in 
the  complex  Hamiltonian  that  was  linear  in  d  and  had  odd  parity.  If  such 
a  term  existed,  the  cross  term  of  it  and  the  Stark  term  in  second  order 
perturbation  theory  would  lead  to  an  energy  shift  proportional  to  Ed; 
i.e.,  an  effective  first  order  Stark  shift.  In  fact,  such  a  term  did 
exist  in  the  Hamiltonian  for  a  complex  before  we  transformed  it  to  zero 
via  Eq.  (7).  If  we  had  not  transformed  this  dipole  term  to  zero,  every 
energy  level  would  have  been  shifted  by  a  term  proportional  to  Ed. 

However,  the  multiplicative  coefficient  would  be  the  same  for  all  states 
and  thus  the  term  would  cancel  out  for  any  energy  difference. 

A  term  linear  In  E  whose  coefficient  depends  on  the  charge  configura¬ 
tion  could  be  useful  In  trying  to  explain  anomalous  linewidths  from 
different  donor  species.  However,  our  argument  does  not  rule  out  such 

Q 

a  term.  It  merely  rules  It  out  for  a  medium  described  by  only  an  isotropic 

g 

dielectric  constant.  Terms  depending  on  other  less  isotropic  properties 
of  the  media  may  yield  such  terms. 
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2.3  CALCULATIONS 

The  energy  shift  of  a  shallow  donor  state  due  to  the  electrostatic 
effects  of  a  complex  is  given  by  Eqs.  (9)  in  the  reduced  units  discussed  at 
the  beginning  of  Sect.  II.  For  GaAs,  lengths  are  measured  in  terms  of  the 

O 

Bohr  radius  aQ  *  99  A,  energies  are  measured  in  terms  of  the  Rydberg 
R  =  5.77  meV  =  46.5  cm“\  and  y  is  the  dimensionless  magnetic  field  strength 
where  B  (in  kG)  =  65.6  y.  In  this  section  we  present  and  discuss  the  results 
of  our  calculations  for  the  energy  difference  of  the  Is  and  2p_  states 
only  since  this  is  the  most  important  spectroscopic  transition. 

The  wave  functions  of  the  Is  and  2p_  states  for  the  zero  order 
Hamiltonian,  Eq.  (2),  describing  an  hydrogenic  atom  in  a  magnetic 
field  were  obtained  by  variational  calculations.  The  trial  wave 
functions  used  were 

*  A  exp(-ar-b2p2  -  c2z2)  ,  (11a) 

*2p  ^  *  Ap  exp(-ar-b2p2  -  c2z2  -  i<*>)  ,  (lib) 

where  a,  b,  and  c  were  varied  in  each  case  to  minimize  the  energy.  One 
reason  for  the  choice  In  Eqs.  (11)  is  that  all  radial  Integrals  can  be 
expressed  In  terms  of  the  parabolic  cylinder  function  and  only  the  6 
Integrals  had  to  be  performed  numerically.  Thus  the  energies  and  all 
derivatives  -with  respect  to  a,  b,  and  c  are  easily  computed  and  one  can 
minimize  the  energies  with  very  few  iterations.  The  energies  obtained 
are  not  quite  as  low  as  have  been  obtained  with  better  trial  wavef unctions. 
For  example  we  obtain 


I 


Eh(Y)  -  el(y)  \ 

c  rm —  p/a\ 


<  2  x  10"3 


(12) 
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for  all  values  of  y  less  than  10.  In  this  expression  E^(y)  is  our  computed 
energy  as  a  function  of  y  for  the  Is  state,  E^(y)  is  the  energy  as  a  function 
of  Y  as  obtained  by  Larson^ using  better  trial  wave  functions,  and  E(0)  is 
the  Is  energy  in  zero  magnetic  field. 

Using  the  above  wavefunctions  we  have  calculated  ap(x)  where  x  refers 
to  the  Is  and  2p_  states.  Since  it  is  the  energy  shift  between  the  Is  and 
2p_  states  that  is  important,  we  write 

AE(y)  -  6(y)Q(^). 

8(y)  *  aQ(ls)  -  ag(2p_)  (13) 

where  AE(y)  is  the  shift  in  E2p  -Els  due  to  the  complex.  Q(n.),  given  by 
Eq.  (9b)  in  reduced  units,  depends  only  on  the  structure  of  the  complex 
and  its  orientation.  At  y  *  0, . 

6(0)  =  1/60  .  (14) 

The  quantity  6(y)  is  plotted  vs  y  in  Figure  1.  For  purposes  of  comparison, 
the  magnetic  field  dependence  of  the  central  cell  correction  E  (y)  is  also 
plotted  on  Figure  1  where 

Ec(y)  *  l/|^ls(r«0) |2  .  (15) 

There  are  two  features  that  distinguish  the  complex  shift  from  the 

•  • 

central  cell  shift.  First  of  all,  the  complex  shift  changes  much  more 
rapidly  with  magnetic  field.  Secondly,  and  more  important,  the  complex 
shift  depends  on  the  orientation  of  a  complex  with  respect  to  the  magnetic 
field.  As  can  be  seen  from  Fig.  1,  the  effects  of  a  complex  as  a  function 
of  field  orientation  can  best  be  observed  at  highest  magnetic  fields. 
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As  an  example  consider  a  complex  consisting  of  two  charges  Z-j  *  2 
and  Z2S-1  separated  by  a  dimensionless  distance  a.  Using  Eqs.  (7)  and 
(9c)  we  obtain 

Q(fl)  *  a2(3cos20  -1)  (16) 

where  6  is  the  angle  between  the  magnetic  field  and  a  line  connecting  the 
two  charges.  Suppose,  for  instance,  such  a  complex  could  be  oriented  along 
any  one  of  the  four  independent  <111>  axes.  If  $  was  pointing  along  one 

i 

of  the  three  independent  <100>  axis,  Q  would  be  zero  in  all  cases.  On  the 

other  hand.  If  3  were  pointing  along  one  of  the  <111>  axes,  complexes 

2 

pointing  along  that  axis  would  have  Q  *  2a  and  complexes  pointing  along 

2  2 

the  other  three  axes  would  have  Q  *  -ya  .  Numerically,  for  GaAs  in  a 
magnetic  field  of  65.6  kG,  the  above  complex  with  length  of  a  equal  to 
one  lattice  spacing  could  give  rise  to  a  shift  that  moves  over  a  range 
of  0.0816  cin~^  as  a  function  of  magnetic  field  angle. 

There  Is,  of  course,  no  guarantee  that  a  given  complex  of  net  charge 
one  is  a  shallow  donor.  The  question  of  whether  It  Is  or  not  is  beyond 
the  scope  of  this  paper  and  will  not  be  discussed  further.  Besides  sub¬ 
stitutions  of  two  or  more  atoms  to  form  a  complex  there  could  be  configura¬ 
tions  where  a  single  substitution  plus  a  shift  in  position  forms  one. 

For  example,  suppose  a  C  atom  replaces  an  As  atom  but,  because  of  its 
small  size,  the  C  atom  moves  slightly  from  the  normal  As  position.  Such 
a  complex  might  be  described  by  a  charge  Z  at  the  C  site  and  a  charge 
1-Z£  at  the  old  As  site  where  Zc  is  an  effective  charge  for  the  C. 
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Figures  2  and  3  exhibit  photoconductive  response  data  of  doped  GaAs 
epitaxial  layers.^  The  arrows  in  the  figures  point  to  structure  that 
could  be  interpreted  as  resolved,  partially  resolved,  and  unresolved 
splittings due  to  a  complex.  This  data  exists  only  at  one  orientation 
and  therefore  it  is  not  known  whether  the  structure  shifts  with  magnetic 
field  angle.  Thus  one  could  also  obviously  Interpret  the  data  as  structure 
due  to  different  chemical  species. 
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Fig.  2.  High  resolution  photoconductivity  spectra  of  a  GaAs  epitaxial 
layer  at  a  magnetic  field  of  50  kG  and  a  temperature  of  42  K. 
(See  ref.  1.)  The  arrows  point  to  structure  that  could  be 
interpreted  as  resolved  and  partially  resolved  splittings 
due  to  complexes. 
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High  resolution  photoconductivity  spectra  of  a  GaAs  epitaxial 
layer  at  a  magnetic  field  of  50  kG  and  a  temperature  of  4.2  K 
(See  ref.  1.)  The  arrow  points  to  structure  that  could  be 
interpreted  as  an  unresolved  splitting  due  to  a  complex. 
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COULOIR  FLUCTUATIONS  FROM  NON-RANDOM  DONOR  DISTRIBUTIONS 


For  years  It  has  been  known  that  a  random  distribution  of  equal 

concentrations  of  fixed  positive  and  negative  charges  will  yield  Infinite 

fluctuations  In  the  Coulomb  potential  In  the  limit  of  an  infinite  volume J 

These  fluctuations  can  be  damped,  of  course,  by  screening  due  to  additional 

mobile  charges  or  by  a  rearrangement  of  the  fixed  charges  themselves.  In 

light  of  this,  consider  a  high-purity  n-type  semiconductor  at  low  temperatures. 

As  the  temperature  Is  lowered,  the  number  of  free  carriers  Is  reduced  toward 

2  • 

zero,  and  thus  the  screening  or  Debye  radius  tends  toward  Infinity.  Since 
one  does  not  normally  think  of  the  charged  donors  and  acceptors  as  being 
mobile,  one  would  naively  expect  the  Coulomb  fluctuations  to  become  huge 
and  finally  infinite. 

However,  there  is  another  available  mechanism  for  partially  screening 
these  fluctuations  that  consists  of  the  spatially  selective  filling  of 

*  l 

donor  states  rather  than  the  usual  picture  of  filling  these  states  randomly. 

That  Is,  the  Coulomb  fluctuations,  and  of  course  the  free  energy,  will  be  | 

lowered  by  selectively  filling  the  charged  donor  states.  The  purpose  of 

j 

this  section  is  to  theoretically  Investigate  this  process  Including  the  j 

magnitude  and  consequences  of  the  remaining  finite  fluctuations.  We  shall 
limit  our  discussion  and  arguments  to  n-type  III-V  semiconductors  with  low 
Impurity  concentrations  of  shallow  donors  and  acceptors.  For  the  purposes 
of  this  report  low  Concentrations  mean  "d"  ao  «  1  where  nd  Is  the  number  j 

density  of  donors  and  aQ  Is  the  Bohr  radius  of  a  shallow  donor.  Thus,  for 
GaAs  with  aQ*99  A,  our  limit  would  Include  samples  with  nd  <  1015  cm’3. 

These  concentrations  are  low  enough  so  that  Impurity  band  conduction  [ 

3 

becomes  very  difficult. 

i 


Given  a  random  distribution  of  donors  and  acceptors,  finding  the 
actual  distribution  of  filled  and  empty  donors  from  first  principles  is 
an  extremely  difficult  problem  that  will  not  be  attempted  here.  Instead 
we  approach  the  problem  In  a  much  more  phenomenological  (and  less  rigorous) 
way  as  follows.  We  assume  that  each  acceptor  can  be  paired  with  an  unfilled 
or  charged  donor.  This  pairing  radius  Is  described  by  the  probability 
distribution  function  for  being  able  to  construct  a  charge  neutral  volume 
around  an  acceptor  with  that  radius.  The  distribution  of  Coulomb  potentials 
can  then  be  obtained  and  used  to  deduce  various  properties  of  the  system. 

The  picture  that  emerges  from  our  analysis  . of  the  system  at  low 
temperatures  Is  as  follows.  There  must,  of  course,  be  fluctuations  In 
the  Coulomb  potential  due  to  the  charged  Impurities.  However,  these 
fluctuations  are  quite  different  from  those  which  one  would  obtain  from 
a  random  distribution  of  charged  Impurities.  This  non-random  distribution 
of  fluctuations  due  to  the  selective  filling  of  donors  has  two  Important 
aspects.  In  the  first  place,  the  average  potential  energy  that  electrons 
feel  from  charged  donors  and  acceptors  Is  negative  In  regions  of  the 
crystal  near  filled  donors  even  though  the  average  Coulomb  potential 
throughout  the  crystal  Is,  of  course,  zero.  Secondly,  there  are  finite 
fluctuations  In  the  potential  energy  of  an  electron  about  this  mean  value 
that  are  smaller  than  the  average  fluctuations  In  the  crystal  as  a  whole. 

Further,  our  analysis  provides  at  least  an  approximate  explanation 

of  a  number  of  puzzling  phenomena  observed  In  relatively  pure  n-type  GaAs 

at  low  temperatures.  These  Include  the  apparent  decrease  In  E^,  the  donor 

4  5 

binding  energy,  with  Increasing  nd  as  observed  In  Hall  measurements  * 
and  the  decrease  In  mobility  of  many  samples  below  about  10  K.®  It  also 
Includes  the  sharpness  and  temperature  dependence  of  the  ls-2p  transition 


along  with  the  lack  of  a  sharp  1 s-conduction  band  edge  observed  in  photo- 

5  6 

conductivity  experiments.  Finally,  it  gives  an  explanation  for  the 

fact  that  experimentally  observed  photoconductivity  lineshapes  are  narrower 

7  8 

than  those  predicted  theoretically  *  using  the  familiar  Holzmark  distribu¬ 
tion  of  electric  fields  in  a  sample. 

There  have  been  a  number  of  previous  explanations  of  some  of  the  above 

5 

phenomena  that  depend  on  the  banding  of  impurity  states.  However,  as  will 
be  discussed,  the  Coulomb  fluctuations  are  much  greater  than  typical  overlap 
integrals  at  low  concentrations  of  impurities  which  should  preclude  banding. 
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3.1  BASIC  MODEL 

First,  It  Is  Instructive  to  pursue  the  concept  of  the  random  filling 
of  charged  donor  states  In  somewhat  more  detail.  Thus,  we  consider  the 
distribution  of  electric  potentials  at  a  given  point  due  to  a  random 
distribution  of  point  charges  with  a  density  n^.  The  potential  at 
the  origin  due  to  single  positive  charge  at  r  is  the  screened  Coulomb 
potential 

<Kr)  *  (e/eQr)exp(-r/rs)  ,  (1) 

where  e  Is  the  static  dielectric  constant  of  the  medium  and  r.  Is  a 
o  s 

screening  radius.  We  assume  that  the  system  Is  charge  neutral  so  that 
the  number  of  positively  charged  Impurities  Is  equal  to  the  number  of 
negatively  charged  Impurities  plus  the  number  of  conduction  electrons. 
Finally,  we  take  the  low-concentration  limit  where  n^  Is  much  smaller 
than  the  density  of  lattice  points. 

The  calculation  of  the  distribution  of  potentials  for  the  above 
system  Is  easily  obtained  as  a  special  case  of  the  calculation  de  - 
scribed  below.  For  our  purposes  here  it  Is  sufficient  to  characterize 
the  Cdttlomb  potential  fluctuations  by  their  second  moment  and  one  obtains 

*  <$2>  ■  2irntrs  e2/e2  •  (2) 

In  order  to  apply  this  to  an  n-type  semiconductor  we  take  r$  to  be  the 
Oebye  screening  radius, 

rs  ■  (kTe0/4irne2)'*  ,  (3) 

where  n  Is  the  density  of  conduction  electrons  (assumed  nondegenerate) 
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and 

nt  “  ^2na  +  *  (4) 

where  n.  Is  the  density  of  charged  acceptors  and  n  +  n  Is  equal  to  the 

a  a 

density  of  charged  donors.  At  this  point  one  can  see  that  as  the  density 
of  conduction  electrons  tends  toward  zero,  the  Debye  radius  and  <frr  tend 
toward  Infinity. 

To  proceed  further  we  assume  a  crude  model  where  the  charged  donors 
have  only  one  bound  state  with  binding  energy  E^  and,  using  effective  mass 
theory, 

Ed  *  ®2/2aoeo  *  (5) 

Since  the  potential  energy  fluctuations  at  different  donor  sites  span  a 
width  In  energy  of  order  e$_  we  approximate  these  fluctuations  by  a  band 

r 

of  width  UQ  where 

Ufl  *  ae*p  ,  (6) 

and  a  Is  a  constant  of  order  one.  Thus,  assuming  that  the  conduction 
band  edge  lies  at  the  potential  energy,  the  bound  donor  states  have 
energies  with  respect  to  the  lowest  part  of  conduction  band  edge  of 

E(x)  -  -Ed  +  U0  x  .  (7) 

where  x  Ties  between  zero  and  one.  From  Eqs.  (2)  through  (6)  we  obtain 

(U0/Ed)  •  a[(87rkT/Ed)(n^/n)r*  ,  (8) 

Now,  as  n  decreases,  UQ  Increases.  However,  If  UQ  becomes  too  large, 
many  of  the  charged  donor  states  will  no  longer  bind  an  electron  because 
E>0  In  Eq.  (7).  This  will  self-cons  1  stent ly  stabilize  the  Coulomb 
fluctuations  at  Ug^Ej.  One  can  obtain  quantitative  results  from  this 


model  by  performing  the  standard  calculation  for  n  using  Eq.  (7)  with  x 
taking  on  all  values  between  zero  and  one  with  equal  probabilities.  The 
results  of  this  calculation  show  that  If  »  kT,  the  number  of  conduction 
electrons  will  not  decrease  exponentially  as  the  temperature  Is  lowered  but 
will  be  proportional  to  kT.  Since  this  argument  does  not  depend  on  the 
density  of  donors  or  acceptors.  It  predicts  Impurity-Induced  conduction 
electrons  at  low  temperatures  for  any  concentration  of  impurities.  Further 
n  will  be  proportional  to  kT  and  thus  the  effective  measured  In  Hall 
measurements  would  be  zero. 

At  least  for  small  enough  concentrations  of  charged  Impurities  the 
above  scenario  does  not  occur.  In  our  picture  we  still  assume  that  the 
original  distribution  of  donors  and  acceptors  Is  random  since  at  the 
temperatures  at  which  materials  are  usually  synthesized  there  are  plenty 
of  conduction  electrons  available  to  screen  the  fluctuations.  However, 
as  the  temperature  Is  lowered  and  donor  states  start  to  fill,  they  will 
not  be  filled  at  random.  Instead  they  will  fill  In  such  a  manner  as  to 
lower  the  Coulomb  fluctuations  and  also  the  energy  of  the  system. 

We  now  consider  In  more  detail  a  system  with  a  density  nrf  of  shallow 
donors,  a  density  n.  of  charged  acceptors,  and  a  compensation  ratio. 


K  -  na/nd  ,  (9) 

less  than  one.  Both  types  of  charged  Impurities  are  distributed  at 
random.  If  at  low  temperatures  all  of  the  conduction  electrons  are 
frozen  out,  then  there  must  be  n,  unfilled  or  charged  donors  and  (K“*-l)n_ 

a  « 

filled  or  neutral  donors.  The  large  fluctuations  discussed  earlier  are 
due  to  the  long-range  nature  of  the  Coulomb  potential  and  arise  from 
regions  of  the  sample  that  are  not  charge  neutral.  They  are  In  no  way 
due  to  the  divergence  of  the  Coulomb  potential  at  the  origin. 
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In  order  to  understand  the  length  scale  of  these  fluctuations  we  ask 
the  following  question:  Given  a  charged  acceptor  at  the  origin,  what  Is 
the  probability  P(r)  of  being  able  to  enclose  that  acceptor  In  a  charge* 
neutral  sphere  of  radius  r.  Besides  the  acceptor  at  the  origin  the  sphere 
may  enclose  k  other  acceptors  and  n  donors.  Since  the  number  of  charged 
donors  must  be  equal  to  the  number  of  acceptors,  all  Integral  values  of 
n  and  k  with 

n-1  >  k  >  0  (10) 

are  allowed.  This  Is  a  well  defined  mathematical  problem  whose  solution 
will  yield  some  Information  about  how  small  the  Coulomb  fluctuations  can 
be  made. 

The  solution  to  this  problem  Is  obtained  below  and  the 
result  Is 

P(r)  -  1  -  Qtr)  ^ 

Q(r)  •  e"x  +  (ZxIC^e"*  J  dz  l}(2xz  K^)exp(-z2x  K)  ,  (11) 

3  o 

x  -  (r/rd)3  . 

where  Is  the  modified  Bessel  function  and  rd  Is  the  average  distance 
between  donors 

rd  -  (3/4imd)1/3  .  (12) 

Plots  of  Q(r)  for  various  values  of  K  are  given  in  Figure  1  and,  for  very 
large  values  of  r  one  can  easily  derive  the  asymptotic  expression 

Q(r)  *  (K16ir2z2)^  e"z  , 

(r/rd)3(l-KM2 


z  ■ 


(13) 


z 


3 


r/rd 


Fig.  1.  Q(r),  the  probability  of  not  being  able  to  enclose  an 

acceptor  in  a  charge  neutral  volume  of  radius  r,  vs.  r  in 
units  of  an  average  interdonor  spacing  rd.  Q(r)  is  plotted 
for  various  values  of  the  compensation  ratio  K  =  na/n  .. 
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For  small  values  of  r/r^  or  for  small  values  of  K,  P(r)  Is  dominated  by 
the  probability  of  enclosing  a  single  donor  In  the  volume.  However,  as 
K  increases.  It  Is  Increasingly  difficult  to  enclose  more  donors  than 
acceptors  and  thus  the  probability  of  large  positively  charged  volumes 
increases. 

On  one  hand  the  above  results  are  exact  and  do  yield  useful  Informa 
tlon  about  the  distribution  of  Coulomb  potentials.  On  the  other  hand 
they  do  not  give  any  quantitative  measurable  properties  of  the  system. 
Therefore,  In  order  to  make  quantitative  calculations,  we  assume  that 
at  low  temperatures  each  acceptor  can  be  associated  with  an  unfilled 
or  charged  donor  and  that  the  probability  density  p(r)  that  the  two 
are  separated  by  a  distance  r  Is  related  to  P(r)  by  the  equation 

4mr2  p(r)  ■  dP(r)/dr  .  (14) 

Further,  we  assume  that  a  filled  donor  cannot  be  closer  to  an  acceptor 
than  the  acceptor's  partner  Is  because  otherwise  the  acceptor  would  have 
paired  with  the  filled  donor.  Although  the  actual  situation  Is  far  more 
complex  than  our  simple  model  suggests,  our  model  does  correctly  reflect 
the  basic  physics  discussed  above. 

With  this  model  we  can  calculate  the  distribution  of  potentials 
In  the  sample  due  to  the  acceptors  and  charged  donors.  This 
Is  done  below.  The  distribution  function  itself  can  only  be 
expressed  In  terms  of  several  Integrals  so  we  have  computed  the  first 
two  moments  of  the  distribution.  These  results  are  plotted  In  Figure  2 
In  the  form 

<U>/Ed  •  -  (U,/Ed)  -  -ys  f,(K)  . 

<(U-<U>)J>%/Ed  «  (Uj/Ed)  •  ya  f2(K)  , 


(15) 


where  the  functions  and  f ^  are  defined  below. 

Further.  U-j  and  Ug  are  both  positive,  is  given  by  Eq.  (5),  ya  Is  a 
dimensionless  measure  of  the  average  spacing  between  acceptors, 

>,  *  V. 

r4  ■  (3/«ims)1/3  ,  (16) 

and  K  is  the  compensation  ratio.  The  quantity  -U^  Is  the  average 
potential  energy  of  an  electron  near  a  filled  donor  due  to  charged 
donors  and  acceptors  and  Ug  Is  the  rms  potential  energy  fluctuation 
about  -U-j  for  the  electron.  The  average  potential  energy  of  an 
electron  anywhere  In  the  sample  Is,  of  course,  zero  and  the  rms 
fluctuations  about  this  value  are  V?  U2*  Thus,  because  of  the 
spatially  non-random  distribution  of  charged  donors  and  acceptors, 
electrons  In  regions  of  the  sample  near  filled  donors  have  a  lower 
(more  negative)  average  potential  energy  with  smaller  rms  fluctuations 
than  electrons  at  arbitrary  positions. 


3.2  APPLICATIONS 


Based  on  the  Ideas  and  calculations  of  the  previous  section  we  can 
now  present  an  approximate  picture  of  a  low* temperature  n-type  III-V 
semiconductor  with  a  low  concentration  of  Impurities  and  compare  this 
picture  to  experimental  observations  on  GaAs.  When  discussing  shallow 
donors  we  shall  use  effective  mass  theory  and,  unless  discussing  transitions  * 
between  hydrcgenlc  states,  we  shall  Ignore  any  excited  states  of  the  donors, 
Including  other  hydrogenlc  states  Is  not  warranted  by  the  accuracy  of  our 
model . 

As  discussed  above,  the  potential  energy  of  an 

electron  due  to  charged  donors  and  acceptors  Is  zero  wlthanrms  fluctuation 
equal  to  Ji  Ug.  Because  of  these  fluctuations  the  bottom  of  the  conduction 
band  Is  not  a  well  defined  quantity.  However,  we  view  It  as  a  spatially 
varying  quantity  with  minimum  value  denoted  by  EQ, 

Ec  ■  a  U2  »  (17) 

where  a  Is  a  constant  of  order  one.  Electrons  that  are  bound  at  donor 
sites  feel  an  average  potential  energy  of  -U-|  wlthanrms  fluctuation  of 
Ug.  Thus  we  take  the  binding  energy  of  these  electrons  to  be 

Eb(x)  -  Ed  «■  U,  +  oxU2  (18) 

where  x  varies  between  plus  and  minus  one  and  E^  Is  the  effective  mass 
binding  energy.  Therefore  the  binding  energy  of  an  electron  with  respect 
to  the  minimum  of  the  conduction  band  Is  given  by  E(x)  where 

E(x)  •  Ed  +  U1  ♦aU2(-/?+x)  ,  -1<x<1  .  (19) 

That  Is,  the  energy  of  a  bound  electron  Is  an  amount  E(x)  below  the 
minimum  of  the  conduction  band. 
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Thus  we  predict  a  spread  or  distribution  of  donor  binding  energies 

with 

E„  *V“i  (20) 

*  <*+1>aU2 

where  E#v,  and  Em1n  are  the  maximum,  average,  and  minimum  binding 
energies  respectively.  At  temperatures  such  that  Ug  Is  of  order  kT  but 
not  much  less  than  kT  we  would  expect  E.u  to  control  the  temperature 
dependence  of  n  In  a  Hall  measurement  For  Ug  »  kT  we  would  expect  E^n 
I  to  be  the  controlling  factor  assuming  that  >  0.  If  E^  <  0  the 

whole  picture  breaks  down  although  there  may  be  a.  temperature  range  where 
the  analysis  is  partly  valid. 

Of  course  the  potential  energy  distribution  of  electrons  averaged 

over  the  whole  crystal  or  averaged  over  regions  near  filled  donors  does 

not  cut  off  sharply  In  our  model .  The  cut-off  at  energies  of  order  the 

ras  fluctuation  Is  taken  as  a  convenience.  Since  the  distribution  of 

potential  energy  Involves  simplifying  assumptions  anyway,  we  feel 

little  is  lost  by  taking  this  convenience.  However  we  note  that  our 

•  * 

distribution  predicts  an  exponentially  small  probability  for  large 
potential  energies  and  we  believe  that  the  exact  distribution  would 
also  have  this  property.  This  means  that  a  few  electrons  could  never 
be  bound  to  donors  although  they  could  be  trapped  In  regions  of  very 
small  potential  energy.  We  shall  not  pursue  this  fact  In  the  remainder 
-  of  this  report. 


We  now  consider  a  qualitative  and  quantitative  comparison  of  our 

theory  with  some  experimental  observations  on  GaAs.  First  we  consider  the 

dependence  of  the  donor  binding  energy  observed  In  Hall  measurements  on 

the  concentrations  of  donor  and  acceptors.  Table  I  contains  values  of 

Impurity  concentrations,  observed  donor  binding  energies,  and  computed 

values  of  K,  U-j ,  and  for  various  samples  of  high-purity  GaAs.  According 

to  Eq.  (20)  and  the  ensuing  discussion  the  observed  donor  binding  energy, 

Ej„,  should  be  E„„  which  Is  related  to  U.,  U0,  and  the  effective  mass 
ao  av  l  z 

donor  binding  energy,  Ed  by  the  equation 

CEd  *  U1 '  Edo>/Ed  ■ n  “  VEd  •  <21> 

In  order  to  check  this  and  to  determine  the  parameter  a  discussed  above 
we  have  plotted  (E^  +  U-j  -E^J/E^  vs  U2/Ed  *n  F*9“re  3  us1n9  the  values 
In  Table  I  and  Ed«5.8  meV. 

The  experimental  points  fall  on  a  remarkably  straight  line  with  a 
slope  corresponding  to  a  value  for  a  of  1.26.  The  fact  that  a  line 
through  the  points  does  not  pass  through  the  origin  Indicates  a  value 
for  E^  of  6.6  meV.  This  Is  somewhat  higher  than  the  5.8  meV  effective 
mass  energy  and  the  difference  Is  larger  than  any  expected  central  cell 
correction.  However,  In  view  of  the  straight  line  fit,  this  discrepancy 
appears  to  be  Independent  of  charged  Impurities.  Since  U-j ,  which  Is  not 
an  adjustable  parameter,  varies  between  20%  and  50%  of  E^-E^,  we  regard 
the  excellent  fit  as  evidence  that  our  model  has  a  reasonable  quantitative 
validity. 


Table  I 


The  quantities  na,  nd>  EdQ,  and  K  are  the  observed  acceptor 
concentrations,  donor  concentrations,  observed  donor  binding  energies, 
and  compensation  ratios.  The  quantities  U«|/Ed  and  were  computed 
using  Eqs.  (IS). 


Ref.  nax10“^2  cm3  ndx10’12cm3 


21.3  48.0 


40.7 

204 

136 

502 

327 

1060 

200 

860 

390 

490 

350 

610 

37.6  72.5 

104  489 


d01mnev 

K 

VEd 

VEd 

5.52 

0.44 

0.046 

0.103 

5.09 

0.20 

0.022 

mm 

4.51 

0.27 

0.045 

0.166 

3.88 

0.31 

0.069 

mm 

4.3 

0.23 

0.043 

0.182 

4.2 

0.80 

0.641 

0.414 

3.8 

0.57 

0.192 

0.293 

5.30 

0.52 

0.190 

0.132 

4.59 

0.21 

0.031 

0.143 

a)  G.  E.  Stillman,  C.  M.  Wolfe,  and  J.  0.  Dlmmock,  Proc.  3rd  Int.  Conf. 
Photocond.,  Stanford,  1969,  p.  265,  Pergamon,  Oxford,  1971. 

b)  Ref.  11. 

c)  G.  E.  Stillman  and  C.  M.  Wolfe,  private  communication. 
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Fig.  3 


Plot  of  (Ed  +  U1  - 
samples  in  Table  I. 


In  photoconductivity  experiments  on  high-purity  GaAs  at  low 

temperatures  one  observes  a  very  sharp  1s-2p  transition  whose  magnetic 

field  dependence  Is  In  excellent  agreement  with  effective  mass  theory 

g 

with  no  corrections  for  fluctuations.  On  the  other  hand,  the  Is-conductlon 
band  edge  Is  either  non-existent  or  very  diffuse.  Photoconduction  from 
the  1s-2p  transition  Is  believed  to  be  thermally  activated  from  the  2p 
level  to  the  conduction  band.  However,  the  activation  energy  of  the 
strength  of  the  transition  Is  less  than  either  EdQ/4  or  E^/4  as  would 
be  expected  from  effective  mass  theory.10 

Our  picture  Is  In  agreement  with  all  of  these  observations.  Since 
the  Coulomb  fluctuations  take  place  on  a  length  scale  much  greater  than 
aQ,  they  will  have  little  effect  on  the  low  lying  hydrogenlc  states  and 

o 

thus  £|S"E2p  k®  affected  as  Is  observed.  However,  the 

conduction  band  edge  Is  spatially  varying  on  a  length  much  greater  than 
aQ  and  thus  the  ls-conductlon  band  transition  will  have  a  spread  of  order 
2aU2  which  will  smear  It  considerably  as  Is  observed.  For  the  first 
sample  In  Table  I,  2aU2  is  1.5  meV.  This  Is  quite  close  to  the  observed 
smearing  for  this  sample.11  Further,  the  2p-conduct1on  band  energy  gap 
will  be  given  by  Eqs.  (19)  and  (20)  with  Ed  replaced  by  Ed/4.  Using  the 
same  analysis  as  that  leading  up  to  Eq.  (21)  we  obtain 

Eav(2s)  •  (Ed/4)  +  ut  -  v7oU2  .  (22) 

For  the  first  sample  In  Table  I  this  yields  an  activation  energy  of 
0.65  meV  which  Is  reasonably  close  to  the  measured  value10  of  0.47  meV 
and  much  closer  than  the  effective  mass  value  of  1.45  meV. 

Next,  we  examine  what  we  consider  a  most  puzzling  and  most 
underrated  piece  of  evidence  obtained  from  high-purity  GaAs  at  low 
temperatures.  It  Is  apparently  widely  believed  that  the  Holtzmark 


distribution  describing  the  distribution  of  electric  fields  from  a  random 

distribution  of  charged  Impurities  explains  the  lineshapes  observed  in 

photoconductivity  experiments.  However,  in  the  few  careful  comparisons 

between  theory  and  experiment* It  has  been  noted  that  while  the  lineshape 

predicted  by  the  theory  Is  fairly  good,  the  number  density  of  charged 

Impurities  necessary  for  quantitative  agreement  is  considerably  smaller 

7  8 

than  that  obtained  from  transport  measurements.  ’  This  is  particularly 
disturbing  since  the  observed  line  Is  narrower  than  the  predicted  line 
and  thus  additional  broadening  mechanisms  cannot  be  Invoked  to  explain 
the  discrepancy. 

Our  picture  provides  a  qualitative  and  even  a  semi quantitative 
explanation  of  this  behavior.  Although  the  charged  acceptors  do  constitute 
a  random  distribution,  the  charged  donors  do  not  because  they  are  more 
likely  to  lie  near  an  acceptor  than  one  would  deduce  from  a  random 
distribution.  Further,  the  filled  donors,  which  photoconductivity 
experiments  sample,  are  not  randomly  placed  either  but  are  likely  to 
be  further  from  an  acceptor  than  one  would  deduce  from  a  random  distribution. 
Both  of  these  effects  tend  to  decrease  the  electric  fields  at  filled  donor 
sites.  Preliminary  calculations  on  the  distribution  of  electric  fields 
have  been  performed  with  the  following  results.  As  K  approaches  one, 

*  <m 

the  distribution  of  electric  fields  Is  exactly  a  Holtzmark  distribution 
with  the  effective  number  of  Impurities  reduced  by  a  factor  of  four. 

For  values  of  K  nearly  one  the  distribution  of  electric  fields  Is 
nearly  Holtzmark  with  the  effective  number  of  impurities  reduced  by 
somewhat  more  than  four. 

Of  course  for  small  values  of  K  the  distribution  of  electric 

fields  will  more  closely  resemble  a  distribution  from  dipoles  with 

i/j 

a  density  n#  and  a  dipole  moment  of  rd  •  K  ra.  We  consider  these 


llneshape  effects  to  be  the  strongest  evidence  In  support  of  our  picture 
because  they  most  clearly  depend  on  the  non-randomness  of  the  distribution 
of  charges. 

Finally,  we  briefly  consider  the  mobility  of  samples  where  »  kT 
but  at  temperatures  high  enough  so  that  hopping  conduction  can  be  Ignored. 
At  temperature  kT^llg  and  below,  the  fluctuations  In  the  Coulomb  potential 
should  constitute  a  significant  scattering  mechanism  for  the  conduction 
electrons  that  Is  not  Included  In  the  potential  scattering  due  to  charged 
Impurities.  In  fact  when  kT  Is  considerably  less  than  Ug  there  must  be 
significant  volumes  of  the  sample  that  are  virtually  Inaccessible  to  the 
electrons.  Since  U,  depends  on  K  for  fixed  n  ,  this  mechanism  should  be 
more  pronounced  as  K  approaches  one.  In  fact  this  sort  of  behavior  Is 
observed,  for  example,  with  sample  P128(b)  and  P121(a)  of  ref.  11  which 
have  almost  Identical  values  of  n  .  However,  the  above  argument  only 
Indicates  that  our  theory  has  the  correct  trend.  Actual  calculations 
on  this  effect  would  be  very  useful. 


3.3  PROBABILITIES 
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In  this  section  we  sketch  the  solution  to  the  following  mathematical 
problem.  We  are  given  a  lattice  with  probabilities  c^  and  c^  of  having 
particles  of  type  1  (donors)  and  type  2  (acceptors)  respectively  on  a 
given  site  where 


c2  <  ci  «  1  *  (23) 

We  wish  to  find  Q(N),  the  probability  that  a  volume  of  N  lattice  sites 
does  not  contain  more  particles  of  type  1  than  of  type  2.  Thus, 


Q(N) 


nkj  N-m,  m_  N-m,. 

c^O-c,)  1  c2*(l-c2)  2 


(24) 


(n  \  u  m 
[m  J  ml  (N-m)  1  * 

The  summation  over  m2  can  be  expressed  In  terms  of  the  Incomplete  beta 
function  yielding 


N 


Q0(N)  »  (1-c,) 


m  N-m 

H,(»)  •  ej(7  -e, )  ■ 


(25) 


c2 

]2|  t-’o-t)"-"  dt  , 


sum 


m  >  1  . 

The  summation  over  m  can  be  performed  exactly  by  constructing  the 

N  N 


to  ♦*.”)(,  ♦«-,v  ■  £  t  (:)(* 

m^«0  n^O  '  l/\  2 


m.+m2  ie(m,-m2) 
x  £  e  1  £ 


(26) 


Integrating  8  from  0  to  2tt,  and  taking  a  derivative  to  obtain 


«2ir  oo 

7-  77  d0^  +  X2  +  2x  COS  0]N  ■  r 

dx  J  m.1 


N  m  v  2m 

m  mx 


Thus  Me  can  write 


o  o 

x  [1  +x2  +  2x  cos  0]N  , 


x  -  (Clt/(l-C|)(l  -t)r 


Since  c-|  and  c ^  are  much  less  than  one  we  can  expand 

[1 +x2  +  2xcos8]N  ■  exp(N  ln(l+x2  +  2x  cos  8))  *  exp(2Nxcos  0)  , 


the  expansion  being  valid  if  x*N  «  1.  One. can  easily  check  that  values 

14 

of  N  violating  this  restriction  would  constitute  10  atoms  for  concentra 

15-3 

tlons  of  Impurities  of  order  10  cm  .  For  N  this  large  Q(N)  Is 
ridiculously  small.  Further,  by  similar  arguments, 

(1  -  c)N  -  exp(-Nc)  , 

(1  -  t)N  =  exp(-Nt) 

x  a  (c<jt)^  .  (30) 

When  Eqs.  (28)  through  (30)  are  combined,  the  8  integration  can  be 
expressed  exactly  in  terms  of  a  Bessel  function  of  imaginary  argument  and 


one  obtains 


9 

*  e  [l  +  f  2  Tre'^lC.t)*  H  I,[2(C,t)*  Nil 


Equation  (11)  can  be  obtained  by  a  simple  change  of  variables  and  by 
noting  that 

(Nc,)  -  (r/rd)3  (32) 

for  a  sphere  of  radius  r  containing  N  lattice  points. 
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3.4  DISTRIBUTION  FUNCTION 

In  this  section  we  derive  the  distribution  function  and  first  two 
noments  for  the  Coulomb  potential  of  an  electron  due  to  the  non-random 
distributions  of  charged  donors  and  acceptors.  We  consider  a  random 
distribution  of  a  small  concentration  of  acceptors  where  U  (r)  Is  the 
potential  energy  of  an  electron  at  the  origin  due  to  an  acceptor  at  r. 
Further,  for  each  acceptor^  there  Is  a  donor  with  a  probability  density 
p(r)  that  the  donor  Is  at  a  position  r  away  from  the  acceptor.  The 
contribution  to  the  potential  energy  of  an  electron  at  the  origin  from 
a  donor  at  r  Is  Ud(r).  By  a  straightforward  generalization  of  a  method 
used  earlier,  one  can  easily  obtain  F(U)dU,  the  probability  that  an 
electron  at  the  origin  has  a  potential  energy  between  U  and  U  +  dU. 

F(U)  .  f  £  .,ut  .-«*>  . 

I(t)  *  nj,d3rd3r'  p(r-r'Hl -•xp[1t(U#(r)  +  U(J(r’))]J  ,  (33) 

where  n#  Is  the  density  of  acceptors. 

In  this  case, 

Ua(r)  -  -Ud(r)  -  e2/eQr  .  (34) 

For  an  electron  at  an  arbitrary  position  Eqs.  (33)  are  correct  as  they 
stand.  However,  for  an  electron  bound  (or  near  to)  a  filled  donor,  the 
acceptor  must  be  closer  to  Its  paired  donor  than  to  the  filled  donor  and 
thus  there  Is  a  restriction  on  the  r  and  r*  Integrations  such  that 

|r|  >  I r  -  r' |  . 


(35) 


The  moments  of  the  distribution  can  be  obtained  by  expanding  I(t)  In  a  power 
series  to  obtain 


For  k 


Mk- na(e2/e0)k  Jd3r  d3r'  p(r')(r-k-|r-rTk)  ,  r>r'.  (36) 

*  1,2,  one  of  the  Integrals  In  (36)  can  easily  be  performed  yielding 
M.j  ■  -(2imae2/3e0)J2  *  <u> 

Nj  *  (2naiie4/E^)J1  -  <(U  -  <U>)2> 

Jk  .  JVrrStf)  •  <37) 


By  using  Eq.  (14),  Integrating  by  parts  once,  and  changing  variables 

3 

from  r  to  N  ■  4irpr  /3  where  p  Is  the  density  of  lattice  points  we  obtain 


Jk  ■  (3/4irp)k/3  (k/3)fVk"3)/3Q(N)dN 


(38) 


h 


where  Q(N)  Is  In  Section  3.3.  By  using  Eq.  (31)  for  Q(N),  the  N  Integral 
Is  a  tabulated  Laplace  transform  of  a  Bessel  function  and  we  obtain 


Jk  -  (3/4irp)k/3(k/3)rr(k/3)c’k/3  +  r(2  +  k/3)^c^  J^dx 

L  1° 

P li/3^  ^  C2X)/ “  ^2x)  )j 


Aci-CgX)-1-"3  D'1 


(39) 


where  c1  and  c2  are  defined  In  Section  3.3,  r(z)  Is  the  gamma  function,  and 
Pk/3  Legendre  function.  By  combining  Eqs.  (39)  with  Eqs.  (37)  and 
making  another  change  of  variables  one  obtains  Eqs.  (15)  with 


f,(K)  ■ 

(2K2/3/3)r(2/3)[l  +  (10/9)L2(K)] 

(f2(k)>2- 

(2K1/3)r(l/3)[l  +  (4/9)L1(K)] 

Lk0()  ■ 

K  (x+l)(k“3)/3  F(-k/3,l+k/3;2,-x)dx 

(40) 

0 


where  F(a,b;c,z)  Is  the  hypergeometric  function.  The  Integral  Is  performed 
numerically  to  obtain  Figure  2.  . 
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4.  PIEZOELECTRIC  STRAIN  SCATTERING  FROM  NEUTRAL  IMPURITIES 


It  Is  well  known  that  static  strains  In  piezoelectric  materials  give 
rise  to  an  electric  field  or  an  electric  potential.  Thus  In  piezoelectric 
semiconductors,  such  as  I I I-V  compounds,  static  strains  can  have  an  effect 
on  the  electrical  properties  of  the  materials.  In  this  report  we  discuss 
some  of  these  electrical  effects  that  are  due  to  a  distribution  of  strains 
generated  by  a  random  distribution  of  point  defects. 

More  particularly,  we  calculate  the  mobility  due  to  a  concentration 

of  charge  neutral  defects  that  each  produce  a  given  strain  field.  The 

electric  potential  due  to  a  strain  producing  point  defect  In  a  piezoelectric 

crystal  Is  quite  similar  to  the  potential  due  to  a  point  dipole  and  both 

will  produce  a  mobility  that  Is  proportional  to  T1*  In  the  effective  mass 

approximation.  Except  for  scattering  by  Ionized  Impurities,  this  is  the 

only  mechanism  that  produces  a  mobility  that  decreases  as  the  temperature 

decreases.  Crude  estimates  of  the  strain  field  associated  with  a  single 

defect  are  estimated  from  llnewldth  measurements  on  a  SI  sample  with 

known  concentrations  of  specific  Impurities.  These  estimates  Indicate 

that  the  piezoelectric  static  strain  scattering  mechanism  may  contribute 

significantly  to  the  mobility  of  electrically  rather  pure  semiconductors 

below  room  temperature  when  neutral  Impurity  concentrations  are  greater 
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than  10  cm  .  Further,  it  could  be  a  dominant  meci.anism  in  determining 
the  mobility  of  I I I-V  semiconducting  alloys  In  some  regimes.  Other 
electrical  effects,  such  as  donor  lineshapes  as  measured  In  photoconductivity 
experiments,  are  also  qualitatively  discussed. 


4.1  POINT  DEFECT  MODEL 


According  to  elastic  continuum  theory^  a  defect  will  produce  a  lattice 
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displacement  proportional  to  r  and  a  strain  field  proportional  to  r  at 

distances  r  that  are  large  compared  to  the  size  of  the  defect.  In  this 

paper  we  shall  make  use  of  a  "model  defect"  defined  so  that  u(r),  the 

displacement  of  a  material  point  at  r  due  to  a  defect  at  the  origin,  is 

given  by  the  equation 

u(r)  »  b^/r3  ,  (1) 

where  b  has  the  units  of  length.  This  gives  rise  to  a  strain  field 

Uij(r)  *^[(3Ui/3rj)  +  (du^/dv^)]  -  (r25i^  «3rirj)b3/r5  .  (2) 

Several  cautionary  and  explanatory  remarks  are  appropriate  at  this 
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point.  The  model  defect  described  by  Eq.  (1)  is  a  fairly  common  model 
because  it  Is  rather  easy  to  manipulate  analytically  and  it  does,  of 
course,  possess  the  correct  long  range  behavior  which  determines  the 
dominant  features  of  our  results.  However  Eq.  (1)  Is  not  to  be  taken 
very  seriously  at  distances  within,  a  few  atomic  spacings  of  the  defect. 

In  fact  the  displacement  near  the  defect  is  a  very  difficult  problem 
that  is  largely  Irrelevant  to  the  present  problem  because  the  behavior 
of  u4<  at  distances  of  a  few  atomic  spacings  will  not  contribute  to  any 


expressions  in  this  report  Further,  although  Eq.  (1)  implies  a  dilatation 
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or  volume  change  of  order  b  ,  this  Is  incidental  to  this  report.  The  piezo 
electric  coupling  in  I I I-V  semiconductors  occurs  only  through  the  shear 
strain  components  u^  (i^j)  which  describes  a  trigonal  shape  change  of  a 
material  element  and  not  a  volume  change.  An  appreciable  shear  strain  in 
a  crystal  may  or  may  not  be  accompanied  by  any  volume  change.  Thus,  while 
one  expects  b  to  be  related  to  the  dimensions  of  a  defect,  it  cannot  be 
related  to  a  volume  change  or  a  nearest  neighbor  displacement.  In  this 
report  b  will  be  estimated  by  examining  experimental  determinations  of 
strain  fields  In  samples  with  known  defect  concentrations.  We  also  note 
that  Eq.  (1)  describes  a  spherically  symmetric  displacement.  This  is 
almost  surely  not  the  case  In  any  real  material.  Even  in  an  isotropic 
continuum  no  finite  number  of  force  pairs  will  produce  a  spherically 
symmetric  displacement  and  the  situation  is  worse  in  a  nonisotropic 
material.  In  general 

Uij(r)  »  f^CflJ/r3  (3) 

where  f^(n)  is  a  very  complicated  function  of  angles.  Thus,  any  detailed 
angular  dependence  predicted  by  our  model  strain  cannot  be  viewed  as 


4.2  DERIVATIONS 


In  this  section  we  derive  expressions  for  the  quantities  used  in  the 
rest  of  this  report.  First  we  obtain  an  expression  for  the  electric  potential, 
4,  due  to  the  model  defect  described  by  Eq.  (1).  We  use  the  notation  and 
basic  equations  from  ref.  (3)  which  includes  cgs  units  and  the  summation 
equation  for  repeated  Indices.  In  the  presence  of  a  strain  field  (and  to 
first  order  in  the  strains)  the  equation  connecting  the  electric  and  dis¬ 
placement  fields  is 

Di  ■  coEl  •  4lrei,Jkujk  •  (4) 

where  we  assume  an  isotropic  dielectric  constant  cq.  For  the  cubic 
sphalerite  structure  the  only  non-zero  element  of  the  piezoelectric 
coupling  constant  e^  is  e14  If  t,  j,  and  k  are  all  distinct  and  is 
zero  otherwise.  The  quantity  e-j4  has  the  units  of  charge/ length2.  The 
equations  »  0  and  ?  »  -v*  together  with  Eq.  (4)  yield 

V2*(r)  *  -4trp(r)  (5) 

where  p  Is  an  effective  charge  density, 

p(r)  «  (3/3xi)  Cei  ,jkujk^)/co^  ' 

Using  Eq.  (2)  for  our  model  defect  yields 

p(r)  »  (90b3e14/eQ) (xyz/r7)  ,  (7) 

where  spatial  directions  refer  to  the  crystalline  axes. 


Equations  (5)  and  (7)  can  be  solved  in  a  variety  of  different  ways  and 
we  have  found  Fourier  transforming  the  equations  to  be  most  convenient. 

Thus  all  functions  f(r)  have  a  Fourier  transform 


f(t)  •Jd3rf(r)e'i*‘r  , 

f(f>  -  f  d3|t,  fltle1*-*  , 

(2^)3 

(8) 

and  Eq.  (5)  can  be  written  as 

k^$(lc)  3  4ttp(1<)  .  (5') 

By  using  standard  methods  one  can  easily  find  that 

♦<t)  -  (96it21eub3/e<))(^kykz/k4)  . 

♦(f)  -  (36irb3e14/e0)  (xyz/r5)  .  (9) 

The  effects  of  Debye  screening  can  also  easily  be  Included  and  the  results 
are 

♦(f)  >  (96^ie14b3/c0)(kxkyy(k2(k2  +  k3))) 

♦(r)  -  (-24TO]4b3/30) (3/3x) (3/3y)(3/3z) (rp/r) (1  -exp(-r/rD)) ,  (9s) 

where  rQ  is  the  Debye  screening  length  and  kg  *  l/rD<  Thus,  <t>(r)  is 
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proportional  to  r  if  r  «  rQ  and  proportional  to  r  if  r  »  rQ. 


Given  the  electrical  potential,  the  calculation  for  the  elastic 

scattering  relaxation  rate  and  thus  the  mobility  is  perfectly  straightforward, 

at  least  in  the  Bom  approximation.  The  calculation  is  virtually  identical 
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to  the  one  described  by  Rode  for  ionized  impurities  except  that  the  Coulomb 
potential  is  replaced  by  the  piezoelectric  strain  field  potential  and  only 
one  point  needs  further  clarification.  Where  the  |$(lc)|2  enters  the 
calculation  we  replace  it  by  its  angular  average  analogue  |<t>(k)|  where 

|*(k)|2  *  <|*(£)|2>=  (3.210TT4/35)(e14b3/eo)2/k2  .  (10) 

for  the  unscreened  version.  This  simplifying  approximation  is  well  within 
the  spirit  of  neglecting  any  detailed  angular  dependence  discussed  above. 

The  relaxation  rate  in  the  effective  mass  approximation  is 

v  -  (3«29-tr3/5*7)(ee14b3/e0)2m*n/f»3k  .  (11) 

If this  Is  the  only  scattering  mechanism  under  consideration  then  the  mobility 
Is  easily  calculated  to  be 

y  a  (5-7/26-Tr3.32)(e01i/em*  e14b3)2(e/n)(2m*kT/TT)Ji  (12) 

and  the  Hall  factor  rH  is 

rH  *  r(7/2)r(5/2)/r2(3)  *  1.10  .  (13) 

Equations  (11)  through  (13)  were  obtained  by  ignoring  screening.  The  effects 
of  screening  can  easily  be  added,  if  necessary,  by  starting  with  Eq.  (9s) 
for  ♦{$)  Instead  of  Eq.  (9). 


In  order  to  make  contact  with  other  effects  from  the  strains  we 


consider  the  distribution  of  strains  due  to  a  random  distribution  of  model 
defects.  In  the  limit  where  the  number  of  defects  is  a  small  fraction  of 
the  number  of  lattice  sites  the  calculation  is  straightforward.  If  p(e)de 
is  the  probability  that  the  strain  at  a  given  point  is  between  e  and  e  +  de 
then 

p(e)  3  (e0/ir)[e2  +  e2]-1 
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eQ  3  (4ir/3)nb  ,  for  e  =  uxy 

eQ  3  (87r2/9/3]nb3  ,  for  e  =  uxx  .  (14) 

Finally  we  note  that  if  there  are  concentrations  ni  of  several  impurities 
with  associated  values  of  b^  then 


eQ  *  Zi  n,  b3 

v  'v-  p_1  ^  E1  ni  b® 


(15) 


These  equations  can  also  Include  the  effects  of  more  extended  defects  such 
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as  complexes  or  dislocation  loops  in  which  b  is  roughly  proportional  to 
the  volume  of  the  loop.  At  distances  far  from  the  defect  the  same  r 
dependence  will  be  obtained^  and  thus  the  k  dependence  of  v  and  the  T 


dependence  of  u  will  remain  the  same. 


4.3  EXPERIMENTAL  COMPARISONS 


One  major  result  of  Sec. 4. 2  Is  that  the  strains  generated  by  point 

defects  in  a  piezoelectric  semiconductor  produce  a  mobility  that  is 

proportional  to  T2*  and  a  Hall  factor  of  r^  *  1.10.  These  are  the  same  as 

would  have  been  produced  by  a  distribution  of  point  dipoles.  Except  for 

scattering  by  ionized  Impurities,  these  are  the  only  mechanisms  that  we 

know  of  that  decrease  the  mobility  as  the  temperature  decreases* 

The  obvious  question  is  whether  the  strength  of  the  static  strain 

piezoelectric  scattering  mechanism  is  large  enough  to  cause  a  measurable 
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effect  In  mobilities.  Concentrations  of  neutral  defects  of  order  10  cm 
are  common  in  almost  all  materials.  Further  EPR,  infrared,  and  optical 
measurements  on  defects  In  many  materials  yield  splittings  or  inhomogeneous 

line  broadening  of  about  1  cm"1.  Since  splittings  are  typically  of  order 

4  -1  -4 

10  cm  per  unit  strain,  this  implies  typical  random  strains  of  order  10 

Measurements  of  neutral  Impurity  concentrations  in  well  characterized  III-V 

semiconductors  are  rare.  There  are,  however,  reports  that  concentrations 

of  order  101®  cm"3  are  quite  common  even  in  electrically  rather  pure  GaAs.6 

However,  we  know  of  no  quantitative  analysis  of  strains  in  any  III-V 

semiconductors.  There  Is  one  quantitative  determination  of  the  strains 

due  to  oxygen  Impurities  In  Czochralski  grown  silicon  that  was  obtained  by 

their  effect  on  the  resonance  llneshape  of  deep  In  acceptors.  Mozurkewich, 

using  a  backward  wave  phonon  spectroscopy  technique,  measured  linewidths 

-4 

corresponding  to  strains  of  about  0.5x10  in  a  sample  with  an  oxygen 
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content  of  0.5x10  cm  .  Using  Eq.  (14)  this  implies  a  value  of 
b3^0.3xl0"22  cm-3. 


-50- 


8  From  these  numbers  we  estimate. that  n*10'8  cm"3  and  b3  =  0.3  x  TO"22  cm"3 
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are  at  least  not  unreasonable.  For  GaAs  with  eo=12.5,  e^4*4.7xl0  esu/cm  , 
and  m*  =  0.0665  m,  this  yields  a  mobility  of  1.5  xlO5  cm2/V  sec  at  a 
■  temperature  T  *  20°K.  This  is  lower  than  the  mobility  of  good  samples 

whose  mobilities  are  published.  However,  our  numerical  estimates  for  b 
are  uncertain  to  at  least  an  order  of  magnitude  which  will  lead  to  a  change 
in  the  mobility  of  a  factor  of  one  hundred.  We  feel  that  our  estimates  are 
certainly  no  better  than  that  but  that  they  do  show  that  the  mechanism  is 
worth  considering. 

We  have  included  the  possibility  of  static  strain  piezoelectric 
scattering  in  analyzing  the  mobility  of  one  well  studied  very  good  sample 
of  GaAs.  Using  reasonable  parameters  the  fit  could  be  made  slightly  better 
|3  than  without  the  mechanism.  However,  the  fit  could  also  be  improved  by 

changing  the  number  of  donors  and  acceptors  by  about  10%  and  thus  we  regard. 

» • 

this  attempt  as  Inconclusive.  Unfortunately,  there  is  probably  a  tendency 
|  for  data  on  well  characterized  samples  with  only  the  highest  mobilities  to 

reach  the  literature  and  these  samples  are  the  worst  candidates  for  the 
effect. 

One  might  expect  strain  effects  to  be  more  important  in  III-V  semi¬ 
conducting  alloys  than  In  the  pure  materials.  Our  analysis,  of  course, 
is  valid  only  for  rather  dilute  alloys  although  we  expect  the  qualitative 
features  to  be  present  at  all  concentrations.  We  would  also  expect  much 
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smaller  values  for  b  in  alloys  than  for  many  impurities  in  purer  substances. 

There  have  been  a  number  of  mobility  measurements  on  Ga^_xAlxAs  at 

i  temperatures  low  enough  so  that  static  strain  piezoelectric  scattering 
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might  be  detectable.  That  Is,  at  low  enough  temperatures  the  various 

a 

lattice  scattering  mechanisms  should  have  become  quite  ineffective  leaving 


i 


only  scattering  by  charged  impurities,  alloy  scattering,  and  possibly  static 

strain  piezoelectric  scattering.  One  expects  a  T  1  temperature  dependence 

in  the  mobility  for  what  Is  usually  called  alloy  scattering.  Although  the 
-k 

T  really  obtains  from  a  weak  scattering  limit,  it  is  difficult  to  see 
how  scattering  from  short  range  fluctuations  due  to  alloying  could  lead  to 
a  mobility  that  decreases  as  the  temperature  does.  That  is,  as  the 
temperature  is  decreased  the  conduction  electron's  average  momentum  is 
lowered  and  these  electrons  are  less  affected  by  spatially  small  potential 
variations. 

GaAs  and  AlAs  have  almost  Identical  lattice  parameters  and  thus  one 
might  expect  rather  small  effects  for  Ga^_xAlxAs.  However,  most  mobility 
measurements  do  exhibit  a  low  temperature  regime  where  the  mobility  decreases 
with  decreasing  temperature  and  with  increasing  alloy  concentration  x. 

This  Is  usually  Interpreted  as  an  Increasing  number  of  donors  and  acceptors 
as  x  Increases  even  though  n  (the  number  of  conduction  electrons  at  77°K  or 
at  room  temperature)  Is  not  correlated  with  the  composition.  We  suggest 
that  at  least  part  of  the  decrease  in  mobility  with  decreasing  temperature 
and  Increasing  x  may  be  due  to  static  strain  piezoelectric  scattering. 

In  order  to  conclusively  verify  this, Independent  determinations  of  NA  and 
Nq  would  have  to  be  made. 

We  have  crudely  analyzed  some  of  the  mobility  data  of  Chandra  and 
Eastman10  In  order  to  see  If  static  strain  piezoelectric  scattering  could 
be  a  dominant  mechanism.  The  analysis  was  limited  to  temperatures  between 
25° K  and  45°K  so  that  lattice  scattering  mechanisms  could  be  safely  ignored. 
Assuming  that  inverse  mobilities  add  and  that  the  number  of  charged 
Impurities  was  largely  independent  of  x,  one  can  obtain  concentration  or  x 
dependent  mobilities  by  subtracting  inverse  mobilities  of  different  samples. 


This  was  done  with  the  higher  concentration  sample  pairs  (F14.F16), 

(F14.F18),  (F16,F18)  and  (F15,F19).  The  results  yielded  mobilities  that 
varied  by  less  than  5%  over  the  temperature  range.  Either  a  T^  or  T"3* 
temperature  dependence  would  have  given  a  mobility  that  varied  by  more 
than  30%  over  this  range  although  the  right  combination  of  T3*  and  T"^ 
would  be  consistent  with  our  analysis.  Further,  all  pairs  gave  a  y-1 
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that  was  roughly  proportional  to  x  with  a  value  of  b  about  one  hundred 
times  smaller  than  discussed  earlier.  We  regard  this  analysis  as  suggestive 
but  certainly  not  conclusive. 

Finally  we  wish  to  make  a  few  qualitative  comments  on  the  effect  of 
strain  produced  electrostatic  potentials  on  the  donor  lineshapes  measured 
In  photoconductivity  experiments.  First  we  note  that  the  ls-2p-lineshapes 
measured  In  alloys  are  much  broader  and  more  symmetric  than  in  pure  compounds. 
The  typical  narrow  asymmetric  line  obtains  because  the  electric  field  from 
Impurities  separated  by  distances  much  greater  than  a  Bohr  radius  contribute 
only  to  second  order  In  the  ls-2p  energy  difference.  This  will  not  be  true 
for  the  electric  potential  due  to  strain  centers  that  are  separated  by 
distances  small  compared  to  the  Bohr  radius.  Secondly  we  note  that  different 
donors  are  characterized  by  different  llnewidths  (or  even  lineshapes)  In 
pure  I I I -V  semiconductors.  A  possible  explanation  of  this  Is  a  strain 
generated  potential  from  the  donor  defect  itself  j contributes  to  the 
llneshape.  These  Idea*  are  presently  being  pursued  quantitatively. 
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5.  Cr  COMPLEXES  AND  BULK  REDISTRIBUTION 
DURING  DIFFUSION  AND  ANNEALING 

Cr-doped  GaAs  has  been  studied  for  some  20  years.  The 
experimental  data  that  have  been  generated  over  this  period 
fall  mainly  into  three  categories:  spectroscopic  results, 
electrical  properties,  and  redistribution  phenomena.  Many 
properties  of  the  GaAs:Cr  system  are  still  not  well 
understood.  In  this  report  we  present  a  model  consistent 
with  a  wide  spectrum  of  experimental  data  from  the  three 
classes . 

We  start  by  reviewing  indications  of  Cr  interacting 
with  other  defects  in  GaAs,  and  evidence  of  the  existence 
of  highly  mobile  interstitial  Cr  donors.  We  can  then  give 
a  unifying  picture  of  the  effect  of  implantation  on  Cr,  and 
reinterpret  the  compensation  and  redistribution  mechanisms 
of  Cr  in  GaAs. 

5 . 1  DEFECT  INTERACTIONS 

It  is  becoming  increasingly  clear  that  one  cannot  in 


general  model  Cr  Independently  from  other  defects  that  are 
present.  We  will  therefore  start  with  a  short  survey  of 
experimental  data  indicating  defect  interaction  with  Cr. 


White  [la]  has  assigned  exci tonic  recombination  at  an 
isoelectronic  Cr-complex  (CrGaDAa)°  as  the  cause  of  the 
observed  0 . 84-eV  photoluminescence  line  earlier  attributed 
to  an  internal  Cr  transition  5E-5T2.  This  assignment  is 
supported  by  the  Zeeman  anisotropy  [lb]  and  by  the  angular 
dependence  of  the  optically  detected  magnetic  resonance 
[lc].  Picoli  et  al.  [Id]  proposed  an  interesting 
alternative:  an  interstitial  Cr  coupled  with  an  acceptor  A 
on  an  arsenic  site.  This  would,  according  to  the  authors, 
explain  the  0.575/0.535  eV  transition. 

Favennec  et  al.  [2]  observed  trapping  of  Cr  in  GaAs 
implanted  with  oxygen.  The  dose  was  rather  large,  which 
could  indicate  guttering  due  to  damage.  However, 
implantation  of  neon  under  the  same  conditions  (same  dose 
and  energy)  did  not  cause  a  build-up  of  Cr  in  the  implanted 
region.  The  authors  give  complexing  of  Cr  and  0  as  a 
possible  explanation.  Favennec  and  L’Haridon  [3,le] 
implanted  Se  and  Zn  into  GaAs:Cr.  Their  results  do  not 
indicate  any  complexing;  i.e.,  no  trapping  of  Cr  was 
observed  in  the  implanted  region.  Asbeck  et  al.  [4] 
implanted  Se  and  Kr.  No  trapping  of  Cr  in  either  case  was 
reported.  However,  Evans  et  al.  [5]  did  see  Cr  trapping  in 
the  Se  implanted  region.  They  used  a  dose  3  orders  of 
magnitude  larger  than  in  the  two  aforementioned  papers.  In 


the  latter  case  the  Cr  pile-up  disappeared  and  formed  again 
at  the  surface  after  further  annealing.  An  interesting 
case  was  reported  by  Magee  et  al.  [6].  They  annealed  GaAs 
implanted  with  high  doses  of  B  (isoelectronic  with  Ga). 
This  increased  the  outdif  fusion  of  Cr,  compared  to 
annealing  without  implantation,  and  resulted  in  a  stable  Cr 
depletion  channel,  A  small  tendency  for  Cr  pile-up  at  the 
implant  peak  was  observed  for  high  Cr  content.  Ne 
implantation  [2]  also  enhances  Cr  outdif fusion. 

Tuck  et  al.  [7]  observed  that  S  doping  of  an  epitaxial 
layer  grown  on  a  Cr  doped  substrate  decreases  the  Cr 
outdif fusion  more  than  an  undoped  layer.  This  indicates 
some  interaction  between  Cr  and  S. 

It  is  possible  that  Cr  is  sensitive  to  lattice  strain. 
Clegg  et  al.  [8]  point  out  that  Cr  doping  of  GaAs  reduces 
the  lattice  constant.  In  regions  of  large  concentrations  of 
As  vacancies  (e.g.  at  the  surface)  the  lattice  constant  is 
larger.  The  gettering  of  Cr  in  such  regions  may  decrease 
the  lattice  strain  energy.  Strain  due  to  interfaces 
between  encapsulants  and  GaAs  may  also  interact  with  Cr.  Eu 
et  al.  [9]  suggest  that  the  surface  pile-up  of  Cr  would  be 
caused  by  precipitation  of  Cr  at  dislocations  generated  at 
the  interface.  They  attribute  the  generally  smaller 
accumulation  of  Cr  at  unencapsulated  surfaces  to  smaller 
dislocation  densities.  This  contrasts  with  the  behavior  in 
the  bulk.  Simondet  et  al.  [If]  observe  that  Cr  prefers  to 


migrate  rather  than  to  precipitate.  They  perform 
implantation  of  Cr  and  study  the  redistribution  during  the 
post-implantation  encapsulated  anneal.  The  Cr  peak  flattens 
out  and  Cr  accumulates  at  the  interface  and  in  the 
encapsulant  (silicon  nitride).  A  higher  dose  of  Cr  causes 
more  damage,  and  if  Cr  had  a  tendency  to  precipitate  in 
damaged  regions,  the  Cr  peak  would  disappear  at  a  less 
rapid  rate,  if  at  all.  However,  at  higher  doses  the  Cr 
shows  the  same  tendency  to  leave  the  region  of 
implantation. 

Mullin  et  al.  [10]  observed  that  the  deep-level-forming 
species  in  their  samples  were  present  at  one  percent  or 
less  of  the  total  Cr  concentration  in  the  crystal.  This  led 
them  to  propose  that  Cr  gives  rise  indirectly  to  the 
species  that  controls  the  Fermi  level.  They  mentioned  a 
Cr-0  complex  as  one  possibility. 

Brozel  et  al.  [11]  studied  the  electrical  compensation 

mechanism  of  Cr  in  GaAs  doped  with  Cr  and  Si.  Using 

localized  vibrational  mode  absorption  measurements  they 
*  * 

concluded  that  there  were  no  signs  of  near-neighbor  Si-Cr 
pairs.  However,  they  did  observe  an  increase  of  the  C!r 
concentration  when  the  Si  concentration  was  increased. 
They  actually  observed  a  geometrical  adjustment  of  the  Cr 
concentration  to  the  Si  concentration.  This  behavior 


certainly  indicates  some  kind  of  interaction  between  Cr  and 
Si. 

Hobgood  et  al.  [12]  observed  that  the  post- implantation 
activation  of  Si  is  more  efficient  when  the  Cr  background 
is  higher.  Although  this  suggested  to  them  a  Cr-Si 
interaction,  it  does  not  seem  likely  that  this  is  due  to 
CrGa  forming  pairs  with  Si.  This  is  because  most  of  the  Si 
atoms  as  donors  reside  on  the  Ga  sublattice  where  Cr  is 
also  thought  to  reside.  This  would  indicate  that  a  Cr-Si 
pair  is  not  sis  favorable  as,  for  example,  a  Cr-0  pair. 

5 . 2  INTERSTITIAL  DONOR 

It  has  been  suggested  that  Cr  can  exist  as  an 
interstitial  donor  in  GaAs  [13].  In  a  paper  to  be  discussed 
in  detail  later.  Tuck  and  Adegboyega  [14]  model  their 
diffusion  data,  which  exhibit  exceedingly  rapid  diffusion 
of  Cr,  by  invoking  a  very  mobile  interstitial  Cr-species 
( Cr) .  We  will  assume  the  existence  of  very  mobile 
interstitial  Cr  donors .  There  have  been  experimental  data 
reported  which  support  this.  Deveaud  and  Favennec  [15] 
observed  a  new  PL  line  in  substrates  that  had  been 
Cr-implanted  and  annealed,  and  also  in  contaminated 
epitaxial  layers;  i.e.,  layers  into  which  Cr  had 
outdif fused  from  the  substrates .  A  very  interesting 
feature  of  the  latter  is  that,  although  the  part  of  the 


had  outdif fused  contained  almost  as 


layer  into  which  Cr 
much  Cr  as  the  substrate  itself,  it  was  not 
semi-insulating .  Furthermore,  the  new  PL  line  was  much 
stronger  than  those  usually  assigned  to  Cr.  The  authors 
suggested  that  Cr  may  be  located  non-substitutionally. 
Asbeck  et  al.  [4]  observed  a  spurious  n-type  layer  just 
below  the  GaAs-Si3N4  interface  after  encapsulated 
annealing,  and  in  this  region  Cr  is  always  observed  to 
accumulate  at  concentrations  above  the  expected  solubility 
[3,16].  Andre  and  Le  Due  [17]  actually  concluded  that  Cr 
incorporated  into  LPE  layers  form  shallow  donors,  contrary 
to  Cr  incorporated  during  bulk  growth.  As  mentioned  in 
section  2.1,  Picoli  et  al.  [Id]  argue  the  formation  of 
(Cr^A) -pairs.  With  our  assumption  this  seems  very  likely  to 
be  enhanced  by  an  electrostatic  interaction. 

5.3  ION  IMPLANTATION  EFFECTS 

The  spectroscopic  models  invoking  Cr  complexes  appear 
well  justified.  The  existence  of  these  complexes  would 
certainly  affect  both  redistribution  of  Cr  and  the 
electrical  properties  of  Cr  doped  GaAs .  Actually, 
successful  reinterpretation  of  both  the  redistribution  and 
compensation  mechanisms  in  terms  of  Cr  complexes  is 
possible,  as  we  will  see  in  the  following  sections.  In 
this  section  we  will  focus  on  the  cause  of  Cr 
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redistribution  about  implants,  and  present  a  unifying 
qualitative  interpretation. 

The  redistribution  of  implanted  Cr  discussed  above, 

indicates  that  Cr  is  much  more  prone  to  outdiffuse  and 

precipitate  at  the  interface,  than  to  be  gettered  in 

damaged  regions.  If  Cr  actually  does  accumulate  in  regions 

of  implantation  of  a  different  ion,  it  is  therefore 

presumably  not  because  Cr  is  attracted  by  the  damage,  but 

instead  due  to  a  direct  interaction  between  Cr  and  the 

implanted  species.  The  observed  redistribution  of  Cr  upon 

post-implantation  annealing  therefore  suggests  that  Cr 

tends  to  complex  with  the  column  VI  elements  (probably  most 

strongly  with  O),  but  that  the  barriers  to  be  overcome  in 

the  complexing  reaction  are  so  large  that  the  material  has 

to  be  nearly  amorphous  before  this  takes  place;  Simondet  et 
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al.  [If]  point  out  that  doses  of  about  10  cm  effectively 

disorder  the  material,  and  it  is  in  this  range  that 

trapping  is  observed.  The  dependence  on  the  dose  is  well 

illustrated  in  [lg]  where  implanted  S  traps  Cr  for  a  dose 
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larger  than  about  10  cm  .  The  proposed  mechanism  also 
conforms  well  with  the  depletion  of  Cr  in  regions  of  B 
implantation.  The  damage  will  undoubtedly  produce  a  large 
amount  of  interstitial  Cr,  which  will  have  to  compete  with 
B  for  Ga  sites.  Excess  interstitial  Cr  will  not  accumulate 
in  damaged  regions,  but  instead  outdiffuse  rapidly,  by  a 
mechanism  that  will  be  discussed  in  detail  in  section  6. 


In  summary,  Cr  does  not  appear  to  be  attracted,  or 
repelled,  by  bulk  damage  caused  by  implantation.  Any 
accumulation,  or  depletion,  of  Cr  in  implanted  regions  is 
likely  to  occur  by  an  attraction,  or  effectively  a 
repulsion,  by  the  specific  implanted  ions.  The  damage  only 
promotes  the  processes. 

5.4  COMPENSATION  MECHANISM 

Our  model  of  Cr  in  GaAs  relies  on  the  existence  of 
( Cr^D^ )  complexes  and  rapidly  diffusing  interstitial  Cr 
donors.  Even  though  the  (Cr^A)  complexes  were  proposed  as 
an  alternative  to  the  (CrGaDAs)  complexes,  it  appears 
likely  that  both  kinds  can  occur.  This  would  offer  an 
alternative  interpretation  of  the  compensation  in  Cr-doped 
GaAs.  As  was  first  observed  by  Cronin  and  Haisty  [18],  the 
semi-insulating  properties  of  bulk  grown  GaAs : Cr  are 
remarkably  independent  of  the  Cr  concentration  added  to  the 
melt.  This  behavior  is  typically  explained  [11]  by 
assuming  that  the  deep  acceptor  level(s)  of  CrGa  pin  the 
Fermi  level  close  to  the  center  of  the  gap.  Excess  Cr  can 
then  precipitate  [18]. 

We  suggest  the  following  compensation  mechanism:  At  the 
high  bulk-growth  temperatures  the  residual  donors  DAg  are 
chemically  compensated  by  Cr  through  the  formation  of  the 
isoelectronic  complex  (CrGaDAs)°.  The  excess  Cr  is  in  the 


form  of  isolated  Cr^  and  Cr^.  As  the  ingot  is  cooled  the 
Cr^  donors  flush  out  because  of  limited  solubility  and 
large  diffusion  constant.  Since  the  acceptors  and  the 
mobile  Cr^  donors  are  oppositely  charged,  electrostatic 
attraction  between  the  two  species  is  expected.  This  will 
enhance  the  formation  of  (CriA)  complexes  and  contribute  to 
good  compensation. 

In  addition  to  giving  an  alternative  interpretation  of 
the  semi-insulating  properties  of  GaAs  doped  with  Cr  during 
bulk-growth,  this  offers  an  explanation  of  the  enhanced 
apparent  activation  of  Si  in  the  presence  of  Cr.  Si  is 
amphoteric  and  the  Si-acceptors  can  be  neutralized  by 
rapidly  diffusing  Cr^ -donors.  The  same  effect  can  explain 
the  increase  of  Cr  concentration  with  increasing  Si 
concentration  that  was  observed  by  Brozel  et  al.  [11]  even 
locally. 

The  compensation  mechanism  gives  a  first  hint  to  an 
understanding  of  the  redistribution  properties  of  Cr  in 
GaAs  since  we  no  longer  can  expect  Cr  incorporated  during 
high-temperature  bulk  growth  to  redistribute  in  the  same 
way  as  Cr  introduced  into  the  crystal  at  lower 
temperatures.  The  rest  of  this  paper  will  be  devoted  to 
the  redistribution  problem. 


5.5  REDISTRIBUTION  MECHANISM 


The  foregoing  discussion  suggests  that  the  problem  of 
Cr  redistribution  in  GaAs  is  quite  complicated.  Cr 
apparently  occurs  in  several  different  forms:  isolated 
interstitial  Cr^,  isolated  substitutional  CrGa,  complexed 
interstitial  (Cr^A),  and  complexed  substitutional 
(CrGaDAs)‘  Furthermore,  there  are  a  number  of  mechanisms 
involved  in  the  redistribution:  diffusion  in  concentration 
gradients,  drift  of  Cr  ions  in  a  built-in  electric  field, 
quasi-chemical  reactions,  and  gettering  in  the  strain 
fields  of  surfaces  and  dislocations .  We  expect  the 
following  reactions  to  be  dominant: 


Cri  +  VGa:  ^Ga' 


Cri  +  VGa  +  °As  *  (CrGa°As) ' 


(1) 

(2) 


and 


Cri  +  A  t  (Cr^A) . 


(3) 


Three  assumptions  appear  physically  well  motivated: ' 

1.  The  actual  transport  of  Cr  occurs  interstitially.  The 
diffusion  and  drift  of  the  substitutional  and  complexed 
Cr  are  negligible.  This  is  typically  assumed  for 


interstitial-substitutional  diffusion. 

2.  At  the  usually  relatively  low  processing  temperatures 
there  is  no  formation  of  (CrGaDAs)  complexes  since  this 
is  a  third  order  process  and  seems  very  unlikely  to 
occur  after  the  original  formation  of  the  crystal. 

3.  The  complexes  are  considerably  more  strongly  bound  than 
the  substitutional  Cr  ions  to  their  sites. 

The  major  technological  problem  with  Cr  doped 
substrates  is  their  apparent  irreproducibility.  Capless 
anneal  sometimes  leads  to  surface  build-up  [lh,9],  at  other 
times  not  [lh,  19] ,  depending  on  temperature  and  As 
overpressure.  Capped  anneal  generally  leads  to  thin  surface 
pile-up  [3,9,16].  Apart  from  the  surface  behavior,  however, 
the  different  results  do  not  appear  too  dissimilar;  i.e., 
they  are  characterized  by  a  thin  (l-2um)  surface  depletion 
of  Cr,  while  the  bulk  concentration  remains  unaffected.  The 
really  astonishing  deviation  from  most  of  these  results  is 
reported  in  the  experiments  conducted  by  Tuck  and  coworkers 
[7,14].  They  observed  a  remarkable  penetration  and 
diffusion  of  Cr  in  GaAs  leading  to  a  uniform  Cr 
concentration  throughout  hundreds  of  microns,  and  a  thick 
(10-20 vim)  surface  pile-up  on  both  indiffusion  and  anneal. 
The  uniform  bulk  concentration  was  observed  to  vary 
exponentially  in  time.  In  this  section  we  will  focus  on 
this  apparent  inconsistency. 
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As  was  mentioned  in  the  last  section,  our  model 
indicates  that  the  redistribution  depends  strongly  on  how 
the  Cr  was  incorporated.  We  will  therefore  study  the  two 
cases  of  post-growth  incorporated  Cr  as  represented  by  the 
data  of  Tuck  and  coworkers,  and  Cr  incorporated  during  high 
temperature  bulk  growth  as  represented  by  the  data  of 
Kasahara  and  Watanabe.  The  latter  cure  typical  of  the 
behavior  usually  encountered.  The  results  lack  the  surface 
build-up,  presumably  due  to  stable  conditions  set  by  the  As 
overpressure  and  temperature.  In  both  cases  the  vacancies 
will  affect  the  redistribution,  and  we  will  therefore 
discuss  the  expected  vacancy  concentration  profiles  in 
GaAs.  In  the  analysis  we  will  ignore  the  drift  of  Cr  in 
electric  and  strain  fields.  This  means  that  we  do  not  model 
the  immediate  vicinity  of  the  surface.  However  we  will 
begin  by  giving  a  separate  qualitative  discussion  of  the 
drift  of  Cr£  in  electric  fields  due  to  dopant  gradients  in 
the  bulk. 

5.5.1  Interstitial  Drift 

Drift  in  built-in  electric  fields  can  typically  be 
neglected  except  at  surfaces  and  interfaces.  At  free 
surfaces  it  is  difficult  to  distinguish  between  field 
induced  and  strain  induced  build-up.  However  at  an 
interface  between  two  doping  regions  in  a  structurally 
homogeneous  crystal,  an  anomaly  in  the  Cr  distribution  is 
more  likely  to  be  caused  by  the  electric  field. 
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It  has  been  observed  [20]  that  Cr  outdif fusing  during 
epitaxy  from  a  semi-insulating  GaAs:Cr  substrate,  through 
an  undoped  GaAs  buffer  layer,  into  a  n-GaAs  active  layer 
(electron  concentration  n),  piles  up  at  the  buffer  layer 
side  of  the  interface  between  the  epitaxial  layers  for 
n=1017cm”3,  but  not  for  n=5xl016cm“3.  At  the  growth 
temperature  of  750°C  the  intrinsic  concentration  has  been 
measured  [21]  to  be  6xl016cm~3.  That  means  that  for  the 
case  where  pile-up  is  observed  there  is  an  electric  field 
in  the  interfacial  region  between  the  active  layer  and  the 
buffer  layer.  This  is  mostly  confined  to  the  undoped 
buffer  layer  and  is  directed  into  this.  The  interfacial 
pile-up  of  Cr  resembles  those  of  charged  impurities 
reported  in  [22].  These  were  caused  by  the  electric  field 
due  to  doping  gradients.  The  observed  Cr  pile-up  is  such 
that  the  mobile  Cr  would  have  to  be  positively  charged.  It 
is  therefore  likely  that  the  interstitial  Cr  actually  is  a 
donor  as  was  argued  above.  The  Cr  profiles  observed  by  Linh 
et  al.  [li]  in  MBE  layers  conform  qualitatively  with  this 
assignment. 

5.5.2  Vacancy  Profiles 

In  the  redistribution  model  we  have  outlined,  the 
vacancies  play  an  important  role  primarily  through  reaction 
(1).  One  possible  mechanism  for  Ga  vacancy  production  is 
the  process, 


. 


G*Ga  *  1  Z  VGa  +  Gai' 


(4) 


where  1  is  an  interstitial  site.  A  similar  reaction  would 
apply  to  As  vacancies.  It  is  also  conceivable  to  have 
vacancy  production  at  dislocations.  In  either  case  the 
continuity  equation  for  the  vacancy  concentration  Cv  can  be 
written  (neglecting  possible  drift), 

3c  32c 

*T  -  Dv^  -  *cv  +  *'  <5> 

where  Dv  is  the  diffusion  coefficient.  It  is  not  clear  what 
the  diffusion  mechanism  is,  but  it  seems  likely  that  it  is 
similar  to  that  suggested  by  Swalin  [23]  for  vacancy 
diffusion  in  group  IV  semiconductors.  The  .mechanism  would 
then,  involve  transfer  of  an  atom  to  a  nearby  vacancy  (for 
III  -V  semiconductors,  at  a  next-nearest  neighbor  site). 
a-6c  is  the  net  generation  of  vacancies  by  the  dominant 
mechanism.  At  the  surface  the  concentration  of  vacancies 
need  not  be  the  same  as  in  the  bulk  since  it  can  be  more 
readily  influenced,  for  instance,  by  the  As  overpressure. 
In  the  steady  state  the  vacancy  concentration  profile  is 

4o)(x)  -  C^b)  +  [C^a)-C^b)]exp(-x/Lv),  (6) 

where  C^a)is  the  surface  concentration,  C^b)  the  bulk 
concentration,  and  Lv  the  vacancy  diffusion  length  given  by 
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5.5.3  Radiotracer  Diffusion 

We  believe  that  Tuck,  Adegboyega  and  coworkers  [7,14] 
have  reported  some  of  the  few  experiments  that  involve 
essentially  only  reaction  ( 1 ) .  In  these  experiments 
radiotracer  Cr  was  diffused  into  bulk  grown  GaAs  (Cr-doped 
or  n-type )  at  temperatures  well  below  the  growth 
temperature .  Under  these  circumstances,  reaction  (2)  is  not 
expected  to  go  to  the  right,  and  reaction  (3)  is  secondary 
due  to  the  small  number  of  isolated  acceptors.  The 
indif fusing  Cr,  which  is  the  only  Cr  that  is  detected,  is 
thus  not  expected  to  interact  strongly  with  the  residual 
impurities . 

Tuck  and  Adegboyega  (TA)  account  for  the  observed 

exponential  time  dependence  of  the  uniform  bulk 

concentration  by  means  of  a  fairly  simple  model  [14] 

involving  the  non-equilibrium  chemistry  of  substititional 

and  interstitial  Cr  and  of  Ga  vacancies.  The  width  of  the 

surface  peak  permits  them  to  estimate  the  diffusion 
•  * 

constant  of  Ga  vacancies.  It  is  assumed  that  Cr  can  exist 
in  GaAs  as  a  highly  mobile  interstitial  species  Cr^  or 
substitutionally  on  Ga  sites  as  immobile  Crg.  The  reaction 
converting  each  species  into  the  other  is  reaction  (1). 
Only  Cr3  is  incorporated  or  desorbed  at  the  surface.  For  Cr 
to  penetrate  to  the  interior,  it  must  first  be  converted  by 


reaction  (1)  to  Cri.  It  then  diffuses  rapidly  and  can  again 
react  with  VGa  to  form  Crs.  The  limit  chose*  ?y  TA  to 
interpret  the  experiments  assumes  that  the  surface 
incorporation  and  conversion  reactions  are  in  equilibrium, 
that  the  interstitial  diffusion  is  extremely  rapid  (in 
effect,  instantaneous  on  the  time  scale  considered),  and 
that  the  bulk  conversion  is  slow.  In  addition  it  is  assumed 
that  no  vacancies  are  generated  or  annihilated  in  the  bulk 
except  through  reaction  (1).  The  surface  concentration 
gradients  are  qualitatively  accounted  for  by  the  diffusion 
of  VGa.  This  vacancy  concentration,  depleted  by  reaction 
(1)  during  indiffusion  of  Cr,  is  replaced  by  diffusion  from 
the  surface,  resulting  in  a  surface  Cr  peak,  &s  observed. 
In  the  outdiffusion  experiment  from  a  homogeneously 
vacancy-depleted  sample,  vacancies  also  diffuse  in  from  the 
surface,  "pulling"  Cr  from  the  interior  to  supply  reaction 
(1),  and  again  resulting  in  a  surface  peak. 

We  adopt  a  similar  model  but  with  a  substantially 
modified  picture  of  the  vacancy  dynamics.  Unlike  TA,  we 
assume  that  the  vacancy  concentration  is  always  able  to 
reach  its  steady  state,  on  the  time  scale  of  the 
experiments  under  study,  leading  to  the  vacancy 
concentration  profile  (6).  This  appears  to  us  more 
consistent  with  a  rather  large  vacancy  diffusion  constant, 
as  inferred  by  TA,  than  the  large  departure  from 
equilibrium  implied  by  their  interpretation.  We  could 
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indeed  adopt  TA's  viewpoint  and  generalize  it  to  . 

Cv(x,t)  ■  C^o) (x)-[Cs(x,t)-ca (x,0) ] ,  (8) 

where  we  assume  equilibrium  initial  distribution 
Cv(x, 0)«C^°* (x),  but  we  find  this  less  plausible.  Either 
approach,  however,  permits  us  to  extend  the  validity  of  the 
treatment  closer  to  the  surface,  and  we  are  able  to 
estimate  the  Ga  vacancy  diffusion  length  from  TA's  data.  We 
start  by  examining  rapid  diffusion  and  find  that  we  can 
deduce  a  condition  on  the  diffusion  constant  necessary  to 
keep  the  concentration  uniform.  From  this  we  can  estimate 
the  lower  limit  of  the  diffusion  constant  for  interstitial 
Cr,  again  using  TA's  data.  We  study  the  limit  considered  by 
TA  and  also  the  limit  where  the  surface  conversion  reaction 
is  the  bottleneck.  TA's  indiffusion  data,  which  exhibit  a 
large  constant  surface  concentration  on  indiffusion, 
indicate  incorporation  equilibrium,  so  this  must  be 
discarded  as  a  possible  bottleneck.  The  entire  process  of 
incorporation,  surface  conversion  of  crg  into  cr^ 
diffusion  of  Cr^  generation  of  Ga  vacancies,  and 
conversion  between  Cr .  and  Cr  in  the  bulk  is  illustrated 

1  9 

in  Figure  1. 

5.5. 3.1  Rapid  Diffusion  Leading  to  Uniform  Concentration 

We  study  an  atomic  species  that  redistributes  in  a 
sample  of  thickness  i  by  a  simple  diffusion  mechanism.  The 


with  boundary  conditions  set  by  the  surface  flux, 

J(0't)  *  ~D§§(0't)  "  v[kca-c(0,t)],  (10a) 

and 

JU,t)  -  -dH<a, t)  -  -v[kCa-C(A,t)].  (10b) 

In  these  equations  C(x,t)  and  Ca  are  the  concentrations  in 
the  sample  and  in  the  ambient  atmosphere  respectively,  D  is 
the  atomic  diffusion  constant,  k  the  segregation 
coefficient  and  v  the  surface  diffusion  velocity.  We  wish 
to  determine  the  condition  on  D  that  will  keep  the 
concentration  practically  uniform  at  all  times, 

C(x,t)  -  c* (t),  (11) 

assuming  the  consistent  initial  condition 

C(x,0)  -  C* (0) .  (12) 

With  these  initial  and  boundary  conditions,  the  solution  of 
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the  diffusion  equation  is  [24] 

C(x,t)-kCa  -  [C*  (0)-kCa]  [cosa^x+^^sina^x]exp  (-Da^t) , 
where  a1  is  the  first  positive  root  of 


(13) 


tanai.  -  4°Liy/D>  .  (14) 

o'4-  (v/D) 2 

For  C(x,t)  to  be  considered  independent  of  x  we  must  have 
<^*-<<1  and  v/o1d<<1.  If  these  conditions  are  fulfilled, 
(13)  becomes 


C’(t)  -  kC  +[C‘ (O)-kC  ]exp(-t/T), 


where 


1  .  ~2_  ..  d 
t  alD  J7 


and  =  2v/Dl  so  that  1/x  =  2v/l. 


(15) 


i 


i 

(16) 


Thus,  if  the  concentration  is  observed  to  relax 
uniformly  with  the  characteristic  time  t,  (16)  sets  a  lower 
limit  on  the  diffusion  constant, 


D  »  i2/ T. 


(17) 


5. 5. 3. 2  Bulk  Cr  Concentration. 


The  continuity  equation  for  the  Cr-.  concentration  Ci  is 
not  a a  simple  as  (9)  since  Cri  is  also  involved  in  reaction 
(1).  The  diffusion  of  Cr^  is  assumed  to  be  fast  enough  to 
maintain  a  uniform  concentration  and  much  faster  than  the 
rate  of  the  reaction.  Under  these  conditions,  the 
continuity  equation  becomes 


3ci 

Tt~ 


TE“' 


(18) 


where  the  bar  indicates  the  spacial  average  (l/i)fdx.  That 
is,  because  of  the  fast  diffusion  of  Cr^,  it  is  only  the 
global  effect  of  the  reaction  that  matters. 

Neglecting  diffusion  of  Crg,  the  continuity  equation 
for  the  Cr8  concentration  Cg  is  governed  entirely  by  the 
rate  of  the  reaction  (1), 


3Cs 

Jt~ 


cs 

Ttet* 

Tsi 


K&) 

+  7ETcv°)  {X)CJ 
Tsi 


It) 


(19) 


This  is  an  alternative  to  and  an  extention  of  TA's 
continuity  equation  where  we  adopt  the  equilibrium  vacancy 
picture  (6),  allow  for  the  time  dependence  of  the  Cri 
concentration,  and  superscribe  the  relaxation  time  xgi  and 
the  equilibrium  constant  <ig  with  "h"  for  ''bulk'' .  This  is 
done  to  distinguish  these  quantities  from  those  associated 
with  the  formally  identical  surface  conversion  reaction, 


for  which  the  parameters  may  be  different.  This  equation 
will  lead  to  an  ’'outer"  solution  [25],  which  does  not 
satisfy  the  boundary  condition  that  Cg  be  constant  at  the 
surface.  At  the  surface  there  is  a  boundary  layer  in  which 
the  outer  solution  connects  with  an  "inner"  solution  which 
falls  off  rapidly  due  to  finite  substitutional  diffusion 
and  effects  due  to  strain  and  fields .  This  will  be 
neglected  here  but  kept  in  mind  when  we  estimate  the 
vacancy  diffusion  length  later. 

The  boundary  condition  for  the  Cr^  concentration  is 
prescribed  by  the  flux  as  in  the  last  section.  This  flux 
equals  the  net  surface  production  of  Cri  by  the  surface 
conversion  reaction.  There  is  a  surface  generation  °9/r^ 
of  Cr  - ,  where  a  is  the  surface  concentration  of  Cr  and 

X  S  9 

T^®*is  the  relaxation  time  analogous  to  in  (19)  but 

9X  SX 

with  superscript  "s"  referring  to  the  surface,  o  is 
constant  since  we  have  assumed  incorporation  equilibrium. 
Thus, 


:sCada' 


(20) 


where  ks  is  the  segregation  constant  for  Crs,  Ca  is  the  Cr 
concentration  in  the  ambient  atmosphere,  and  da  is  the 
thickness  of  an  effective  surface  layer.  There  is  also 
surface  annihilation  of  Cr^;  i.e.,  Cri  becomes  Crs.  This 
process  is  described  by  a  similar  relaxation  term, 
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Ci(t)di/T^g)  where  d£  is  analogous  to  da  and  is  the 
relaxation  time  for  this  process.  We  find  then 


B 


B 


j(Oft)  -  -jUft)  -  v(kica“ci(t)]  ' 


(21) 


where  we  have  defined  the  surface  diffusion  velocity  v  and 
segregation  constant  ki  for  Cr^  as 


(22) 


and 


k,  *  k. 


da 

37  7737' 


si 


(23) 


to  point  out  the  similarity  with  the  boundary  conditions 
(10)  for  the  problem  in  the  last  section.  Any  direct 
incorporation  of  Cr  will  just  add  to  v  in  (22)  and  kA  in 
(23). 

Equation  (18)  in  its  present  form  is  not  similar  to 
equation  (9),  but  it  will  be  in  the  following  two  limiting 

cases : 

1.  Bulk  conversion  bottleneck:  If  conversion  from  CrA  to 

Jr~ 

Crg  in  the  bulk  is  so  slow  that  we  can  neglect  s  in 
(18),  it  becomes  identical  with  (9)  by  putting  D*Di. 
The  estimate  (17)  of  Di  is 


! 
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1  Di  >>  i2/^.  (24) 

However,  the  observed  time  constant  is  presumably 
which  is  associated  with  the  slow  bulk  conversion 
reaction.  Thus,  Ti<<T3i^  and  the  estimate  of  the 

required  magnitude  of  tends  to  get  larger.  On  the 
experimental  time  scale  then, 

ci<t)  -  kiCa  -  CiQ,  (25) 

and  the  solution  to  (19)  is 

Cs(x,t.>  -  x^,cioC^>  (x)  + 

+  [C>(x,0)-x^)CloC^o)  (xHexpI-t/T^1).  (26) 

The  total  Cr  concentration  is 

c(x,t)  -  cio+ca (x,t) .  (27) 

2.  Surface  conversion  bottleneck:  The  bulk  conversion  is 
assumed  to  be  in  equilibrium, 

Cs(x,t)  -  <ia)C^o) (x)ci(t) ,  (28) 


and  the  space  average  is 
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(29) 


Equation  (Id)  then  becomes, 

Dj 
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and  for  uniform  C 
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(30) 


(31) 


In  this  case  the  observed  relaxation  time  is  but  one 
usually  expects  <  to  be  much  larger  them  unity 
since  it  is  essentially  the  ratio  of  Cs  to  Ci.  So  again 
the  estimate  of  would  tend  to  get  larger  than  that 
associated  with  the  observed  time  constant.  The 
concentration  of  Cri  becomes 


C^t)  -  kiCa+(Ci(0)-kica]exp(-t/Ti), 


(32) 


and  the  total  Cr  concentration  is 


C<x,t)  -  cA(t)  H+<i3)C^o)  (x)] .  (33) 

5. 5.3. 3  Estimate  of  Ga  Vacancy  Diffusion  Length 

No  matter  which  of  the  conversion  reactions  is  the 
bottleneck,  a  plot  of  ln[C(x,t)-C(«,t)]  vs.  x  of  TA’s  data 


— ■  * — • 
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]  B  should  yield  a  straight  line  with  a  slope  -L^1.  If  the 

assumption  of  vacancy  equilibrium  is  correct  the  slopes 
should  be  independent  of  time  t.  The  data  for  such  a  plot 
!  0  from  TA's  Figures  1,  3,  4  and  6  are  displayed  in  our  Figure 

2.  From  this  it  is  apparent  that  the  data  points  for  x>10um 
can  be  fitted  with  straight  lines  except  possibly  for  the 
!  _  data  recorded  with  the  900°c,  4  hour  indiffusion  (line  g). 

The  750°C,  1  hour  outdiffusion  data  (line  i)  fall  on  a 

straight  line  all  the  way  to  the  surface.  This  is  not 

* 

„  surprising  since  there  is  no  large  Crs  surface 

concentration  producing  a  large  gradient  in  a  boundary 
layer.  All  the  lines  except  e  and  f  yield  roughly  the  same 
l  ft  value  of  Lv, 

Lv(Ga)  -  (11  t  3) inn.  (34) 

!& 

The  two  cases  e  and  f  are  responsible  for  effectively  the 
entire  standard  deviation  of  3ym.  It  may  be  argued  that 
this  value  of  Lv  is  not  precise.  However,  there  are  two 
important  points  to  be  made: 

1.  Lv  is  essentially  time  independent  (lines  a,b,c,d  and 
e ) ;  and 

2.  Lv  is  essentially  temperature  independent  (lines  b,f,g 
and  h) . 


The  time  independence  supports  the  assumption  of  vacancy 
equilibrium.  The  temperature  independence  is  not  unexpected 


C(oo,t)xcm 
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Figure  2.  Some  of  TA's  diffusion  data  replotted  in  order  to 
estimate  the  Ga  vacancy  diffusion  length.  The 
numbers  in  parenthesis  refer  to  the  figures  in 
[14]  used.  . 

a  i  au  e  ^  aaOa  t  1  /  A  \ 


a.  •  lOh,  110Q°C,  L»ll. 7ym  (4) 

b.  O  4h,  1100°C,  l3ll.7ya  (1,4) 

c.  B3h,  1100°C,  LXall.lym  (4) 

d. T2h,  1100°C,  &ll.lum  (4) 

e.  A lh,  1100°C,  rZ»14.7ym  (4) 

f .  □  4h,  100Q°C,  iZsS.lym  (1) 

g. V4h,  900TC,  l39.1ym  (1) 

h.  A  4h,  800rc,  &12.2ym  (1) 

i. Olh,  750°C,  L^a10.9ym  (6) 
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if  one  considers  the  vacancy  production  to  occur  through 
the  reaction  (4),  and  the  diffusion  mechanism  to  be  of  the 
Swalin  kind  as  discussed  earlier.  In  this  case  both  Dv  and 
6  in  expression  (7)  for  Lv  contain  an  exponential  of  the 
activation  energy  to  remove  an  atom  from  its  substitutional 
site.  If  an  additional  vacancy  production  mechanism  were 
present,  for  instance  due  to  a  high  dislocation  density, 
one  would  expect  Lv  to  decrease  (8  increase).  We  would  also 
expect  temperature  dependence  and  some  time  dependence 
since  the  properties  of  dislocations  can  change,  for 
instance  by  motion.  This  may  be  a  reason  for  the  anomalous 
cases. 

5. 5. 3. 4  Estimate  of  Interstitial  Cr  Diffusion  Constant 

During  indiffusion  we  have 

C(x,t)  *  C(x,«»)  ll-exp(-t/x)3 ,  (35) 

where  the  timeconstant  x  is  x^J*  if  the  bulk  conversion  is 

the  bottleneck  and  xi  if  the  surface  conversion  is  the 

bottleneck.  In  the  former  case  we  have  neglected  the 

concentration  of  Cr^  compared  to  the  equilibrium 

concentration  of  Cr_.  The  condition  on  D,  is 

s  i 

»  A2/t,  (36) 

Plotting  ln[l-C(x,t)/C(x,«) J  vs.  t  should  yield  a  straight 


•r*: 
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line  with  slope  -t”1.  in  Figure  3  we  show  such  a  plot  for 
x*~,  i.e. ,  in  the  uniform  bulk,  where  the  fit  is  expected 
to  be  best.  We  used  TA's  Figure  4  for  indiffusion  at 
1100°C.  The  estimated  straight  line  to  fit  the  data  points 
has  been  drawn,  keeping  in  mind  that  the  relative  errors 
increase  as  C(«o,t)  approaches  C(«,«).  It  yields  1*2.1x10  s. 
In  TA's  work  i*0.05  cm  so  A2/t  =10“7cm2/s.  A  rough  lower 
limit  for  would  be 

D± (1100°C)  >  10“Scm2/s.  (37) 

It  is  interesting  to  compare  this  with  the  largest 

conceivable  diffusion  constant.  This  occurs  when  the 
diffusion  is  driven  by  optical  phonons  with  a  jump 

frequency  equal  to  the  phonon  frequency.  For  interstitial 
diffusion  in  GaAs  the  optical  phonon  frequency  is 

fs»8xl012s“1.  The  interstitial  jump  distance  is  a/4, 

where  a  is  the  lattice  constant  (5.65  A),  and 

D^max)  -  a2f/8  *  3xl0“3cm2/s.  (38) 

•  * 

The  Cr±  diffusion  constant  is,  thus,  very  large  and  appears 
to  be  quite  close  to  the  phonon-driven  limit.  It  is 

interesting  at  this  point  to  recall  the  theory  developed  by 
Weiser  [26]  for  interstitial  diffusion  in  the  diamond 

lattice.  He  found  that  one  would  expect  some  ions  of 
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2 

t/h 


bulk  concentrat 
a  function  of 


stimate  the  associated  time 


intermediate  size  to  diffuse  interstitially,  virtually 
without  any  potential  barrier  between  sites.  He  proposed 
that  this  may  be  the  case  for  Cu  in  GaAs  which  would 
explain  its  diffusion  behavior. 

5.5.4  Bulk  Growth  Outdif fusion 

The  contrasting  experimental  results  that  we  set  out  to 
study  are  those  represented  by  the  data  of  Kasahara  and 
Watanabe  (KW)  [19].  These  authors  perform  annealing 
experiments  and  study  the  outdiffusion  of  Cr  from  the 
sample.  They  model  the  result  assuming  substitutional 
diffusion  and  a  finite  surface  diffusion  velocity.  They  get 
good  fit  far  from  the  surface  with  an  increasing 
discrepancy  between  theory  and  experiment  as  the  surface  is 
approached. 

The  present  model  offers  an  alternative  explanation 
which  appears  capable  of  accounting  for  the  results  closer 
to  the  surface.  In  the  experiments  the  measured  Cr  was 
incorporated  during  bulk  growth.  That  means  that  most  of 
the  Cr  is  presumably  bound  in  relatively  strong  complexes, 
mostly  as  ( CrGaDAa ) .  We  suggest  that  the  dominant 
redistribution  mechanism  is  the  dissociation  of  these 
complexes.  Note  that  there  is  no  uniform  reduction  of  the 
Cr  concentration.  If  the  dissociation  had  occurred 
throughout  the  sample,  one  would  expect,  in  view  of  the 
fast  diffusion  of  Cr^.  inferred  in  the  last  section,  that 
the  reduction  of  the  Cr  concentration  would  extend  further 


into  the  sample.  The  dissociation  appears  to  be  catalyzed 
close  to  the  surface.  The  defect  responsible  for  this  is 
likely  to  be  native  to  the  GaAs  crystal  since  the  surface 
depletion  appears  in  a  variety  of  samples.  The  region  of 
the  outdif fusion  is  much  thinner  than  the  Cr  peak  in  TA's 
data,  which  indicates  that  Ga  vacancies  are  not  directly 
involved.  We  therefore  suggest  that  the  arsenic  vacancies 
are  responsible  for  the  catalysis.  An  As  vacancy  adjacent 
to  a  (C7gaDas)  complex  could  very  well  increase  the 
probability  of  dissociation  which  would  then  be  described 
by 

+  VA.  :  Cri  +  +  VGa  +  VA«-  (39) 

By  a  technique  similar  to  the  one  used  to  analyse  TA's 
data,  we  can  estimate  the  As  vacancy  diffusion  length  from 
KW's  data.  We  assume  that  the  As  vacancy  equilibrium  is 
maintained.  This  appears  reasonable  since  no  VAS  are 
consumed  in  the  reaction  (with  the  possible  exception  of 
the  formation  of  di vacancies ) .  We  will  approximate  the 
equilibrium  vacancy  profile  (6)  by 

C^o)  -  <4#)exp(-x/Lv),  (40) 
since  the  surface  is  typically  very  VAg  rich.  In  view  of 
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the  rapid  diffusion  of  Crif  the  Cr^  formed  will 
"immediately"  result  in  a  uniform  distribution.  The 
concentration  of  will  be  exceedingly  small  and  we  can 
neglect  this  and  the  capture  by  VGa*  With  these  premises 
the  continuity  equation  for  the  total  Cr  concentration  will 
be 

■  -Y (x)c,  (41) 

where  y(x)  is  proportional  to  the  As  vacancy  concentration. 

With  a  uniform  initial  concentration  C  the  solution  is 

o 

C(x,t)  »  CQexp(-Y(x)t) ,  (42) 

with  approximation  (40)  for  the  As  vacancy  concentration. 

A  plot  of  ln[(lnC(x,t)-lnC0)/t]  vs.  x  should  yield  a 
straight  line  with  slope  -L^1.  For  different  times  the 
lines  should  coincide;  i.e. ,  the  position  and  slope  should 
be  time  independent.  We  plot  the  data  of  KW's  Figure  2  in 
the  described  manner  in  our  Figure  4 .  In  each  case 
(different  t)  a  straight  line  fits  the  points  very  well  and 
we  can  calculate  the  Average  estimated  diffusion  length  and 
the  standard  deviation, 

L  (As)  -  (0.5  ±  O.Durn.  (43) 


Figure  4.  KW’s  850 °C  anneal  data  ([19]  Figure  2)  replotted 
in  order  to  estimate  the  As  vacancy  diffusion 
length. 

a.  0  15  min,  L=0.38ym 

b. D30  min,  lZ=0 .45um 
C.O60  min,  lZ=0.50ym 
d.Al20  min,  l^-0 .  63  um 


The  standard  deviation  represents  a  reasonable  error.  The 
drift  of  slope  and  intercept  may  indicate  the  extent  to 
which  the  vacancy  equilibrium  fails  to  be  maintained. 

5.6  CONCLUSIONS 

We  have  proposed  a  model  for  Cr  impurities  in  GaAs  that 
is  consistent  with  a  wide  variety  of  experimental 
observations,  including  spectroscopic,  redistribution  and 
diffusion,  and  electrical  data.  This  model  does  not  appear 
to  be  in  conflict  with  any  reported  experiments.  In 
particular,  we  can  account  in  detail  for  the  excellent 
compensation  that  is  achieved  with  Cr  doping,  and  we  can 
reconcile  the  apparently  incompatible  but  equally  reliable 
diffusion  and  annealing  data  that  have  been  reported  by 
different  workers.  It  is  also  capable  of  interpreting  the 
observed  Cr  redistribution  about  implants  during  annealing. 

By  using  this  model  to  reinterpret  experimental  work 
reported  in  the  literature,  we  have  also  established  new 
estimates  of  several  useful  parameters  for  both  the 
impurity  and  the  host  crystal. 

We  hope  that  the  model  also  will  be  helpful  in 
analyzing  other  experimental  data  on  Cr  in  GaAs,  such  as 
thermal  conversion  of  GaAs:Cr  substrates  and  surface  Cr 
pile-up  during  heat  treatment,  and  those  associated  with  Cr 
outdif fusion  into  epitaxial  layers. 
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6.  Cr  SURFACE  REDISTRIBUTION  DURING  ANNEALING 


There  has  been  considerable  recent  interest  in  the  redistri¬ 
bution  of  Cr  in  GaAs  during  annealing  processes.  A  prominent 
feature  of  this  redistribution  is  often  an  accumulation  of  Cr 
at  the  surface.  This  pileup  has  been  observed  under  a  variety 

of  conditions,  and  attributed  to  As  vacancies1'”'*  or  encapsulant 
1  4-9 

strains  '  .  To  examine  these  possibilities  it  is  useful  to 

consider  the  characteristic  lengths  associated  with  Cr  buildup. 

In  many  cases  the  cause  of  a  physical  phenomenon  can  be  deter¬ 
mined  from  the  various  length  scales  associated  with  different 
mechanisms.  For  example,  the  estimated  diffusion  length  of 
As  vacancies  is  about  O.Sym1-®.  Although  the  characteristic 
lengths  associated  with  surface  strain  can  change  with  condi¬ 
tions,  it  is  also  a  long  range  phenomenon  compared  to  the 
smaller  regions  of  Cr  pileup  observed  at  the  surface.  In  this 
paper  we  examine  the  shortest  length  scale  associated  with  Cr 
accumulation,  which  is  apparently  due  to  the  electrostatic 
field  caused  by  surface  states. 

6 . 1  CHARACTERISTIC  LENGTHS 

The  characteristic  lengths  associated  with  surface  Cr 
buildup  as  determined  from  SIMS  data  in  recent  literaturi”^ ' 
and  from  our  own  annealing  experiments  are  plotted  in  Fig.  1 
as  a  function  of  annealing  temperature.  The  characteristic 
lengths,  L^ ,  were  extracted  from  the  SIMS  data  with  the  ex¬ 
pression1-^, 


N(r)  =  n.  coth2  +  a), 

i  Li 


where  N(r)  is  the  impurity  profile,  n^  is  the  intrinsic  carrier 
concentration,  r  is  the  distance  from  the  surface,  and 


Characteristic 


lltaratura 
our  data 
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Temperature  (K) 
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Comparison  of  characteristic  lengths,  L.  obtained 
from  recent  literature®”^ / 11-14  and  our  own  experi 
ments  with  the  intrinsic  Debye  length,  Lei.  Open 
symbols  indicate  encapsulated  samples  ana  closed 
symbols  indicate  unencapsulated  samples. 
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This  is  the  equilibrium  profile  expected16  for  positively 
charged  impurities  in  the  presence  of  an  electric  field  due 
to  surface  states.  Since  Cr  doped  substrates  are  nearly  in¬ 
trinsic  at  typical  annealing  temperatures17,  the  intrinsic 
Debye  length, LD, was  also  plotted  as  a  function  of  temperature  for 
comparison.  As  can  be  seen  in  Fig.  1,  all  the  data  points  lie  above 
or  near  the  intrinsic  Debye  length  curve.  This  suggested  to 
us  that  electric  fields  due  to  surface  states  is  the 
limiting  factor  in  the  redistribution  of  Cr  near  the  surface. 

6 . 2  THERMODYNAMIC  MODEL 


To  examine  this  possibility  in  more  detail,  an  equilibrium 
thermodynamic  model  was  developed  for  the  behavior  of  Cr  near 
the  surface.  Basically,  this  model  consists  of  a  fixed  surface 
charge  and  charged  Cr  -species  that  can  redistribute.  Space 
charge  neutrality  is  assumed  so  that  the  surface  charge  is 
exactly  compensated  by  the  net  charge  within  the  crystal.  We 
also  assume  a  surface  area  that  is  large  with  respect  to  the 
sample  thickness,  so  that  the  problem  is  one  dimensional.  The 
free  energy  of  the  system  is  then  minimized  to  produce  an 
equilibrium  distribution  of  Cr. 

We  first  consider  the  case  of  only  one  mobile  positively 
charged  Cr  ion.  The  equation  for  the  Helmholtz  free  energy, 

F,  of  the  system  is  given  by, 

F  -  -  j  V(r)N+(r)dr  +  1/2  J[N+(r)-N+]  v(r-r ' ) [N+ (r * ) -N+]  dr 'dr 
+  kT  (  N+ (r) lnN+ (r)dr  +  f  yN+(r)dr, 


(3) 
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where  V  is  the  potential  due  to  surface  states,  v  is  the  cou¬ 
lomb  potential,  y  is  the  chemical  potential,  r  is  the  distance 
from  the  surface,  N+(r)  is  the  concentration  of  positive  Cr 
ions,  and  N*  is  the  bulk  concentration  of  positive  Cr  ions. 

The  first  integral  on  the  right  hand  side  of  Eq.  (3)  is  the 
energy  due  to  surface  states,  while  the  second  term  corresponds 
to  the  coulomb  interaction  among  all  charges.  The  third 
integral  is  the  entropy  contribution  to  the  free  energy,  and 
the  last  term  is  the  chemical  potential  energy. 

To  determine  the  equilibrium  situation,  Eq.  (3)  for  the 
free  energy  was  minimized  everywhere.  This  resulted  in  a 
nonlinear  differential  equation  which  has  an  exact  solution. 


\  j 


_dy_ 


yV  y-l-ln  y 


where 


VT  In 


is  the  normalized  distance  from  the  surface. 


x.hVi 

No 


is  the  normalized  Cr  concentration,  and 


xQ  -  x (0) 


is  the  normalized  surface  Cr  concentration. 


(4) 


(5) 


(6) 


(7) 


The  free  energy  of  a  system  with  two  Cr  charged  species 
capable  of  redistributing  was  also  examined  and  is  given  by, 

F  »  -  V(r) [N+(r)  -  N- (r ) ]  dr 

•  * 

+  1/2  { [N+ (r)  -  N*]  - [N_ (r)-N~  ]  }v(r-r' ) { [N+(r ' )-N*  ] 

-  [N^fr')”  N~  ]}*dr'dr  +  kT  j[N+(r)ln  N+(r)  +  N~(r)ln  N~(r)]dr 
+  |  utN+(r)  +  N~ (r) ]  dr,  .  (8) 

where  N~(r)  is  the  concentration  of  negative  Cr  ions.  The 
terms  in  this  equation  have  the  same  meaning  as  in  Eq.  (3) . 

The  solution  for  Eq.  (8)  is. 


Equations  (4)  and  (9)  were  solved  numerically  with  the  results 
indicated  in  Fig.  2.  The  solution  for  two  Cr  ions  is  shown 
for  R  *  0.10.  As  can  be  seen,  in  both  cases  there  is  a  pileup 
of  N+(r)  at  the  surface.  The  profile  for  two  Cr  ions,  however, 
has  a  dip  in  total  concentration  due  to  the  depletion  of  N  (r) 
near  the  surface. 

6 . 3  ANNEALING  EXPERIMENTS 

To  determine  how  well  this  model  corresponds  to  an  experi¬ 
mental  situation,  unencapsulated  GaAs  substrates  were  annealed 
in  an  As  overpressure  using  a  boat  of  InAs  held  at  the  same 
temperature  upstream.  The  InAs  has  about  an  order  of  magnitude 


higher  dissociation  pressure  than  GaAs  and  has  been  previously 
used  to  stabilize  the  GaAs  surface  during  annealing.  SIMS 
analyses  done  on  two  of  these  samples  are  shown  in  Fig.  3.  One 
of  the  interesting  features  of  this  experimental  data  is  that 
the  characteristic  length  at  900K  is  much  larger  than  that  at 
1200K.  This  behavior  is  opposite  to  that  expected  for  As  vacan¬ 
cies  or  surface  strain,  but  in  agreement  with  the  temperature  de¬ 
pendence  of  the  intrinsic  Debye  length.  The  characteristic 
length  at  900K,  however,  is  about  a  factor  of  seven  too  large. 

Both  samples  show  a  dip  below  the  bulk  level  in  their  profiles 
which  is  in  qualitative  agreement  with  the  behavior  for  two  Cr 
ions  as  indicated  in  Fig.  2.  The  total  depletion  in  the  1200K  data, 
however,  is  too  large  to  be  explained  by  surface  states  alone. 

The  dashed  lines  in  Fig.  3  indicate  an  attempt  to  fit  Eq.  (9) 

to  the  experimental  data.  The  ratio  R  in  the  model  determines 

the  magnitude  of  the  dip  when  only  two  charged  Cr  species  are 

present.  If  there  are  neutral  Cr  atoms,  N°,  in  the  system, 

however,  they  will  tend  to  mask  the  dip  since  SIMS  analysis 

shows  the  total  Cr  concentration.  The  presence  of  neutral 

Cr  would  have  no  effect  in  Eq.  (9) ,  since  it  acts  only  as  an 

added  amount  of  Cr  distributed  evenly  over  the  sample.  It 

does  give  us,  however,  an  added  parameter  that  can  be  adjusted 

to  fit  the  experimental  data.  Also,  the  condition  of  charge 

neutrality  allows  us  to  determine  the  surface  charge  density, 

Ng.  The  fit  shown  in  Fig.  3  for  the  900K  SIMS  data  corres- 

-4 

ponds  to  an  R  of  5.8  x  10  ,  a  surface  charge  density,  N  , 

12  —2  0 
of  2.8  x  10  cm  ,  and  a  neutral  Cr  concentration,  N  ,  of 


101* 


1200K 


Fig.  3.  SIMS  data  for  900K  and  I200K  curve  fitted  to  the 
two  charged  species  model. 
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5.9  x  lO^cm  The  fit  for  the  1200K  SIMS  data  corresponds 

to  anRof  2.0  x  10-^,  an  N  of  5.5  x  10 ^cm"^ ,  and  an  of 

s 

7.4  x  10 The  values  of  Ng  obtained  in  this  process 
can  be  compared  to  a  previously  determined^  value  at  1000K 
of  5.1  x  lO^cm”^. 

As  shown  in  Fig.  3  the  agreement  between  the  model  and  the 
experimental  data  is  only  qualitative.  The  reference  point  used 
for  the  curve  fitting  process  was  located  at  the  point  where 
the  Cr  concentration  equaled  that  of  the  bulk.  The  calculated 
curve  obtained  by  this  method  indicates  that  the  surface  on 
both  samples  should  begin  further  in  the  material.  This  differ¬ 
ence  in  surface  location  could  be  due  to  the  presence  of  an 
oxide  layer  on  these  samples.  The  presence  of  surface  arti¬ 
facts  is  known  to  cause  matrix  effects  that  affect  the  accuracy  of 

19  20 

SIMS  analysis  near  the  surface  '  .  Part  of  the  cleaning  pro¬ 

cess  for  these  samples  involved  a  5:1:1  H2S04:H202 :H20  etch, 
which  invariably  leaves  a  thin  oxide  layer  on  GaAs  surfaces. 

To  support  this  conjecture,  the  data  indicate  a  thinner  oxide 
layer  for  the  sample  annealed  at  1200K  which  could  be  explained 
by  greater  oxide  reduction  in  a  H2  ambient  at  higher  temperatures. 

6.4  APPLIED  VOLTAGES 

To  investigate  further  the  role  that  electric  fields  play 
in  the  redistribution  of  Cr  at  the  GaAs  surface,  samples  were 
annealed  with  different  applied  voltages.  Due  to  the  surface 
states  present  there  is  a  built-in  potential  at  the  surface 
which  causes  band  bending.  By  applying  a  voltage  to  the  sur¬ 
face,  the  built-in  electric  field  can  be  altered  and  the  amount 
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of  band  bending  reduced  or  increased  accordingly.  If  Cr 
redistribution  is  electric  field  dependent,  then  larger 
electric  fields  and  band  bending  due  to  negative  bias  should 
cause  larger  positively-charged  Cr  pileups  at  the  surface. 
Conversely,  positive  applied  voltage  should  reduce  the  amount 
of  Cr  pileup. 

For  this  purpose  an  MOS  structure  was  used  to  apply  various 
voltages  to  different  areas  of  the  substrate  during 
annealing.  Anneals  were  performed  from  800  to  1000K  for  two 
different  time  periods  at  each  temperature.  The  SIMS  data 
shown  in  Fig.  4  are  typical  of  the  results  obtained.  As  can 
be  seen,  a  negative  bias  on  the  MOS  structure  produced  a  larger 
build  up  of  Cr  near  the  surface,  while  a  positive  bias  reduced 
the  build  up  of  Cr.  This  behavior  is  in  agreement  with  the 
expected  behavior  for  positively-charged  Cr.  In  all  the  samples, 
however,  the  Cr  depletion  beyond  the  surface  is  larger  than  can 
be  attributed  to  the  effect  of  surface  states. 

In  conclusion,  our  results  indicate  that  the  electric 
field  caused  by  surface  states  is  a  limiting  factor  in  the 
redistribution  of  Cr  at  the  GaAs  surface.  This  is  supported 
by  the  following  observations:  (1)  The  lower  limit  on  the  charac 
teristic  lengths  of  the  surface  regions  is  very  close  to  the  in¬ 
trinsic  Debye  length.  (2)  A  thermodynamic  model  for  the 
redistribution  of  Cr  in  the  presence  of  surface  states  predicts 
the  general  shape  of  the  Cr  profile.  (3)  The  application  of 
electric  fields  to  the  surfaces  during  annealing  changes  the 
magnitude  of  the  Cr  build  up. 


Fig.  4  .  Cr  redistribution  as  a  function  of  applied  voltage. 
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7.  CHARGED  IMPURITY  REDISTRIBUTION  DURING 
EPITAXIAL  GROWTH 


The  incorporation  of  impurities  in  epitaxially  grown 
layers  of  semiconductors  and  their  subsequent 
redistribution  often  leads  to  impurity  and  and  carrier 
density  profiles  that  are  unforeseen  and  that  may  be 
unfavorable  for  certain  applications.  Examples  are:  the 
formation  of  high  resistance  or  inverted  polarity  layers 
in  Gunn  devices  and  the  occurence  of  high  conductivity 
layers  close  to  the  semi- insulating  substrate  in  FET's  and 
IC's. 

A  model  has  been  proposed  by  Wolfe  and  Nichols  [1,2] 
to  account  for  some  of  these  effects.  In  this  model  the 
impurities  redistribute  by  diffusion  in  concentration 
gradients,  drift  in  the  built-in  electric  field  due  to 
surface  states  and  the  doping  discontinuity  at  the 
layer-substrate  interface,  and  finally  due  to  the  growth 
itself.  In  Figure  l  the  geometry  of  simple  planar 
epitaxial  growth  is  shown.  The  incorporation  and 
redistribution  of  impurities  are  indicated. 

In  order  to  model  these  processes  analytically  it  was 
necessary  to  make  several  restrictive  assumptions,  some  of 
which  were  known  to  be  unrealistic. 


v(t) 


INCORPORATION  OP  IMPURITIES  AT  MOVING 
6R0WTH  SURFACE 


Figure  1  'Geometry  of  an  epitaxial  layer  growing  with  a 
velocity  v(t)  on  a  thick  substrate. 


In  particular,  it  was  assumed  that  the  impurity 
distributions  were  stationary,  that  no  impurity  transport 
crossed  the  substrate-epitaxial  layer  interface,  and 
calculation  was  restricted  to  epitaxial  layer  impurity 
concentrations  so  small  as  not  to  affect  the  interfacial 
electric  fields.  Therefore,  although  insight  could  be 


gained,  especially  into  the  usually  neglected  effect  of 
the  field,  no  quantitative  conclusions  could  be  drawn.  In 
the  present  work  we  improve  the  model,  renounce  the  less 
realistic  assumptions  and  provide  for  the  simulation  of  a 
greater  variety  of  growth  conditions. 

7.1  NONEQUILIBRIUM  MODEL 

With  typical  ionic  diffusion  constants  and  mobilities, 
and  with  epitaxial  growth  rates  of  technological 
importance,  the  impurity  distribution  does  not  reach  a 
steady  state  during  epitaxy.  Transport  of  impurities 
across  the  epitaxial  layer-substrate  interface  can  be  an 
important,  even  dominant  feature  of  impurity 
redistribution.  If  the  doping  of  the  epitaxial  layer 
exceeds  the  intrinsic  carrier  concentration, 
redistribution  affects  the  field.  These  points  are 
accommodated  in  the  new  model.  Moreover,  although  the 
specific  simulations  we  have  carried  out  assume  constant 
initial  and  intentional  doping  and  uniform  growth  rate, 
the  model  is  designed  to  accommodate  arbitrary 
nonuniformities  in  both  space  and  time. 
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7.1.1  Assumptions 

The  following  assumptions  are  made: 

1.  The  material  is  structurally  homogeneous 
(homoepitaxy) . 

2.  The  material  is  nondegenerate,  and  the  electrons  and 
holes  sure  in  thermal  equilibrium.  This  implies  that 
the  electron  and  hole  concentrations  (n  and  p 
respectively)  satisfy 

np  -  nj,  ( 1) 

where  nA  is  the  intrinsic  carrier  concentration. 
From  the  absence  of  carrier  current,  the  electron 
concentration  and  electric  field  E  are  related  by 

E  -  |jln(n) .  (2) 

3.  The  substrate  is  effectively  semi-infinite,  and  the 
potential  and  field  deep  in  the  substrate  go  to  zero. 

4.  The  growth  surface  is  characterized  by  a  fixed  surface 
state  concentration.  These  states  are  always  fully 
ionized  and  the  field  at  the  surface  (x»0)  is 
constant,  and  determined  by  the  surface  charge  density 

V 

eE(0,t)  -  aM,  (3) 
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where  e  is  the  total  dielectric  constant. 

5.  The  ambient  phase  is  ideally  stirred,  so  that  there  is 
no  diffusive  impurity  flow  there. 

6.  There  are  no  temperature  gradients. 

7.  All  impurities  jure  fully  ionized.  This  means  that  the 
total  doping  Q  is 


Q  'E  w  l4) 

k 

where  Nk  is  the  concentration  of  the  impurity  Ik, 
and  zk  the  number  of  elementary  charges  it  carries. 
zk  is  >1  for  a  simple  donor  and  -1  for  a  simple 
acceptor. 

8.  The  impurities  diffuse  by  a  simple  non-reactive 
mechanism.  This  implies  that  the  flow  Jk  of  Ik  is 

Jk  -  "W* +  “k™*-  (5) 


where  the  ionic  mobility  for  Ik  is  given  by  the 
Einstein  relation; 


zkq 

VeJV 


(6) 


Bk  is  a  numerical  factor  which  depends  on  the 
diffusion  mechanism;  it  is  1  for  interstitial 
diffusion  and  1.27  for  a  vacancy  mechanism  [3].  The 


continuity  equation  for  1^  is 


3Nk  3Jk  (7) 

3t  3x  ' 

without  generation  or  recombination  terms. 

9.  The  surface  conductance,  or  diffusion  velocity,  which 
determines  the  outdif fusion  from  the  growth  surface 
[4],  is  much  smaller  than  the  growth  velocity.  The 
implications  of  this  assumption  will  be  discussed  in 
the  next  section. 

10.  The  high  temperature  impurity  distribution  is  frozen 
in  upon  cooling. 

7.1.2  Incorporation  of  Impurities  at  the  Growth  Surface 
We  will  adopt  a  simple  kinetic  model  for  the 
incorporation  of  impurities.  The  values  of  the  different 
parameters  will  depend  on  the  detailed  chemistry  of  the 
incorporation . 

At  the  interface  between  the  semiconductor  and  the 
ambient  phase  (a  gas  or  a  liquid)  there  is  a  dynamical 
process  in  which  atoms  leave  the  solid  and  go  into  the 
ambient,  and  vice  versa.  Here  we  are  interested  in 
impurity  atoms.  We  will  describe  the  situation  by  two 
opposite  flows  as  illustrated  in  Figure  2.  N(x,t)  is 
the  concentration  of  the  impurity  in  the  solid,  and  N  (t) 
the  concentration  of  the  impurity  in  the  ambient  phase 
(for  the'  present  purpose  we  leave  out  the  subscript  for 


Figure  2  incorporation  of  impurities  modeled  by  two 
opposite  flows  at  the  solid-ambient  interface. 
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the  impurity  type).  N_  is  assumed  to  be  uniform 

(assumption  5)  due  to  good  stirring.  This  is  really  an 
equivalent  concentration,  since  impurity  atoms  can  be  tied 
up  in  molecules  in  the  ambient  phase.  The  chemistry  of 
the  incorporation  is  compounded  in  the  flow  v  N  (t).  The 

3S  3 

opposite  flow  is  given  by  v  N(0+,t).  v  and  v  are 

surface  diffusion  velocities,  related  by  the  distribution 
coefficient  k;  v  =kv  .  The  net  flow  just  outside  the 
solid  is 

j(0-,t)  -  v  N  (t)  -  v  N (0+, t)  •  (8) 

AS  CL  S3 

Just  inside  the  solid  the  flow  is  due  to  diffusion  and 
drift; 

J(0+,t)  -  -l>||(0+,t)  +  yE(0+,t)N(0+,t) .  (9) 

If  the  boundary  is  moving  with  a  velocity  v  in  the  -x 
direction  due  to  growth,  and  we  fix  our  frame  of  reference 
in  the  moving  boundary,  we  have  to  add  a  term  vN  (t)  to 
the  right  hand  side  of  (8)  and  a  term  vN(0+,t)  to  the 
right  hand  side  of  (9).  The  current  must  be  continous 
at  x*0,  so  the  boundary  condition  becomes 

-Df5(°»t)  +  (uE(0,t)+v(t)+V„]N(0,t)  - 

ax  sa  (10) 

-  [v(t)+vaa]Na(t). 


» 
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If  there  is  no  diffusion  or  drift  in  the  solid  either, 
v(t)*v,  N(x,  t)=Ng  and  Na(t)=Na,  independent  of  time  and 
position,  and  (2.10)  simplifies  to 


vas+v 

vsa+v 


(11) 


This  is  one  of  the  expressions  for  the  effective 
distribution  coefficient  keff  reviewed  by  Kroger  [5]  in 
connection  with  purification.  In  the  case  of  sufficiently 
fast  growth  (assumption  9),  we  can  neglect  v„_  and  vs  . 
and  keff  approaches  1. 

7.1.3  Critique  of  the  Assumptions 

For  some  of  the  assumptions  the  validity  is  determined 
by  how  the  growth  situation  is  set  up;  for  others  by  what 
types  of  impurities  are  used.  Most  electronically  useful 
materials  are  nondegenerate  especially  at  growth 
temperatures.  Also,  many  impurities  are  fully  ionized, 
even  at  room  temperature.  During  growth  this  is  true  for 
an  even  larger  class  of  impurities.  If  the  chemistry  of 
any  reactive  diffusion  mechanism  present  were  well 
understood,  the  model  could  be  extended  to  include  this  as 
will  be  discussed  below. 

The  assumption  of  constant  surface  field  is  an 
approximation.  In  principle  one  should  fix  neither  the 
potential  nor  the  field  at  at  the  surface,  but  instead 
seek  a  self-consistent  solution  taking  into  account  the 


- 


m 


energy  dependence  of  the  surface  state  density,  a  quantity 
about  which  little  is  known  at  growth  temperature.  We 
shall  see  below  that  for  common  growth  conditions  the 
choice  of  boundary  condition  made  here  is  not  critical. 

There*  is  very  little  information  on  the  surface 
conductance  (or  diffusion  velocity)  of  impurities  in  GaAs. 
It  is  known  to  be  very  small  for  Sb  diffusing  in  Ge  [4], 
and  a  tentative  analysis  of  more  recent  data  on  the 
outdif fusion  of  Cr  from  GaAs  [6]  yields  a  value  of  lum/h, 
consistent  with  our  assumption  under  typical  growth 
conditions . 

In  most  instances  the  combination  of  diffusion 
constants,  their  temperature  dependence,  and  rapid  cooling 
rate  will  make  the  last  assumption  a  very  good  one. 

In  summary,  the  assumptions  made  here  are  not  very 
restrictive,  and  the  model  should  be  useful  to  simulate  a 
wide  variety  of  growth  processes.  It  is  to  be  understood, 
as  indicated  in  Figure  l»  that  we  model  only  planar 
growth.  Lateral  growth  and  edge  effects  are  excluded  from 
consideration. 

7.1.4  Governing  Equations 

We  make  the  following  normalizations: 

1.  The  position  x  is  expressed  in  terms  of  the  intrinsic 


Debye  length; 
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2.  The  potential  V  is  expressed  in  terns  of  the  thermal 
voltage ; 


V  *  — 
VT  q  * 


3.  The  electric  field  E  is  expressed  in  terms  of 


(  13) 


E„ 


T 


(  14) 


4.  The  concentrations  n,  p,  N  and  Q  are  expressed  in 

terms  of  the  intrinsic  carrier  concentration  n . . 

x 

5.  The  time  t  is  expressed  in  some  convenient  unit  tQ. 

6.  A  diffusion  constant  Dk  is  expressed  in  terms  of 


4 


o  t. 


7.  The  growth  velocity  v  is  expressed  in  terms  of 


V°  "  ^i* 

8.  An  impurity  flow  Jk  is  expressed  in  terms  of 


Jo  "  niV 


(  15) 


(  16) 


(  17) 


9.  The  conductivity  o  associated  with  the  mobile  carriers 


is  expressed  in  terms  of 


The  electrostatic  potential  is  obtained  by  solving 
Poisson's  equation,  which  for  a  one  dimensional 
nondegenerate  semiconductor  with  the  mobile  carriers  in 
thermal  equilibrium  reads  [7] 


(19) 


This  will  be  referred  to  as  the  Shockley-Poisson  equation. 
It  has  been  written  with  partial  derivatives  although  it 
is,  strictly,  only  valid  for  a  static  situation.  However, 
the  change  in  time  of  the  total  net  doping  Q  is  so  slow 
that  an  electrostatic  approach  is  valid.  The  electric 
field  is 


E  - 


(20) 


The  boundary  conditions  for  our  semi-infinite  case  are 


The  constant  is  given  by 


*0  -  sinh"1(Q0/2)/  (22) 

where  Q0aQ(°°,t),  independent  of  time. 

From  the  potential  the  electron  and  hole 
concentrations  can  be  obtained  from  the  relations 

n  -  n(«)exp(i|i) ;  p  -  1/n.  (23) 

The  conductivity  is  given  by 


o  -  n+(yp/un)p.  (24) 

Since  we  put  the  origin  (x=0 )  at  the  moving  growth 
surface  the  impurity  flow  is 

3NV 

Jk  -  -°kSir  +  Wk^k  +  vNk-  (25) 


The  continuity  equation  (7)  can  be  written 


32N.  a 

3k jjr  "  zkBkDkk[ENk]  -  v^r* 


(26) 


There  is  an  initial  condition; 


The  boundary  conditions  are 


Nk(-,t)  -  Nk(«,0)  -  Nj^M/  (28) 

and  (10)  with  vaa  and  neglected,  and  subscript  "k" 

do  9d 

for  impurity  type  reintroduced. 

The  electric  field  problem  and  the  redistribution 
problem  constitute  a  system  of  two  coupled  second  order 
differential  equations,  the  former  ordinary  and 

intrinsically  very  nonlinear,  the  latter  partial  and 
nonlinear  in  the  sense  that  the  field  depends  in  a 
complicated  way  on  all  the  impurity  concentrations . 

Without  making  simplifications  this  calls  for  numerical 
treatment.  The  results  will  be  interpreted  analytically 
in  section  7.3. 


7.2  HIGH  TEMPERATURE  RESULTS 

We  have  chosen,  as  an  initial  application  of  the 
model,  to  simulate  two  epitaxial  growth  processes  in  which 
the  effect  of  interfacial  fields  might  be  significant. 
They  are: 

1.  Growth  of  lightly  doped  n-type  GaAs  on  a  heavily  doped 
n-type  substrate;  and 

2.  Growth  of  moderately  doped  n-type  GaAs  on  a 
compensated  semi-insulating  substrate. 
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These  processes  are  interesting  because  they  are 
representative,  respectively,  of  the  technologies  for 
fabricating  Gunn  effect  devices  and  MESFET  material.  We 
have  chosen  realistic  values  of  initial  impurity 
concentrations  and  of  growth  velocity,  and  what  we  believe 
are  reasonable  estimates  of  ionic  mobilities  and  diffusion 
constants  [8].  The  resulting  impurity  profiles  show  some 
striking  inhomogeneities  in  the  distribution  of  minority 
species.  In  the  two  runs  we  have  simulated  in  detail  the 
redistribution  does  not  generate  harmful  conductivity 
regions.  It  is  possible,  however,  to  recognize 
combinations  of  parameters  that  may  prove  troublesome  in 
applications  as  we  will  give  an  example  of  later. 


Before  discussing  the  numerical  simulations  in  detail, 

we  can  point  out  some  features  that  they  have  in  common. 

We  assume  a  high  temperature  acceptor-like  surface  state 

density  as  calculated  by  Wolfe  and  Nichols  [2]  of 
11  —2 

N  =6.3x10  cm  .  The  surface  electric  field  associated 

4 

with  the  corresponding  charge  density  is  E  =-8.6x10  V/cm. 

In  Figures  3  to  16  all  concentrations  are  given  in 

units  of  the  measured  [9]  1000K  intrinsic  carrier 

16  “3 

concentration  of  n^=6xl0  cm  ,  distances  in  units  of  the 

-2 

intrinsic  Debye  length  1-^=3. 2x10  pm,  potential  in  units 
of  the  thermal  voltage  VT=86.2  mV,  electric  field  in  units 
of  ET=2.66xl04V/cm,  and  time  in  minutes.  The  resulting 


units  of  diffusion  constant  and  velocity  are, 

-13  2 

respectively,  D0=l. 75x10  cm  /s  and  v0=1.94um/h. 

We  find  that,  for  any  reasonable  growth  velocity,  the 
growth  surface  rather  quickly  outruns  any  effects 
associated  with  the  layer-substrate  interface.  This  means 
that  the  interfacial  field  rapidly  divides  into  two 
regions.  The  chief  effect  of  the  field  region  associated 
with  the  growth  surface,  after  the  initial  growth  stage, 
is  to  establish  a  surface  concentration  profile.  If,  as 
is  usual,  the  growth  rate  is  held  constant,  this  profile 
does  not  change.  At  the  interface  the  field  will  have 
different  effects  on  minority  and  majority  species.  The 
two  regions  will  be  discussed  in  detail  in  section  7.3 

The  separation  of  the  two  inter facial  field  regions 
has  another  interesting,  and  not  a  priori  obvious  bearing 
on  the  growth  surface  boundary  condition.  We  have  assumed 
a  fixed  surface  field,  and  we  find  that,  after  the  initial 
stage,  the  surface  potential  remains  effectively  constant 
(Figures  3  and  11).  Obviously  the  converse  would  be 
true  as  well.  Thus,  knowledge  of  the  actual  and  poorly 
understood  high  temperature  surface  state  energy 
distribution  which  determines  the  relation  between  the 
surface  field  and  potential,  is  not  essential  for  a 
realistic  simulation. 


For  some  of  the  impurity  profiles  ripples  occur  close 
to  the  layer-substrate  interface.  These  are  numerical 
artifacts  and  are  discussed  in  Chapter  3. 

7.2.1  Growth  on  n+  Substrate  (Table  1) 

When  the  n-type  substrate  is  more  heavily  doped  than 
the  epitaxial  layer,  the  interfacial  field  points  from  the 
substrate  to  the  epitaxial  layer  (Figures  3  and  4). 
The  substrate  is  partially  compensated,  its  donors  (Figure 
5)  experience  a  force  pulling  them  out  of  the  substrate, 
its  acceptors  (Figure  6)  a  force  pushing  them  back. 
Hence  the  electric  force  increases  the  concentration 
gradient  for  the  donors  while  decreasing  it  for  the 
acceptors.  If  the  diffusion  constant  for  the  two  species 
are  roughly  equal,  this  means  that  the  outdif fusion  of  the 
acceptors  lags  behind  that  of  the  donors,  as  can  be  seen 
in  Figure  10.  If  the  acceptors  are  an  order  of 
magnitude  more  mobile  than  the  donors,  we  find  (Figure 
9)  <  that  very  little  charge  separation  occurs. 

The  opposite  effect  is  experienced  by  the  donors 
(Figure  7)  and  acceptors  (Figure  8)  grown  into  the 
epitaxial  layer.  Here  interfacial  field  and  concentration 
gradient  drive  the  acceptors  in  the  same,  the  donors  in 
the  opposite  directions.  The  somewhat  surprising 
concentration  profiles  of  Figures  7  and  8  result  from 
this  interaction.  These  profiles  can  be  understood  when 
it  is  realized  that  the  field  is  mostly  confined  in  the 


Table  1  First  n-GaAs  growth  simulation:  Lightly  doped 
epitaxial  layer  on  heavily  doped  substrate. 
v=10ym/h. 


Impurity  Intended/Initial  D  B 

Ml _  M _ i. _ _  / - "  /.  \ 


No. 

Type 

Concentration 

(cm”3) 

( cm2/s ) 

1 

substrate  donor 

lxlO18 

lxlO"14 

1.27 

2 

substrate  acceptor 

5xl017 

lxlO"13 

1.27 

3 

incorporated  donor 

lxlO15 

lxlO"14 

1.27 

4 

incorporated  acceptor 

2.5xl014 

lxlO"13 

1.27 
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Figure  3 
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The  electrostatic  potential  for  the  growth  case 
listed  in  Table  1,  after  0,  1,  5  and  10 
minutes  of  growth. 
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Figure  4 
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The  electric  field  for  the  growth  case  listed 
in  Table  1,  after  0,  1,  5  and  10  minutes  of 
gr 
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NORMALIZED  DISTANCE  FROM  SURFACE 

Figure  5  The  distribution  of  impurity  1  for  the  growth 
case  listed  in  Table  1,  after  1,  5,  and  10 
minutes . 
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N0RMAUZED  DISTANCE  FROM  SURFACE 

Figure  6  The  distribution  of  impurity  2  for  the  growth 
case  listed  in  Table  1,  after  1,  5  and  10 
minutes . 
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NORMALIZED  DISTANCE  FROM  SURFACE 

Figure  7  The  distribution  of  inpurity  3  for  the  growth 
case  listed  in  Table  1,  after  1,  5,  and  10 
minutes . 


NORMALIZED  DISTANCE  FROM  SURFACE 

Figure  8  The  distribution  of  impurity  4  for  the  growth 
case  listed  in  Table  1,  after  1/  5  and  10 
minutes . 
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NORMALIZED  DISTANCE  FROM  SURFACE 

Figure  9  Superposition  of  the  impurity  profiles  after 
10  minutes  for  the  growth  case  listed  in  Table 
1.  (The  impurity  concentrations  are  plotted 
to  different  scale). 


NORMALIZED  DISTANCE  FROM  SURFACE 

Figure  10  Same  as  Figure  9  but  for  a  slightly, 

different  growth  case  with  D1=D2=D3=D4=10_14cm/s 


more  lightly  doped  epitaxial  region.  Hence,  while  it 
influences  only  the  edge  of  the  substrate  impurity 
profiles,  it  can  cause  a  pile-up  of  donors,  and  "eat  a 
hole"  into  the  distribution  of  acceptors  in  the  epitaxial 
layer.  It  is  the  fact  that  these  impurities  are  minority 
species  that  allows  the  non-diffusion-like  profiles  to 
occur,  as  we  will  see  in  section  7.3.- 

The  dominant  effect  in  this  simulation  is  the 
field-enhanced  outdif fusion  of  donors  from  the  substrate. 
Since  this  species  has  a  concentration  from  one  to  three 
orders  of  magnitude  greater  than  any  other,  it  blots  out 
the  variation  of  all  the  minority  species,  so  far  as  the 
room  temperature  carrier  concentration  is  concerned.  The 
effect  posited  by  Wolfe  and  Nichols  [1,2]  to  account  for 
an  interfacial  high  resistivity  layer  would  be  expected  to 
occur  if  the  substrate  donors  were  effectively  immobile. 
Then,  outdiffusion  of  partially  compensating  acceptors 
could  produce  a  thin  nearly  intrinsic  or  polarization 
reversed  layer.  This  will  be  discussed  in  more  detail 
be  low . 

7.2.2  Growth  on  a  Compensated  Semi -Insulating  Substrate 
(Table  2  ) 

We  have  attempted  to  model  a  Cr-compensated  substrate 

as  follows:  A  concentration  of  8xl016cm“3  slowly  diffusing 
-14  2 

(D=10  cm  /s)  residual  donors  (for  example  Si)  is 

17 

compensated  by  1x10  cm  Cr  of  which  90%  are  totally 


Table  2  Second  n-GaAs  growth  simulation:  Moderately 
doped  epitaxial  layer  on  semi-insulating 
substrate.  v=10um/h.  The  substrate,  is 
initially  perfectly  compensated  by  9xlO±'>cm_J  Cr 
acceptors  which  are  treated  as  fixed. 


Impurity  Intended/Initial  D  B 


No. 

Type 

Concentration 

(cm-3) 

(cm  /s) 

1 

substrate  donor 

8xl016 

lxlO"14 

1. 

2 

substrate  donor 

lxlO16 

lxlO"12 

1. 

3 

incorporated  donor 

lxlO17 

lxlO"14 

1. 

4 

incorporated  acceptor 

2.5xl016 

lxlO"13 

1. 

immobile  substitutional  acceptors,  10%  mobile 
-12  2 

(D=10  cm  /s)  interstitial  donors.  This  is  not  a 
complete  model  since,  for  instance,  conversion  between 
interstitial  and  substitutional  Cr  is  neglected.  It  will, 
however,  illustrate  the  effect  of,  particularly,  the 
surface  field  on  the  redistribution  of  interstitial  Cr. 
In  Section  5  a  more  elaborate  and  and  comprehensive  model 
of  Cr  in  GaAs  is  -  presented.  The  epitaxial  layer 
grown  on  this  substrate  incorporates  lxl017cm~3  slowly 
diffusing  (D-10  14cm2/s)  donors  and  2.5xl016cm  3  an  order 
of  magnitude  more  mobile  acceptors,  giving  a  net  n-type 
doping  of  7.5xl016cm“3.  With  this  concentration,  barely 
in  excess  of  the  intrinsic  carrier  concentration  of 
6xl016cm”3  at  1000K,  the  epitaxial  layer-substrate 
interfacial  potential  (Figure  11)  and  field  (Figure 
12)  are  quite  small  and  modify  the  diffusion  only 
slightly.  The  chief  effect  of  the  impurity  redistribution 
comes  from  the  outdif fusion  of  Cr-donors  (Figure  14)  and 
the  preferential  indiffusion  of  acceptors  into  the 
substrate  (Figure  16),  which  turn  the  substrate  slightly 
p-type  near  the  surface.  As  a  result,  upon  cooling  to 
room  temperature,  the  Debye  shielding  length  is  reduced 
below  what  it  would  be  without  redistribution,  and  the 
carrier  concentration  falls  off  considerably  more  steeply 
as  will  be  illustrated  below,  - 
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NORMALIZED  DISTANCE  FROM  SURFACE 

Figure  11  The  electrostatic  potential  for  the  growth 
case  listed  in  Table  2,  after  0,  1,  5  and  10 
minutes  of  growth. 
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NORMALIZED  DISTANCE  FROM  SURFACE 

Figure  12  The  electric  field  for  the  growth  case  listed 
in  Table  2,  after  0,  1,  5  and  10  minutes  of 

growth . 
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Figure  13  The  distribution  of  impurity  1  for  the  growth 
case  listed  in  Table  2,  after  1,  5,  and  10 

minutes . 


Figure  14  The  distribution  of  impurity  2  for  the  growth 
case  listed  in  Table  2,  after  1,  5  and  10 
minutes . 
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NORMALIZED  DISTANCE  FROM  SURFACE 
Figure  15  The  distribution  of  impurity  3  for  the  growth 
case  listed  in  Table  2,  after  1,  5,  and  10 

minutes . 


Figure  16  The  distribution  of  impurity  4  for  the  growth 
case  listed  in  Table  2,  after  1,  5  and  10 
minutes . 


The  epitaxial  material  for  use  in  MESFET's  and  IC's  is 
frequently  one  to  two  orders  of  magnitude  more  lightly 
doped  than  our  numerical  example.  Such  epitaxial  layers 
are  intrinsic  at  1000K,  and  their  growth  on 

semi-insulating  substrates  would  occur  without  an 

inter facial  field.  Under  these  conditions  preferential 
diffusion  [10]  of  donors  from  a  heavily  compensated 
substrate,  as  shown  in  Figure  13  could  produce  a  high 
conductivity  layer  at  the  interface. 

7.2.3  Room  Temperature  Potential  and  Carrier  Concentration 

As  mentioned,  neither  of  the  two  growth  simulations 
result  in  any  harmful  conductivity  regions.  In  the  second 
case  the  redistribution  actually  makes  the  electron 
concentration  in  the  substrate  fall  off  even  faster  than 
without  redistribution,  as  we  will  see  in  this  section. 

We  calculate  the  potential  distribution  at  room 
temperature  (300K)  after  the  high  temperature  (1000K) 
redistribution  that  we  studied  for  10  minutes  in  the  last 
two  sections.  Our  interest  is  to  see  how  the 
redistribution  itself  affects  the  potential  <Hx),  and  thus 
the  electron  and  hole  concentrations  (23).  We  therefore 
set  the  surface  field  to  zero  (no  charged  surface  states). 
We  use  normalized  units  as  before  but  now  using  the  room 
temperature  values;  position  x  is  expressed  in  terms  of 
(300K)*0 .34cm,  potential  in  terms  of  VT(300K)=26mV, 


and  concentrations  in  terms  of  the  intrinsic  carrier 
concentration  ni(300K)=1.6xl06cm“3. 

In  the  first  growth  case  we  do  not  get  any  harmful 
conductivity  regions  since  the  substrate  donors  and 
acceptors  redistribute  in  a  similar  fashion.  However,  we 
also  calculate  the  potential  that  would  result  if  the 
substrate  donors  were  fixed  but  the  rest  of  the  impurities 
had  the  diffusion  properties  listed  in  Table  '1.  We  do 
not  recalculate  the  redistribution.  Instead  we  use  the 
final  distribution  of  the  other  impurities  computed  during 
the  simulation,  but  assume  a  fixed  substrate  donor 
distribution.  This  is  of  course  not  quite  correct  since  a 
fixed  majority  dopant  concentration  would  result  in  a 
practically  constant  field  at  the  interface  as  opposed  to 
the  first  simulation,  where  the  peak  field  decreases 
(Figure  4).  This  would,  for  instance,  cause  less 
outdif fusion  of  the  substrate  acceptors.  Nevertheless,  it 
will  give  an  idea  of  the  result  of  a  rigorous  calculation. 
The  result  is  shown  in  Figure  17.  The  material  does  not 
become  p-type,  but  there  is  a  considerable  lowering  of  the 
potential  due  to  the  outdif fusion  of  acceptors  from  the 
substrate,  and  thus  an  even  larger  (about  four  orders  of 
magnitude)  lowering  of  the  electron  concentration  and 
conductivity  in  quite  a  wide  region.  A  rigorous 
calculation  would  make  the  effect  less  pronounced. 


Figure  17  Room  temper ature  potential  t  as  a  function  of 
the  distance  x  from  the  growth  surface 
resulting  from  the  first  growth  simulation 
(Table  1).  The  material  would  be  intrinsic 
for  4»-26. 

a.  with  no  redistribution. 

b.  after  redistribution. 
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When  we  do  the  calculation  for  the  potential  at  room 
temperature  resulting  from  the  second  growth  simulation 
we  get  the  peculiar  result  shown  as  curve  b  in  Figure 
18.  Curve  a  is  the  potential  for  no  redistribution  and 
this  falls  off  very  slowly  (note  that  the  abscissae  is 
lg(x  ).  In  the  relatively  moderately  doped  epitaxial 
layer  (lg(x)<-3.31),  where  the  fast  outdiffusion  of 
donors  does  not  affect  the  total  doping,  curve  a  and  b 
coincide  as  expected.  For  -3.3<lg(x)<-3.1  curve  b  falls 
off  much  faster  than  curve  a  due  to  the  uncovered 
acceptors  in  the  substrate.  This  means  that  the 
conductivity  falls  off  faster  than  "expected"  which  is 
advantageous  in  this  case.  Beyond  lg(x)*-3.1  the  result 
is  incorrect  due  to  errors  in  the  concentrations  that  we 
have  allowed  in  the  simulation.  For 
lg(x)<-3.04  the  total  net  doping  should  be  practically 
zero  since  at  lg(x)=-3.04  the  total  net  doping  ceases  to 
be  negative.  However,  the  numerical  treatment  produces  a 
small  positive  doping  in  the  region  -3.04<lg(x)<-2.70  and 
a  small  negative  doping  in  the  region  -2.70<lg(x)<-1.41. 
At  room  temperature  these  errors  get  effectively  much 
larger  since  we  measure  the  concentrations  relative  to  the 
intrinsic  one.  The  reason  that  similar  effects  do  not 
occur  in  the  first  simulation  is  that  there  the  material 
is  more  heavily  doped,  and  the  errors  in  potential  are 
essentially  inversely  proportional  to  the  doping 


Figure  18  Room  temperature  potential  V  as  a  function  of 
the  distance  x  from  the  growth  surface 
resulting  from  the  second  growth  simulation 
(Table  2).  The  material  is  intrinsic  for 
¥»*0. 

a.  with  no  redistribution. 

b.  after  redistribution.  Net  doping  set  to 
zero  for  lgx>-2.70. 

c.  after  redistribution.  Net  doping  set  to 
zero  for  lgx>-3.04. 


Curve  c  shows  the  result  when  we  require  the 
doping  to  be  zero  beyond  the  point  lg(x)=-3 . 04,  where  the 
simulation  says  that  the  doping  becomes  positive  (which  it 
should  not).  This  curve  looks  more  realistic;  it 
coincides  with  curve  b  in  the  substrate  close  to  the 
interface,  where  the  p-conversion  is  most  pronounced, 
while  coinciding  with  curve  a  deep  in  the  substrate  fear 

I 

away  from  the  region  of  redistribution.  The  region  where 
curve  b  and  c  coincide  shows  that  the  potential  actually 
falls  off  much  faster  than  it  would  in  the  case  of  no 
redistribution. 

7.3  QUALITATIVE  ANALYTICAL  INTERPRETATION 

The  numerical  results  clearly  illustrate  that  the 
problem  separates  into  two  practically  independent  pieces, 
one  associated  with  the  surface  and  the  other  with  the 
layer-substrate  interface.  In  this  section  we  attempt  an 
analytical  interpretation  of  the  results  based  on  this 
observation. 

7:^3. 1  Surface  Impurity  Profiles 

In  the  simulations  presented,  using  a  constant  growth 
velocity,  the  surface  profiles  very  quickly  become 
stationary.  The  choice  of  growth  velocity,  ionic  mobility 
and  surface  electric  field  allows  us  to  solve  the 
stationary  problem  approximately.  Assuming  steady  state 
and  applying  the  boundary  condition  (10),  again 
neglecting  vas  and  v3a,  the  continuity  equation  (26) 


becomes  (suppressing  subscripts  for  the  impurity  type) 


dN 

3x  ' 


[zBE+g]N 


-§v 


(29) 


The  solution  of  this  is  straightforward; 

m 

N(x)  -  Nae”zB,<'(x>y*e“5ezB,,'(x'Hv5)dC.  (30) 

0 

Here  we  have  used  the  fact  that  the  surface  profile  is 
practically  independent  of  the  interfacial  impurity 
profile  by  extending  the  upper  limit  of  integration  to 
infinity.  zB^(x+(D/v)C)  goes  to  zero  exponentially  as  £ 
becomes  large.  The  largest  contribution  to  the  integral 
comes  from  small  £.  We  make  the  approximation 


’Mx+^5)  *  ♦  <*>  *  E(x)^£,  (31) 

which  is  expected  to  be  better  the  smaller  zBDE(x)/v  is 
compared  to  1.  With  this  approximation  we  get 

N(x)  *  NaU+^|l2l]“1.  (32) 

This  is  the  same  result  as  the  zeroth  order  approximation 
in  a  perturbation  scheme  of  (29),  in  which  dN/dx  is 
considered  the  perturbation.  In  Figures  19  and  20  the 
surface  profiles  of  the  incorporated  donors  and  acceptors 


(impurity  3  and  4  in  Table  1,  respectively)  are  plotted 
together  with  the  prediction  of  equation  (32)  using  the 
numerical  results  for  the  electric  field.  This  is 
therefore  also  a  check  of  the  consistency.  The  agreement 
is  good  in  both  cases  exept  for  the  donors  in  an  extremely 
thin  surface  layer.  Mathematically  this  is  real  (a  so 
called  boundary  layer  [11]),  but  it  is  of  no  physical 
importance  in  this  case  since  it  is  only  about  a  lattice 
parameter  thick. 

The  treatment  here  is  basically  the  same  as  in  [1], 
but  we  use  the  boundary  condition  (10),  and  make  the 
approximation  (31)  which  appears  to  be  valid  for  large 
growth  velocity. 

7.3.2  interfacial  Region 

After  the  interface  has  separated  from  the  surface  it 
can  be  treated  as  a  moving  interface  between  two 
semi'* infinite  regions,  in  which  case  it  makes  more  sense 
to  fix  the  frame  of  reference  at  the  interface.  The 
continuity  equation  then  reads 

H  -  -  zBd|j[EN]  .  (33) 

dX 

If  N  describes  the  concentration  of  the  majority 
species,  and  charge  neutrality  is  assumed  so  that 


N(x)/N 


Figure  20  Surface  profile  of  impurity  4  for  the  growth 
case  listed  in  Table  1. 
a  numerical  result  for  the  concentration, 
o  (32)  using  the  numerical  result  for  the 
electric  field. 


one  can  [12]  rewrite  (33)  as  a  diffusion  equation; 


If  <35> 

with  the  concentration  dependent  diffusion  constant 
D‘ (N)  -  D[1+z2BN(4+N2)~Z] .  (36) 


We  see  that  the  diffusion  constant  is  smaller  for  smaller 
concentration.  This  is  observed  in  Figure  5  as  the 
concentration  gradient  being  steeper  in  regions  of  lower 
concentration. 

If  N  describes  a  minority  species,  the  field  is 
practically  independent  of  N  and  is  determined  by  other 
species.  We  could  treat  the  drift  term  in  (33)  as  a 
perturbation.  Assuming  constant  initial  doping  Ns  in  the 
substrate  (x>0)  and  constant  intentional  doping  Ne  of  the 
epitaxial  layer  (x<0),  the  zeroth  order  solution  is  [13] 

N„(x,t)  -  £tVV  +  ItV',alerf(57TDt))-  (37) 


The  first  order  correction  Nj  is  obtained  by  solving  the 
equation 


zbd|^[en0]  . 


(  38) 


Just  by  looking  at  the  structure  of  this  equation,  some 
information  can  be  extracted.  (38)  is  a  diffusion 
equation  for  with  a  known  generation  term 
-zBD( 3(ENq)/3x) .  For  the  donors  in  Figure  7  it  will  act 
as  a  source,  while  for  the  acceptors  in  Figure  8  it  will 
act  as  a  sink  on  the  epitaxial  layer  side  of  the 
interface,  leading  to  the  observed  inter facial  structure. 

The  separation  of  the  problem  observed  in  our 
numerical  treatment  suggests  that  one  would  get  a  good 
picture  of  the  redistribution  if  one  treated  the  surface 
and  interface  problem  separately  when  the  situation  is  one 
of  sufficiently  fast  growth  and  constant  intentional 
doping  of  the  epitaxial  layer.  This  is  certainly  a  common 
situation  and  would  allow  a  substantial  reduction  in 
excecution  time  due  to  the  absence  of  moving  boundaries  in 
such  an  approach. 

7.4  QUASI -CHEMICAL  REACTIONS 

It  is  well  known  that  many  impurities  participate  in 
reactions  which  greatly  influence  the  redistribution. 
Since  GaAs  is  of  particular  interest  to  us,  a  good  example 
is  Zn,  which  diffuses  interstitially,  but  can  get  captured 
by  Ga  vacancies.  If  the  reactions  are  in  local 
equilibrium  the  analysis  is  simplified  since  the  mass 
action  laws  give  relations  between  impurity  concentrations 
that  u.  e  valid  locally.  Zn  diffusion  has  been  analysed 
rather  successfully  in  this  way  [8],  leading  to  a 


concentration-dependent  diffusion  constant.  However,  if 
the  reactions  are  not  in  equilibrium,  we  have  to  solve  the 
continuity  equations  with  proper  generation  and 
recombination  terms  included.  This  will  enhance  the 
coupling  between  the  redistribution  of  the  different 
species  that  already  exists  in  our  model  through  the 
electric  field.  Nonlinearities  will  be  introduced  which 
will  necessitate  a  less-than- implicit  numerical  scheme, 

Extension  of  our  model  would,  in 
principle,  be  straightforward,  but  it  is  conceivable  that 
in  practice  one  might  run  into  problems  with  stability  and 
convergence  in  the  numerical  implementation.  The  primary 
problem,  however,  is  to  model  what  reactions  go  on,  and  at 
what  rate.  In  an  attempt  to  do  so,  we  started  out  to 
study  the  diffusion  of  Cr  in  GaAs .  This  system  is  of 
technological  importance  since  it  earn  produce 
semi-insulating  material  useful  as  substrates  for  IC's. 
It  is  particularly  interesting  since  Cr  exhibits 
remarkable  redistribution  behavior  in  GaAs,  unlikely  to  be 
accounted  for  by  a  simple  diffusion  mechanism.  The 
problem  proved  to  be  a  complicated  one,  and  is  at  this 
point  not  sufficiently  well  understood  to  be  studied  by 
our  general  numerical  redistribution  model. 


7.5  CONCLUSIONS 


In  this  work  some  aspects  of  impurity  redistribution 
in  semiconductors  during  processing  at  elevated 
temperatures  has  been  discussed. 

We  have  presented  a  general  nonequilibrium  model  for 
the  field-assisted  diffusion  during  epitaxy  and  annealing. 
It  was  demonstrated  that  even  under  very  simple  epitaxial 
growth  conditions,  ionic  redistribution  can  produce  fairly 
complicated  impurity  concentration  profiles,  particularly 
for  minority  species.  We  also  discussed  how  harmful 
conductivity  regions  could  come  about  as  a  result  of 
outdif fusion  from  the  substrate.  The  model  should  prove 
useful  in  simulating  a  wide  variety  of  situations  since 
the  numerical  implementation  contains  options  such  as 
inhomogeneous  substrate  doping,  depletion  of  the  ambient 
and  changes  of  growth  rate.  In  addition,  it  appears 
straightforward  to  include  reactions,  which  would  be 
particularly  interesting  when  trying  to  simulate 
redistribution  during  post- implantation  annealing.  Under 
some  simple  epitaxial  conditions,  improved  execution 
speed  could  be  accomplished  by  treating  the  surface  and 
interface  regions  separately. 
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8.  INCORPORATION  OF  AMPHOTERIC  IMPURITIES 


The  Group  IV  impurities,  namely  C,  Si,  Ge,  Sn,  and  Pb,  are  expected  to 
behave  either  as  donors  or  acceptors  in  GaAs,  depending  on  whether  they 
substitute  on  Ga  or  As  sites.  It  is  well  known  that  the  relative 
incorporation  as  donors  and  acceptors  is  different  for  different  elements, 
growth  techniques,  and  growth  conditions  (e.g.  growth  temperature,  III/V 
ratio,  etc.)  Much  of  this  information  has  come  from  doping  experiments 
which  rely  on  Hall  effect  data  and  the  tacit  assumption  that  only  changes 
in  dopant  incorporation  are  important  in  determining  changes  in  carrier 
concentration  and  77  K  mobility.  The  direct  measurement  of  the  concentra¬ 
tions  of  individual  impurity  species  would  afford  a  clearer  picture  of 
impurity  incorporation,  particularly  at  low  concentrations.  To  this  end, 
photothermal  ionization  spectroscopy  (described  in  detail  by  Stillman  et 
al.  (1977)),  and  low  temperature  photoluminescence  have  been  used  to 
measure  the  concentrations  of  donors  and  acceptors,  respectively,  in  high 
purity  GaAs  prepared  by  a  wide  variety  of  growth  techniques.  The  relative 
donor  concentrations  were  determined  from  the  1  s-2p(m—l)  peak  amplitudes 
at  high  magnetic  fields  using  the  identifications  of  Uolfe  et  al.  (1976) 
for  Pb,  Sn,  Se,  and  of  Low  et  al.  (1982  a  and  b)  and  Ozeki  et  al.  (1977) 
for  Xi»Si,  X2“S,  and  Xj-Ge.  The  relative  acceptor  concentrations  were 
estimated  from  the  relative  amplitudes  of  the  (D°-A°)  peaks  at  2  K  under 
1  mU  excitation  from  an  Ar+  laser,  using  the  identifications  in  Ashen  et 
al.  (1975).  Absolute  donor  and  acceptor  concentrations  could  then  be 
determined  using  the  total  Nq  and  NA  calculated,  under  the  assumption  that 
the  shallow  impurities  dominate  the  ionized  impurity  scattering,  from  the 
77  K  Hall  effect  data  by  the  method  of  Wolfe  et  al.  (1970). 


8.1  LIQUID  PHASE  EPITAXY 

The  residual  shallow  impurities  in  high  purity  undoped  LPE  samples  from 
seven  different  sources  were  measured  (see  Skromme  et  al.  1983a).  The 
photothermal  ionization  and  photoluminescence  spectra  for  the  Hewlett 
Packard  sample  are  shown  in  Fig.  1  (a)  and  (b),  and  are  typical  of  the 
spectra  of  the  other  samples  measured.  The  characteristic  LPE  residual 
donors,  Pb,  Xj-Si,  Sn,  and  X2“S  are  present  at  concentrations  of  0.28, 
0.22,  0.12  and  1.54  xlO^cm-3  respectively.  The  dominance  of  X2“S  is 
usual  for  both  graphite  (Skromme  et  al.  1983a)  and  silica  (Cooke  et  al. 
1978)  boat  grown  LPE  GaAs.  The  characteristic  LPE  residual  shallow 
acceptors,  C,  Mg,  Si  and  Ge  are  present  at  approximate  concentrations  of 
0.13,  0.18,  1.0  and  0.03  xl014cm-3  respectively. 

The  photothermal  ionization  and  photolumineecence  data  for  the  LPE  samples 
measured  consistently  indicates  that  the  Group  IV  impurities  Pb  and  Sn 
incorporate  preferentially  as  donors  while  Si,  Ge,  and  C  incorporate  pre¬ 
ferentially  as  acceptors.  The  Incorporation  of  Pb  in  GaAs  grown  by  AsCl-j- 
VPE  is  small  (Wolfe  et  al.  1976a)  and  although  Pb  doping  experiments  have 
associated  Pb  with  a  peak  coincident  with  that  of  Mi  deep  acceptors,  it 
seems  likely  that  this  can  be  explained  by  traces  of  Mi  present  in  the  Pb 
dopant  (Bebb  et  al  1972),  and  that  Pb  may  not  Incorporate  as  an  acceptor 
in  GaAs.  It  should  be  noted  however  that  the  Mn  acceptor  peak  was 
observed  in  most  of  the  LPE  samples  measured.  The  Pb  and  Sn  donor  peaks 
were  sometimes  larger  in  the  spectra  of  the  other  LPE  samples  than  that 
shown  in  Fig.  1  but  were  always  dominated  by  sulfur.  In  the  two  Cornell 
LPE  samples  in  which  Sn  was  most  prominent  (nearly  equal  in  concentration 


Fig.  1  (a)  Photothermal  ionization  and  (b)  photoluminescence 

spectra  of  Hewlett  Packard  LPE  sample. 


to  S)  traces  of  Sn  acceptors  were  detected  (via  the  neutral  Sn  acceptor 
bound  exciton  line),  but  no  Sn  acceptors  could  be  detected  in  any  of  the 
other  LPE  samples.  It  was  impossible  to  determine  the  Sn  acceptor 
concentration  so  from  this  data  the  incorporation  ratio  [Snga] / [SnAs]  can 
only  be  said  to  be  much  greater  than  one.  The  Xj«Si  donor  peak  was 
sometimes  absent  from  the  photothermal  ionization  spectra  and  this 
correlated  with  a  much  smaller  level  of  the  usually  dominant  Si  acceptor 
detected  by  photoluminescence.  In  the  three  samples  for  which  both  Si 
donor  and  acceptor  concentrations  could  be  accurately  measured,  (each  of 
which  was  grown  at  700°C)  incorporation  ratios  of  [SiGa] / [SiAs]  -  0.16, 
0.19,  and  0.22  were  determined.  The  X3«Ge  donor  peak  has  never  been 

observed  in  LPE  GaAs  but  small  concentrations  (0.01-0.10  xlO^cm-^)  of  Ge 
acceptors  were  observed  in  most  of  the  LPE  samples  measured.  The 

detection  limit  for  X3»Ge  donors  in  the  sample  with  the  largest  Ge 

acceptor  concentration  can  be  estimated  as  about  0.02x10*^,  and  so 
[GeGa]/[Ge^8]  must  be  less  than  about  0.2.  Early  doping  experiments 
suggested  that  C  donors  may  contribute  a  photothermal  ionization  peak 

coincident  with  the  X3*Ge  peak  (Wolfe  et  al.  1976).  There  is  evidence  to 
suggest  that  C  does  not  incorporate  as  a  donor  at  least  in  ASCI3-VPE  GaAs 
(Ozeki  et  al.  1977),  but  even  assuming  that  C  donors  do  contribute  such  an 
X3  peak,  the  detection  limit  for  X3  donors  together  with  the  largest 
concentration  of  C  acceptors  detected  determine  an  upper  bound  on 
l^Gal/l^Asl  *n  GaAs  of  about  0.02. 

These  observations  of  the  relative  incorporation  of  the  Group  IV 
impurities  [^Qa]/[IV^S]  are  in  qualitative  agreement  with  the  well 
established  behavior  of  Si,  Ge  and  Sn  as  dopants  in  LPE  deduced  from  Hall 
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effect  data,  and  have  been  predicted  quantitatively  for  these  Impurities 

O 

in  terms  of  Ga  and  As  activities  and  the  free  energy  difference  AG^yg 
between  IVGa  and  IVAg  site  occupation  (Teramoto  1972)  as: 

log  ([IVGa]/[IVAs])  -  log  (X*s/x£a)  +  e/T 

°  9  9 

where  e  is  determined  by  AGA/g  and  the  heats  of  solution  AHqs  and  AHAg. 

The  melt  mole  fractions  XAg  and  Xga  can  be  determined  at  a  given  growth 

O 

temperature  from  the  liquidus  data  of  Hall  (1963).  Using  the  AG^/g  values 
for  Si,  Ge,  and  Sn  calculated  by  Teramoto  and  the  typical  growth  tempera¬ 
ture  of  700°C  we  obtain  [SiGa]/[SiAs]  -  0.16,  [GeGa]/[GeAg]  -  0.034,  and 
[SnGa]/[SnAg]  ■  4.7,  which  are  in  good  agreement  with  our  observed  results 
above.  |r 

Following  Teramoto  and  using  the  tetrahedral  covalent  radii  of  Pauling 

O 

(1960)  one  can  calculate  AGA/B  for  both  C  and  Pb  as  +14.75  kcal/mole  and 
-28.33  kcal/mole,  respectively.  For  a  growth  temperature  of  700°C,  we 
obtain  incorporation  ratios  of  [CGa] / f CAsl  -  6.4x10“^  and  [PbGa] / [PbAg]  - 
3.0x10  .  These  numbers  should  be  taken  as  very  approximate  in  view  of 

the  inadequacy  of  equilibrium  thermodynamics  to  describe  the  non-equi- 

O 

librium  process  of  crystal  growth,  and  of  the  way  in  which  AG^g  is 
calculated  and  its  sensitivity  to  choice  of  covalent  radii,  but  they  are 
in  agreement  with  the  observation  of  Pb  as  exclusively  donors  and  of  C 
exclusively  as  acceptors  in  GaAs.  The  physical  interpretation  of  these 
results  is,  since  the  "size"  of  the  Ga  site  is  larger  than  that  of  As,  the 

O 

lattice  strain  energy  contribution  to  AGA/B  then  forces  the  smaller  C  onto 
As  (acceptor)  sites,  and  forces  the  larger  Pb  onto  Ga  (donor)  sites. 
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8.2  AsC13  VAPOR  PHASE  EPITAXY 

The  residual  shallow  impurities  in  several  high  purity  undoped  ASCI3-VPE 
samples  from  seven  different  sources  were  measured.  The  bottom  photo- 
thermal  ionization  spectrum  in  Fig.  2  for  the  Motorola  sample  shows  the 
characteristic  ASCI3-VPE  residual  donors  Xj*Si,  X2“S,  and  X3*Ge,  which  are 
present  at  concentrations  of  0.47,  0.53,  and  0.27  xlO^  cm-'*  respec¬ 
tively.  These  same  residual  donors  were  observed  in  the  other  samples 
although  Si  was  sometimes  more  and  Ge  less  prominent.  Because  of  its  very 
high  purity  the  peaks  in  this  sample  are  extremely  well  resolved,  and  so 
at  a  given  magnetic  field  it  was  used  as  a  reference  for  peak  identifica¬ 
tions  in  many  of  the  other  samples  described  in  this  work.  The  photo¬ 
luminescence  spectrum  for  this  sample  in  Fig.  3  shows  the  typically 
prominent  residual  acceptor  Zn  (0.51  xlO*^  cm-^)  and  smaller  concen¬ 
trations  of  C  and  Ge  acceptors  (0.17  and  0.10  xlO*^  cm--*).  Residual  Si 
acceptors  were  not  detected  in  this  sample  but  were  detected  in  trace 
amounts  in  some  of  the  other  samples  measured.  Residual  Ge  acceptors  were 
frequently  observed  at  approximately  this  same  concentration  and  always 
dominated  the  traces  of  Si  acceptors  when  the  latter  were  present.  The 
observed  preferential  incorporation  of  Si  as  donors  rather  than  acceptors 
is  in  agreement  with  the  general  absence  of  Si  acceptors  in  the  photo¬ 
luminescence  data  of  Ashen  et  al.  (1975)  for  undoped  ASCI3-VPE  GaAs,  and 
is  in  qualitative  agreement  with  their  prediction,  under  the  typical 
A8CI3-VPE  growth  conditions  of  Tq  ■  750°C  and  P(AsCl3)  ■  6  x  10-^  atm, 
that  [SiGa]/[SiAs]  -  1.6  x  10^.  The  arguments  of  Ashen  et  al.  can  be 

O 

extended  to  the  other  Group  IV  impurities,  using  the  AGA/B  values  cal¬ 


culated  in' the  manner  of  Teramoto,  giving  incorporation  ratios  which  scale 


Ill/V  =  0.020 


lll/V  =0.029 


lll/V  =  0.040 


34 

Wavenumber  (cm~' ) 


Photothermal  ionization  spectra  for  the  Toshiba  MOCVD  samples 
(top  three  plots)  in  which  Sn  and  Ge  donor  incorporation  was 
suppressed  by  increasing  the  IIX/V  ratio.  Lower  plot  is  for 
the  Motorola  AsCIj-VPE  reference  sample. 
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Fig.  3  Photoluminescence  spectrum  for  Motorola 
AsCl^-VPE  sample. 
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for  a  given  growth  temperature  with  those  calculated  for  LPE.  •  We  obtain 

then  [CGa]/[CAg]-1.2xlO-2,  [GeGa]/[GeAg]-3.7xlO+2,  [SnGa]/[SnAg]-3.9xlO+4, 

and  [PbGa ] / [PbAg]  ■  1.9xl0+^  for  the  above  AsClj-VPE  growth  conditions. 

of 

The  measured  concentrations^  Si  and  Ge  donors  and  acceptors  indicate  that 
the  calculated  incorporation  ratios  for  these  elements  are  substantially 
too  large.  These  ratios  correctly  predict  however,  the  conduction  type  (n 
or  p)  when  Ge,  Si,  Sn,  Pb  (Wolfe  et  al.  1976a)  and  C  (Ozeki  et  al,  1977) 
are  used  as  dopants  in  ASCI3-VPE,  and  the  dominance  observed  here  of  Si 
and  Ge  donors  over  acceptors.  In  the  samples  measured  which  have  com¬ 
parable  Si  and  Ge  donor  concentrations,  they  also  correctly  predict  the 
observed  dominance  of  Ge  acceptors  over  Si  acceptors. 

8.3  ASH3  VABOR  PHASE  EPITAXY 

The  residual  shallow  impurities  in  high  purity  undoped  ASH3-VPE  samples 
from  three  different  sources  were  determined  and  these  measurements  are 
discussed  in  more  detail  by  Skromme  et  al.  (1983b).  The  three  residual 
donors  Xj*Si,  X2**S,  and  X3»Ge  were  observed  in  the  samples  from  each 
source  with  X2*S  always  dominant,  while  C  and  Zn  were  the  dominant  shallow 
acceptors.  Both  Si  and  Ge  acceptors  were  observed  at  low  concentrations 
in  the  University  of  Illinois  samples  but  were  absent  in  all  the  samples 
from  the  other  two  sources.  The  most  striking  behavior  with  respect  to 
Group  IV  impurities  was  observed  in  the  Han  scorn  AFB  samples  grown  under 
otherwise  similar-  conditions  but  varying  the  III/V  ratio  (i.e. 
P(HCl)/P(AsH3)) .  The  photothermal  ionization  spectra  of  these  samples  in 
Fig.  4  show  the  steady  decrease  in  Group  IV  donor  incorporation  (X^-Si  and 
X3*Ge)  relative  to  Group  VI  donor  incorporation  (X2“S)  as  the  III/V  ratio 
is  increased,  until  the  Group  IV  donors  are  barely  detectable  at  III/V  - 
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Fig.  4  Photothermal  ionization  spectra  for  AsH^-VPE 

samples  showing  suppression  of  Group  IV  donors 
Si  and  Ge  with  increasing  III/V  ratio. 


1.7.  A  similar  but  less  dramatic  trend  was  observed  in  photoluminescence 
for  the  CAs  to  Zn q3  acceptor  incorporation  ratio.  The  ratio  [C^g]/[ZnQa] 
varied  from  3.8  to  11  as  the  III/V  ratio  was  increased  from  0.50  to  1.7. 
These  data  are  consistent  with  the  idea  that  the  gas  phase  stoichiometry 
influences  the  relative  impurity  incorporation  through  the  relative  con¬ 
centrations  of  £a  and  As  vacancies  (incorporation  sites)  on  the  growing 
surface.  Even  in  the  sample  grown  with  the  highest  III/V  ratio,  however, 
neither  Si  nor  Ge  acceptors  could  be  detected. 

8.4  METAL-ORGANIC  CHEMICAL  VAPOR  DEPOSITION 

High  purity  undoped  MOCVD  samples  from  four  different  sources  were 
measured  and  some  of  this  data  is  decribed  in  more  detail  by  Low  et  al 
(1983c).  In  all  the  samples  measured  Xj-Ge  was  the  dominant  residual 
donor  while  smaller  concentrations  of  Xj»Si,  Sn,  X2*S,  and  rarely  Pb 
donors  were  also  observed.  The  dominant  residual  shallow  acceptor  was 
usually  C,  but  Zn  was  usually  prominent  and  occasionally  exceeded  C  in 
concentration.  Traces  of  Ge  acceptors  were  occasionally  detected.  The 
photothermal  ionization  spectra  of  three  samples  from  Toshiba,  all  grown 
under  similar  conditions  (Tq  *  660°C)  but  with  varying  III/V  ratio  (i.e. 
P(TMGa)/P(AsH3))  are  shown  in  the  upper  three  curves  Fig.  2.  The  only  two 
residual  donors  detected  in  these  samples  were  Sn  and  Ge,  both  from  Group 
IV.  The  Nq  values  determined  from  Hall  data  steadily  decreased  from  17.3 
to  11.2  to  6.6  xlO^  cm”^  (while  NA  remained  essentially  constant  at  4.3, 

4.4  and  4.1  xlO^  cm-^)  as  the  III/V  ratio  increased  from  0.020  to  0.029 
to  0.040  respectively.  This  behavior  is  similar  to  that  observed  in  the 
AsHj-VPE  samples  in  that  Group  IV  donor  incorporation  is  supressed  by 
increasing  the  gas  phase  III/V  ratio,  and  might  be  interpreted  in  terms  of 


vacancy  incorporation  site  availability  in  the  same  way.  It  is 
interesting  to  note  that  the  Sn  donor  concentration  seems  to  decrease  more 
rapidly  than  that  of  Ge  with  increasing  III/V  ratio.  No  trend  was 
observed  in  the  ratio  [CAg]/[ZnQa]  (which  assumed  the  values  36,  26,  and 
38,  respectively)  as  the  III/V  ratio  was  increased.  Finally,  the 
concentration  of  Ge  acceptors  decreased  from  6.5  xloH  cm“^  to 
3.6  xlO**  cm-^  to  undetectable  as  the  III/V  ratio  increased  (see  Fig.  5) 
while,  consistent  with  the  expected  relative  vacancy  concentrations  the 
lGeGal/ l®eAs^  ratio  increased  from  2.3  x  lCp  to  3.1  x  10^  to  £  4  x  lCp 
respectively.  The  increase  in  concentration  of  both  Ge  donors  and 
acceptors  with  decreasing  III/V  ratio  (increasing  P^sHj))  is  consistent 
with  the  suggestion  that  the  source  of  the  residual  impurity  Ge  in  MDCVD 
GaAs  may  be  the  ASH3  (possibly  in  the  form  of  GeH^). 

8.5  MOLECULAR  BEAM  EPITAXY 

High  purity  elemental  As  source  grown  MBE  samples  from  four  sources  were 
measured  and  some  of  this  data  is  discussed  in  more  detail  by  Low  et  al. 
(1982  a  and  d) .  The  dominant  and  usually  the  only  residual  donor  detected  in 
these  samples  was  X2*S  but  the  necessity  of  using  Si  or  Sn  dopants  to 
achieve  the  n-type  conduction  required  for  the  photothermal  ionization 

measurements  would  mask  their  presence  as  residuals.  No  donor  peak 

attributable  to  C  was  observed  in  any  of  these  samples  but  C  was  the  only 
shallow  acceptor  detected  (with  one  exception  noted  below) .  Several  Si 
doped  samples  grown  at  Bell  Labs  using  various  growth  temperatures  (Tq  - 
600-6508C)  and  various  Ga  to  As  pressure  ratios  (P^s^Ga  “  15  to  25  with 
Pga  and  ?§£  held  constant)  were  measured,  but  no  strong  trends  were 

observed  in  impurity  incorporation.  Figure  6  is  a  photothermal  ionization 


Photolu  m  i  n  escence 


Pig.  5  Photoluminescence  spectra  for  Toshiba  MOCVD  samples 
showing  decrease  of  Ge  acceptors  with  increasing 
III/V  ratio. 


-162- 


spectrum  which  Is  typical  of  these  samples.  A  snail  peak  not  visible  In 
the  other  Bell  samples  coincident  with  SI  acceptor  (D°-A°)  (see  Fig.  7) 
recombination  was  observed  in  the  sample  grown  with  the  lowest  V^s^Ga’ 
which  is  consistent  with  the  expected  highest  As  site  availability  for  Si 
incorporation. 

8.6  CONCLUSIONS 

In  conclusion,  the  relative  incorporation  of  the  Group  IV  impurities  in 
high  purity  GaAs  prepared  by  a  variety  of  growth  techniques  has  been 
measured.  The  model  of  Teramoto  (1972)  for  LPE  and  its  extension  of 
ASCI3-VPE  by  Ashen  et  al.  (1975)  are  in  qualitative  agreement  with  the 
observations  for  all  the  Group  IV  elements.  In  particular  the  large 

O 

lattice  strain  energy  contributions  to  for  Pb  and  C  are  consistent 

with  the  failure  to  conclusively  observe  C  donors  or  Pb  acceptors  in 
epitaxial  GaAs  grown  by  any  technique.  Finally,  the  III/V  ratio 

dependences  of  Group  IV  incorporation  for  AsHj-VPE,  MDCVD,  and  MBE  are  in 
qualitative  agreement  with  the  expected  influence  on  concentrations  of  Ga 
and  As  vacancy  sites  available  for  impurity  incorporation. 
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9.  IDENTIFICATION  OF  Si  DONORS 

Silicon  is  widely  believed  to  be  an  important  residual  impurity  in  GaAs 
prepared  by  several  grovrth  techniques.  The  results  of  many  experiments  on  th 
growth  of  high  purity  GaAs  have  been  explained  in  terms  of  tha  influence  of 
various  growth  parameters  on  the  incorporation  of  Si.  Most  such  experiments 
have  relied  on  Hall  effect  data  and  on  the  tacit  assumption  that  only  changes 
in  Si  incorporation  are  important  in  determining  changes  in  the  carrier 
concentration.  Since  Si  is  neither  the  only  residual  acceptor^  nor  the  only 
residual  donor  present  in  high  purity  GaAs  grown  by  a  wide  variety  of 
techniques,  it  is  necessary  to  specifically  measure  the  concentration  of 
electrically  active  Si  in  order  to  assess  the  role  of  Si  as  a  residual 
impurity.  Photothermal  ionization  spectroscopy  can  detect  the  donor  species 
present  in  a  high  purity  semiconductor  sample  and  measure  their  relative 
concentrations  from  the  amplitudes  of  the  corresponding  spectral  peaks. 
Carefully  controlled  doping  experiments  have  been  performed  by  several 
research  groups  to  identify  various  donor  species  with  the  associated 
photothermal  ionization  peaks.  Such  experiments  are  difficult  because  of  the 
small  range  of  donor  concentrations  between  that  of  the  purest  GaAs  which  can 
be  grown  by  a  given  technique  and  that  for  which  impurity  interactions8  and 
other  effects^  degrade  the  spectra  by  broadening  and  distorting  the  spectral 
peaks.  In  this  report  we  discuss  photothermal  ionization  measurements  made  on 
Si  doped  GaAs  grown  by  MBE  which  have  resulted  in  the  identification  of  Si 
donors  with  a  spectral  peak  different  from  that  previously  associated  with 
Si.8 

9.1  PHOTOTHERMAL  IONIZATION  SPECTRA 

Si  doped  samples  1  and  2  were  grown  at  Hewlett-Packard  Laboratories  on 
(100)  undoped  semi-insulating  LEC  substrates  in  a  Varian  GEN/II  MBE  system 
equipped  with  liquid  nitrogen  cryoshrouds.  All  source  materials  were 
evaporated  from  pyrolytic  boron  nitride  (PBN)  crucibles,  contained  in  high 


purity  source  furnaces  in  which  the  heated  zones  are  comprised  of  only  Ta  and 
PBN.  An  elemental  As  source  was  used,  and  undoped  MBE  GaAs  prepared  in  this 
way  is  generally  p-type  with  p  =  lO^cm-^.^  Doping  with  elemental  Si  to 
produce  the  n-type  samples  required  for  the  measurements  produces  a  Si  donor 
concentration  which  dominates  the  small  background  concentration  of  other 
donors  and  acceptors  present  in  undoped  material.  The  only  materials  other 
than  Ga,  As,  and  Si,  which  had  been  previously  evaporated  in  this  MBS  system 
were  Al,  and  Be  (a  p-type  dopant).  The  layers  were  grown  with  a  substrate 
temperature  of  585 °C  at  a  rate  of  1  pm/hr.  The  electrical  properties  of  these 
samples  and  the  other  samples  reported  in  this  work,  corrected  for  surface  and 
substrate  depletion*^,  are  shown  in  Table  1. 

Photothermal  ionization  spectra  for  samples  1  and  2  and  for  reference 
sample  1  are  shown  in  Fig.l.  These  spectra  show  the  ls-2p(m=-l)  transitions 
of  the  various  donor  species  present,  and  were  recorded  at  precisely  the  same 
magnetic  field  (B=2.92T)  by  operating  the  superconducting  magnet  in  the 
persistent  mode.  The  three  peaks  in  the  spectrum  of  the  ultra-pure  ASCI3-H2 
grown  reference  sample  1,  labeled  X^,  X2,  and  X3,  correspond  to  the  three 
characteristic  residual  donors  present  in  GaAs  grown  by  this  technique.  At  a 
given  magnetic  field,  the  energies  of  these  peaks  provide  a  reference  for 
identifying  the  peaks  in  the  spectra  of  other  samples. 

The  spectra  of  both  Si  doped  samples  1  and  2  contain  two  peaks  which  have 
the  typically  observed  asymmetric  lineshape  predicted  by  Larsen.^  In  both 
spectra,  the  large  peak  corresponds  with  In  the  reference  spectrum,  while 
the  smaller  peak  corresponds  with  X2»  The  slight  downshift  in  the  peak 
energies  of  the  MBE  spectra  with  respect  to  those  in  the  reference  spectrum  is 
often  observed  in  camples  of  this  doping  level.  The  dominance  of  the  peak 
in  the  spectra  for  these  Si  doped  MBE  samples  by  itself  strongly  suggests  the 


Table  1  Hall  effect  data  and  impurity  concentrations  derived 

from  spectra  and  Hall  mobility  analysis  of  Wolfe  et  al.  . 


Impurity  Concentrations  (xlO^cm-^) 


Sample 

ND( total) 

NA( total) 

ND(Pb) 

ND(X1=Si) 

ND(Sn) 

nd(x2) 

nd(x3) 

1 

9.7 

3.2 

- - r_ 

7.3 

2.4 

2 

11.3 

3.7 

— 

9.0 

— 

2.3 

- - 

3 

2.8 

1.6 

*0.1 

*0.1 

1.3 

1.3 

— , — 

reference  1 

0.89 

0.45 

— 

0.47 

— 

0.36 

0.  06 

reference  2 

0.43 

0.13 

- - 

0.8 

0.14 

0.14 

0.08 

.31  32  33 

Wavenumber  (cm-1) 


Photothermal  ionization  spectra  of  the 
ls-2p (m=-l)  transitions  for  the  Si  doped 
MBE  samples  1  and  2,  and  for  reference 
sample  1. 


association  of  Si  donors  with  the  peak,  but  further  evidence  is  provided  by 
the  absolute  donor  concentrations  in  these  samples.  The  relative 
concentrations  of  and  X2  donors  were  derived  from  the  relative  peak 
amplitudes  in  the  spectra.  Absolute  concentrations  could  then  be  calculated 
from  the  total  donor  concentration  using  the  77K  Hall  effect  data  and  the 
method  of  Wolfe  et  al.*^  In  the  more  heavily  doped  sample  2,  the 
concentration  of  X^  donors  is  1.23  times  larger  than  in  sample  1,  but  the 
concentration  of  X2  donors  is  essentially  the  same  in  both  samples.  This  is 
interpreted  to  mean  that  the  concentration  of  X^=Si  donors  increases  with 
increasing  doping  while  the  background  of  X2  donors  remains  constant.  It  is 
also  interesting  to  note,  in  view  of  the  amphoteric  nature  of  Si  in  Gaks,  that 
the  total  acceptor  concentration  in  the  more  heavily  doped  sample  2  is  1.16 
times  that  of  sample  1.  The  identification  of  the  X^  peak  with  Si  is 
consistent  with  the  work  of  Ozeki  et  al.^ 

The  data  which  led  to  the  earlier  identification  of  peak  X2  with  Si  by 

O 

Wolfe  et  al.  are  shown  in  Fig.  2.  In  these  experiments  the  photothermal 
ionization  spectrum  of  a  Sn  doped  ASCI3-H2  VPE  grown  GaAs  control  sample  was 
compared  to  that  of  a  sample  which  was  prepared  identically  except  that  Si  was 
added  to  the  Ga  melt.  On  the  basis  of  the  larger  amplitude  and  width  of  the 
X2  peak  in  the  Si  doped  sample  the  X2  peak  was  identified  with  Si. 

The  new  identification  of  Si  with  X^  can  be  reconciled  with  the  data  of 

O 

Wolfe  et  al.,  by  examining  the  region  of  the  spectra  of  Fig.  2  in  the 
vicinity  of  X^.  While  the  X2  peaks  have  the  typically  observed  low  energy 
tail  predicted  by  Larsen,^  the  shapes  of  the  X^  peaks  in  both  spectra  are 
anomalous.  The  photoresponse  drops  precipitously  for  energies  less  than  that 
of  the  X^  peaks,  and  after  a  sharp  notch,  forms  the  broad  wing  noted  by  Wolfe 

O 

et  al.  A  careful  comparison  of  the  original  spectra  for  these  samples  with  a 


spectrum  of  reference  sample  2,  in  which  all  peaks  have  the  typically  observed 
shape,  shows  that  the  peak  in  the  reference  spectrum  corresponds  not  with 
the  peaks  labeled  in  the  spectra  in  Fig. 2,  but  with  the  notches  in  the 
photoresponse  just  below  these  peaks.  Photothermal  ionization  peaks  having  a 
sharply  peaked  component  slightly  above  the  usual  peak  energy  and  a  broad  wing 
well  below  the  usual  peak  energy  have  been  observed  previously,  particularly 
for  the  larger  amplitude  peaks  (i.e.  higher  concentration  donors).  Space  does 
not  permit  a  thorough  discuss xon  of  this  effect  and  related  experiments,*^  and 
the  details  of  the  phenomenon  are  not  yet  fully  understood,  but  we  believe  the 
anomalous  peak  shape  is  caused  by  an  absorption  by  the  corresponding  (in  this 
case  X^=Si)  donors  near  the  substrate-epilayer  interface,  which  does  not 
contribute  significantly  to  the  photoconductive  response,  but  reduces  the 
intensity  of  far  infrared  light  within  the  sample  at  the  peak  energy,  so  that 
a  sharp  notch  is  observed  in  the  associated  photothermal  ionization  peak.  In 
the  absence  of  this  interface  absorption,  the  peaks  in  both  spectra  would 
have  the  more  typical  shape  of  X£  and  substantially  larger  amplitudes,  indica¬ 
tive  of  the  actual  relative  concentrations  of  Xj=Si.  The  greater  width  of  the 
X2  peak  in  the  Si  doped  sample,  as  well  as  that  of  the  broad  wing  below  X^, 
which  with  the  present  interpretation  is  the  usual  low  energy  tail  of  X^ 
distorted  by  interface  absorption,  are  the  result  of  additional  Stark  broaden- 

c  * 

ing  caused  by  the  additional  charged  centers  (ionized  acceptors  and 

compensated  donors)  introduced  by  the  Si  doping. 

9.2  RESIDUAL  Si  DONORS 

The  identification  of  Si  with  X^  instead  of  X£  in  photothermal  ionization 
spectra  has  important  implications  with  regard  to  the  role  of  Si  as  a  residual 
donor  in  GaAs  grown  by  liquid  phase  epitaxy  'LPE).  Until  recently  was 
never  observed  in  LPE  GaAs,  and  this  led  to  the  suggestion  that  this  peak 
might  be  due  to  a  stoichiometric  defect.*^  The  only  residual  donor  peaks 


-173- 


previously  observed  la  LPE  materials  were  identified  as 

and  X2,  which  is  often  the  dominant  peak  in  the  spectra,  and  which  has  been 
associated  with  The  earlier  identification  of  Si  with  X2  fit  well  with 

the  expectation  that  reduction  of  hot  Sif^  by  H2  in  the  LPE  reactor  would  lead 
to  incorporation  of  Si  into  the  melt,  and  subsequently  into  the  growing 
epilayer,  *  and  that  this  mechanism  was  the  dominant  contribution  to 
residual  donors  in  LPE  GaAs.  The  present  identification  of  Si  with  X-^, 

o  5  IQ 

together  with  the  previous  photothermal  ionization  data  for  LPE  GaAs,  »  » 
indicate  not  only  that  Si  is  not  the  dominant  residual  donor,  but  that  it  is 
often  present  in  undetectably  small  concentrations  relative  to  the  other 
residual  donor  species.  Figure  3  shows  the  spectrum  cf  the  high  purity  LPE 
GaAs  sample  3,  which  was  grown  at  Cornell  University.  Peaks  corresponding  to 
the  usual  Pb,  Sn  and  X£  donors  are  clearly  visible  in  this  spectrum,  but  a 
small  peak  at  the  energy  of  X^=Si  was  also  reproducibly  observed.  Most  of  the 
photoresponse  at  the  energy  of  is  due  to  the  low  energy  tail  of  the  Sn  peak 
and  so  the  relative  concentration  of  X^  is  quite  small.  The  proximity  of  the 
X^  and  Sn  peaks  in  energy  and  the  fact  that  the  residual  Sn  concentration  is 
usually  much  larger  than  the  concentration  of  X^  donors  observed  here  may 
account  for  the  fact  that  this  is  the  first  report  of  the  X^  peak  in  the 
photothermal  ionization  spectra  of  undoped  LPE  GaAs. 

The  identification  of  Si  with  X^  instead  of  X2  also  has  implications  with 

regard  to  the  role  of  Si  as  a  residual  donor  in  GaAs  grown  by  various  other 

techniques.  In  AsH<j  VPE  grown  GaAs  the  residual  donors  are  X^,  X2,  and  X-j,  as 

in  ASCI3-H2  VPE,  but  the  concentrations  of  X^  and  X3  are  lower,  and  X2  is 
9 

dominant  .  In  MBE  GaAs  grown  with  an  elemental  As  source  and  using  Sn  as  the 
n-type  dopant,  both  X^  and  X2  are  present  but  the  concentration  of  Xj^  is 
low,  while  if  A3H3  is  cracked  in  a  SiC^  furnace  to  provide  the  As,  the 


Photoresponse 


amplitude  of  the  X^=Si  peak  increases  with  increasing  cracking  furnace 

temperature.  In  MOCVD  grown  GaAs  the  relative  concentrations  of  various 

residual  donors  including  X^  and  X2  deduced  from  pliotothermal  ionization 

measurements  have  been  correlated  with  various  growth  parameters  and 

purification  techniques  from  the  metalorganics.^  These  results  indicate  that 

in 

the  residual  Xi=Si  donor  concentration  MOCVD  material  can  be  made  negligible. 

In  conclusion,  on  the  basis  of  photothermal  ionization  measurements  on  Si 
doped  MBE  GaAs  samples,  an  identification  of  Si  with  the  spectral  peak  Xj^  has 
been  made.  This  identification  together  with  previous  photothermal  ionization 
data  indicate  that  Si  is  not  an  important  residual  donor  in  high  purity  LPE 
GaAs,  and  also  necessitates  a  reevaluation  of  the  importance  of  Si  as  a 
residual  donor  in  GaAs  grown  by  a  variety  of  other  techniques. 
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10.  RESIDUAL  DONORS  IN  MOLECULAR  BEAM  EPITAXY 
Molecular  beam  epitaxy  (MBC)  of  I1I-V  semiconductors  is  of  great 
interest  for  semiconductor  device  applications  because  of  its  capability 
for  producing  highly  uniform  layers  with  precisely  controlled  thicknesses 
and  doping  profiles.  For  some  applications,  such  as  FET  buffer  layers 
and  undoped  layers  in  high  electron  mobility  transistors,  it  is  also 
important  to  be  able  to  grow  high  purity  material.  In  this  report  we 
discuss  photothermal  ionization  measurements  of  the  shallow  donor  species 
present  in  high  purity  MBE  GaAs  samples  grown  in  two  different  laboratories 
Photothermal  ionization  spectroscopy^  is  a  technique  which  can  be 
used  to  detect  the  donor  species  present  and  measure  their  concentrations 
relative  to  the  total  donor  concentration  in  high  purity  semiconductor 
samples.  The  extrinsic  photoconductivity  spectrum  of  a  high  purity  n-type 
GaAs  sample  consists  of  sharp  peaks  at  energies  corresponding  to  transition 
energies  of  the  hydrogenic  donors.  Because  different  donor  species  have 
slightly  different  ground  state  energies  but  nearly  identical  excited 
state  energies,  each  hydrogenic  transition  has  associated  with  it  a 
closely  spaced  multiplet  of  peaks.  Each  peak  in  the  multiplet  corresponds 
to  a  hydrogenic  transition  of  a  particular  donor  species.  The  relative 
amplitudes  of  the  peaks  are  a  measure  of  the  relative  concentrations  of 
the  corresponding  donor  species.  The  narrow  linewidth  and  large  amplitude 
of  the  ls-2p(m=-l)  transition  at  high  magnetic  fields  make  it  especially 
suitable  for  resolving  the  closely  spaced  peaks  of  the  individual  donor 
species.  The  ls-2p(m~-l)  transition  energies  change  substantially  with 
changing  magnetic  field  so,  to  permit  comparison  of  the  peaks  observed 
in  a  given  sample  with  donor  identifications  in  the  literature,  the 


spectrum  of  that  sample  is  compared  to  the  spectrum  of  a  well  characterized 

reference  sample  taken  at  precisely  the  same  magnetic  field.  This  is 

accomplished  by  operating  the  superconducting  magnet  in  persistent  mode 

while  both  spectra  are  recorded.  The  reference  sample  used  in  this  work 

is  an  ultra  pure  AsCl^-VPE  GaAs  sample  grown  by  C.M.  Wolfe,  which  has  a 

2 

liquid  nitrogen  temperature  mobility  of  p  ^=201,000  era  /Vs  and  a  carrier 

13  -3  2 

concentration  of  n^=4.5xl0  cm  .  The  three  peaks  present  in  the 

spectra  of  this  sample,  labeled  X^  and  X^,  correspond  to  the  three 

characteristic  residual  donors  present  in  high  purity  GaAs  grown  by  this 
3 

technique.  Since  many  of  the  early  donor  identification  experiments 
were  done  using  AsCl^-VPE  material,  the  energies  of  peaks  corresponding 
to  different  donor  species  relative  to  the  energies  of  X^  X^  and  X^  are 
fairly  well  known. 

10.1  EPITAXIAL  SAMPLES  6  7 

Hall  effect  data  for  the  KBE  samples  described  here  ’  show  them  to 

be  the  highest  purity  GaAs  prepared  by  this  technique  which  has  yet  been 

reported.  Sample  A-147  was  grown  by  H.  Morkoq  and  A.Y.  Cho  who  reported 

2 

a  liquid  nitrogen  temperature  mobility  cf  105,000  cm  /Vs  and  a  carrier 

14  -3  5 

concentration  of  4x10  cm  for  the  27  pm  thick  layer.  The  MBE  system 

used  to  grow  this  sample  used  an  elemental  As  source,  and  undoped  GaAs 

14  -3 

grown  in  this  way  is  generally  p-type  with  p;10  cm  .  In  order  to 
produce  an  n-type  sample  it  was  necessary  to  lightly  dope  the  sample  with 
Sn.  The  details  of  the  substrate  preparation,  growth  procedure  and 

6  8 

precautions  for  obtaining  high  purity  have  been  described  elsewhere.  * 
Samples  MBE-91  and  MBE-92  were  grown  by  A.R.  Calawa^  in  an  KBE  system 
which  used  an  AsH„  cracking  furnace  in  place  of  the  usual  effusion  cell 
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As  source.  Undoped  samples  grown  in  this  way  are  n-type  and  so  no 

intentional  doping  was  required.  The  details  of  the  cracking  furnace 

and  growth  procedure  for  these  samples  have  also  been  described  previously.^ 

The  liquid  nitrogen  temperature  mobilities,  and  carrier  concentrations 

corrected  for  both  surface  and  interface  depletion  in  these  5  pm  thick 

samples  are^*^  p^=90,000  cm^/Vs,  n^=2.7xl0^  cm  ^ ,  and  y^=110,000  cra^/Vs, 
14  -3 

n^=2.4xl0  cm  respectively. 
lQ  2  PHOTOTHERMAL  IONIZATION  SPECTRA 

'  For  extremely  pure  samples  in  which  the  photothermal  ionization  peaks 
are  separated  by  more  than  a  linewidth  in  energy,  the  relative  amplitudes 
of  the  peaks  are  expected  to  correspond  to  the  relative  concentrations  of 
donor  species  present.  In  less  pure  and/or  more  compensated  samples,  as 
the  concentration  of  charged  centers  (ionized  acceptors  and  compensated 
donors  in  an  n-type  sample)  increases,  the  peaks  broaden  due  to  the  electric 
fields  and  field  gradients^  from  these  charged  centers,  until  the  line- 
width  becomes  comparable  to  the  peak  separations.  The  photoresponse  at 
the  energy  of  a  peak  of  a  certain  donor  species  is  then  no  longer  simply 
proportional  to  the  concentration  of  that  species,  but  is  related  also  to 
the  concentration  of  donors  associated  with  adjacent  peaks.  This  is  the 
case  for  the  spectrum  of  MBE-91  in  Fig. 1(c).  Because  the  photoresponse  at 
the  position  of  the  Pb  peak  is  largely  due  to  the  contributions  of  the 
low  energy  tails  of  the  Si  and  Sn  peaks,  the  relative  amplitude  of  the 
Pb  peak  is  substantially  larger  than  the  relative  concentration  of  Pb 
donors.  In  order  to  obtain  a  quantitative  estimate  of  the  concentrations 
of  the  various  donor  species  for  each  sample  reported  here,  a  function  f(e) 

4 

f(e)  “  I  a.P(0[e-e  J) 
i=.l 


of  the  form 


Fig.  1  Photothermal  ionization  spectra  of  the  ls-2p(m*-l) 
transitions  for  the  MBE  samples  and  the  reference 
sample  at  a  magnetic  field  of  2.92T. 
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was  fitted  to  the  envelope  of  ls-2p(m=-l)  transitions  using  a  least  squares 
minimization  technique.  Each  term  in  f(e)  is  a  peaked  function  whose 
amplitude  a  is  interpreted  as  a  relative  concentration  of  the  i*"^  donor 
species.  The  function  P(e)  was  derived  from  a  fit  to  the  ls-2p(m=-l) 
lineshape  calculation  by  Larsen. ^  The  values  of  a ^  and  B  were  adjusted 
for  a  best  fit  while  the  peak  positions  were  determined  from  the  spectra 
of  Fig. 1(a)  and  (d) .  The  width  scale  B  was  constrained  to  be  the  same  for 
each  peak  because  the  inhomogeneous  broadening"^  of  peaks  is  expected  to 
be  the  same  for  each  donor  species  in  a  uniformly  doped  sample,  independent 
of  its  relative  concentration.  The  resulting  best  f(e)  and  the  individial 
terms  of  f(e)  are  plotted  for  MBE-91  along  with  the  experimental  data  in 
Fig. 1(c).  The  total  donor  and  acceptor  concentrations  for  the  samples 
reported  here  were  determined  using  the  Hall  mobility  analysis  of  Wolfe 
et  al."^  The  relative  concentrations  for  the  various  donor  species  present 
in  each  sample  could  then  be  expressed  as  the  absolute  concentrations  shown 
in  Table  1. 

Photothermal  ionization  spectra  for  sample  A-.147  and  for  the  reference 
sample,  both  recorded  at  a  magnetic  field  of  4.97  T,  are  shown  in  Fig. 2. 

As  expected  from  the  intentional  Sn  doping  of  this  sample,  the  dominant  peaks 
in  its  spectrum  occurs  at  an  energy  which  has  been  previously  associated 
with  Sn  by  Wolfe  et_  al.  *  and  Cooke  ej^  _al.  ^  The  smaller  peak  at  higher 

energy  corresponds  to  X^.  The  peak  was  originally  associated  with  Si 

4  12 

by  Wolfe  al^.  but  subsequent  doping  experiments  by  Ozeki  et_  al_.  and 

13  12 

in  our  laboratory  have  indicated  that  X^  is  Si,  while  Ozeki  et  al .  have 

associated  X^  with  S.  The  shoulder  on  the  low  energy  side  of  the  dominant 

Sn  peak  occurs  at  an  energy  which  has  been  associated  with  Pb  by  Wolfe 
4  14 

et  al.  and  Stradling.  In  a  homogeneously  doped  sample,  the  ratio  or 


pe.nk  width  to  peak  height  should  be  the.  same  for  all  peaks.  This  ratio, 


»ri 
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Table  1  Impurity  concentrations  derived  from  the  photothermal 
ionization  spectra  and  the  Hall  mobility  analysis  of 
Wolfe  et  al. 


Sample 

rt 

O 

rt 

& 

M 

N  ( total) 
A 

ND(Pb) 

VXl“Si> 

ND(Sn) 

VV 

nd(x3) 

Reference 

0.89 

0.45 

— 

0.47 

— 

0.36 

0.06 

A-147 

5.8 

1.8 

<0.3 

<0.3 

3.9 

1.3 

<0.1 

MBE-91 

4.9 

3.9 

0.9 

2.5 

0.8 

0.7 

<0.2 

MBE-92 

3.4 

2.8 

0.8 

<0.15 

1.3 

1.2 

— 

Wavenumber  (cm  ’) 


Fig.  2  Photothermal  ionization  spectrum  of  the 

ls-2p(m=-l)  transitions  for  MBE  sample  A-147 
and  the  reference  sample  at  a  magnetic  field 
of  4.97T. 
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taken  at  75%  of  the  full  peak  heights,  is  1.7  times  larger  for  the  Sn  peak 
than  the  X^  peak,  and  this  suggests  that  another  unresolved  donor  (possibly 
X^=Si)  may  be  contributing  to  the  low  energy  tail  of  the  Sn  peak. 

Photothermal  ionization  spectra  for  samples  MBE-91,  MBE-92,  A-147, 
and  the  reference  sample,  recorded  at  a  magnetic  field  of  2.92  T,  are  shown 
in  Fig.l.  The  maximum  photoresponse  in  the  spectrum  of  MBE-91  is  at  the 
energy  of  X^=Si,  but  if  only  Si  were  present  in  this  sample,  the  photoresponse 
would  drop  off  more  sharply  for  energies  above  that  of  the  maximum.  The 
additional  photoresponse  on  the  high  energy  side  of  the  Si  peak  is  probably 
due  to  Sn,  especially  in  view  of  the  presence  of  a  Sn  peak  in  the  spectrum 
of  MBE-92.  At  still  higher  energy  there  is  a  small  peak  coincident  in 
energy  with  X 2  and  possibly  also  with  X^,  while  at  lower  energy  there  is 
a  peak  at  the  energy  associated  with  Pb. 

The  photothermal  ionization  spectrum  for  MBE-92  in  Fig. 1(d)  shows  the 
clear  presence  of  Pb  Sn  and  X^,  and  the  fact  that  the  Sn  peak  is  slightly 
broader  relative  to  its  amplitude  than  either  the  Pb  or  X^  peak,  suggest 
there  may  be  a  slight  amount  of  X^=Si  present.  The  presence  of  Sn  in  the 
unintentionally  doped  samples  MBE-91  and  MBE-92  may  be  due  in  part  to  the 
previous  use  of  Sn  as  an  n-type  dopant  in  the  MBE  growth  chamber.  Several 
growth  runs  were  made  between  the  use  of  Sn  and  the  growth  of  MBE-91  and 
MBE-92,  however,  and  so  passivation  of  the  walls  of  the  growth  chamber 
may  have  minimized  this  contribution  of  Sn.  Samples  MBK-91  and  MBE-92  were 
grown  under  identical  conditions  except  for  the  temperature  of  the  Aslt^ 
cracking  furnace.  Calawa^  correlated  an  increase  in  net  carrier  concentra¬ 
tion  for  samples  grown  with  increasing  cracking  furnace  temperature  and 
attributed  it  to  Si  donors  provided  by  the  hot  quartz  walls  of  the  cracking 
furnace.  This  idea  is  substantiated  by  the  donor  concentrations  derived 
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from  the  fits  to  the  spectra  of  MBE-91  and  MBE-92.  In  MBE-91,  grown  with 

a  cracking  furnace  temperature  of  690°C>  the  concentration  of  X^=Si  donors 

is  over  half  of  the  total  donor  concentration,  while  in  MBE-92,  grown  with 

the  lower  cracking  furnace  temperature  of  610°C,  and  having  a  lower  total 

donor  concentration,  the  presence  of  X^=Si  donors  is  barely  discernable. 

The  concentrations  of  Pb,  Sn  and  X^  donors,  however,  are  similar  for  both 

samples  (see  Table  1.)* 

10.3  CONCLUSIONS 

In  conclusion,  the  photothermal  ionization  data  presented  here  for 

high  purity  MBE  GaAs  samples  grown  in  different  laboratories  and  using 

substantially  different  growth  techniques  show  remarkably  similar  donor 

backgrounds.  Peaks  associated  with  Pb,  Si,  Sn,  and  S  are  present  in  each 

of  the  MBE  samples  measured.  The  X^  donor,  if  present  at  all  in  these 

samples,  has  a  concentration  at  least  a  factor  of  20  smaller  than  the  total 

donor  concentration.  The  absence  of  Xy  which  has  been  associated  with  both 

and  Ge*'**^  in  the  literature,  is  interesting  because  it  is  one  of  the 

characteristic  residual  donor  species  present  in  AsCl^-H^  VPE,  AsCl^-N^  VPE, 

3  10 

and  AsH^  VPE  grown  GaAs  *  and  it  is  the  dominant  residual  donor  species 

3  5  17  .  . 

in  MOCVD  grown  GaAs.  *  *  In  contrast,  the  Pb  donor  is  absent  in  uninten¬ 
tionally  doped  GaAs  grown  by  all  of  the  above  techniques  except  MBE. 
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11.  MULTIPLE  DOPING  EXPERIMENTS 

In  this  section  we  present  the  results  of  a  study  under¬ 
taken  to  determine  the  interactions  between  or  complexing  of 
multiple  dopants  in  GaAs.  The  impurities  initially  studied 
are  sulfur  and  oxygen.  To  examine  possible  interactions, 
experiments  involving  the  growth  and  characterization  of  sam¬ 
ples  homogeneously  and  inhomogeneous ly  doped  with  H20  and/or 
H2S  were  undertaken.  The  inhomogeneous ly  doped  layers  provide 
self-calibrating  samples  for  subsequent  analysis. 

11.1  HOMOGENEOUSLY  DOPED  LAYERS 

All  of  the  homogeneously  doped  layers  were  grown  on  Cr 
doped,  semi-insulating  substrates  oriented  2°  off  {100} 
towards  {110}.  An  extensive  series  of  growth  runs  were  per¬ 
formed  using  H2S(g)  diluted  to  11  ppm  as  the  only  doping 
source.  This  series  of  runs  yielded  the  expected  result  that 
sulfur  is  a  well-behaved  shallow  donor  in  GaAs.  The  variation 
of  carrier  concentration  (given  by  the  effective  distribution 

coefficient)  with  H2S  flow  rate  is  shown  in  Fig.  1.  Doping 

15  -3  17  -3 

levels  from  3x10  cm  to  3x10  cm  were  reproducibly 
obtained. 

Another  series  of  runs  was  performed  using  only  H2Q  as  the 
doping  source.  The  input  partial  pressure  of  H20  was  controlled 
by  a  constant  temperature  bath.  The  dopant  was  transported 
into  the  reactor  by  passing  H2  o^er  a  bed  of  solid  H20  for 
those  runs  performed  with  the  temperature  bath  at  less  than 
0°C,  and  by  bubbling  H2  through  liquid  H20  for  those  runs 
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Figure 


performed  with  the  temperature  bath  at  greater  than  0°C. 

The  results  of  the  H20  doping  for  the  temperature  range  between 
-78°C  and  +24°C  were  quite  striking.  For  each  run  the  result¬ 
ing  layer  was  highly  compensated  with  carrier  concentrations 
ranging  from  lxlO^cm  ^  at  -78°C  to  8xl0^cm  3at  24°C. 

Growth  runs  were  also  performed  in  which  the  layers  were 
simultaneously  doped  with  H2S  and  H20.  A  series  of  runs  was 
performed  to  determine  the  effect  of  H2  flow  rate  over  the 
H20(s).  The  H2S  flow  rate  used  was  7.5  ml/min  and  the  H20  source 
was  maintained  at  a  temperature  of  -78°C.  Fig.  2  shows  the 
effect  that  varying  the  gas  flow  over  the  H20(s)  has  on  the 
carrier  concentration  obtained  from  van  der  Pauw1  measurements 
at  300K  and  77K.  From  these  data,  two  effects  are  demonstrated. 

The  first  effect  to  be  noticed  is  the  increase  in  carrier 
concentration  with  an  increase  in  flow  rate  over  the  H20.  The 
second  effect  can  be  seen  in  the  freeze-out  of  carriers  at 
77K.  The  last  column  of  Fig.  2  shows  the  difference  between 
the  carrier  concentrations  at  300K  and  77K  for  each  of  the  flow 
rates  given.  This  shows  that  the  concentration  of  carriers 
freezing  out  is  linearly  dependent  on  the  flow  rate.  Such 
behavior  is  indicative  of  a  deep  donor. 

The  effect  of  increasing  the  partial  pressure  of  H20  in 
the  system  is  similar  to  increasing  the  flow  rate  of  H2  through 
the  H20,  although  a  definite  upper  limit  to  the  obtainable 
carrier  concentration  is  evidenced.  Layers  were  grown  with 
an  H2S  flow  of  7.5  ml/min  and  a  H2+H20  flow  of  5.0  ml/min. 

The  temperature  of  the  H20  was  varied  from  layer  to  layer 
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over  the  range  from  -78°C  to  24°C.  Fig.  3  shows  the  300K 
carrier  concentration  versus  H20  temperature  and  partial 
pressure  .  At  the  low  end  of  the  temperature  range  (-78°C) , 
the  carrier  concentration  is  about  3xl0^cm-^  which  is  the 
same  as  for  layers  grown  in  the  absence  of  H20.  In  the  high 
temperature  range  the  carrier  concentration  is  relatively 
insensitive  to  changes  in  H20  temperature  and  remains  constant 
at  about  2x10 ^7cm  ^ . 

11.2  SELF-CALIBRATING  LAYERS 

The  self  calibrating  layer  structure  used  in  this  work  is 
shown  in  Fig.  4.  ND1  and  ND2  are  the  dopants  added  to  the 
layer;  in  this  case  ND1  and  ND2  are  due  to  H2S  and  H20, 
respectively.  The  layer  is  structured  into  five  regions. 
Regions  I  and  V  are  the  calibration  regions  in  which  the 
chemical  and  electrical  properties  of  each  dopant  may  be  deter¬ 
mined  independent  of  any  interactions  with  the  other  dopant. 
Regions  II  and  IV  are  the  transition  regions  for  dopants  D2 
and  Dl,  respectively,  which  contain  the  information  on  how 
the  chemical  and  electrical  concentrations  change  as  a 
function  of  concentration  of  the  dopant  being  switched  on  or 
off.  Region  III  is  the  region  in  which  both  dopants  are  at 
their  maximum  levels  and  the  effects  of  any  interactions  or 
complexing  should  be  greatest.  Routine  analysis  of  self¬ 
calibrating  layers  consists  of  SIMS  profiling  to  determine 
the  chemical  composition  of  the  layer  and  C-V  profiling  to 
determine  the  free  carrier  concentration.  The  C-V  analysis  was 
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Figure  4.  Self-calibrating  layer  structure. 
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performed  using  an  electrochemical  cell  so  that  thick, 
heavily  doped  layers  could  be  examined4 ' ^ ® ' 7 ' ® . 

The  self-calibrating  layers  were  grown  on  Si  doped  sub¬ 
strates  oriented  2°  off  {100}  towards  {110}  with  a  carrier 

17  —3  18  “3 

concentration  between  2x10  cm  and  2x10  cm  .  The  H2S 
flow  rate,  when  on,  was  set  for  7.5  ml/min.  The  flow  of 
H2+H20  was  5.0  ml/min  with  the  H20  temperature  ranging  from 
-78°C  to  24°C.  Figs.  5,  6,  and  7  show  the  SIMS  and  C-V  analy¬ 
ses  for  three  samples  grown  with  various  input  partial  pressures 
of  H20.  In  these  figures  n^  is  the  intrinsic  concentration 

g 

at  the  growth  temperature  .  In  the  H2S  calibration  region  of 
all  three  samples,  the  free  carrier  concentration  at  the  growth 
temperature  (740°C)®  is  between  lxlO^cm  1  and  2xl01^cm  ^ 
which  is  an  order  of  magnitude  higher  than  the  sulfur  concentra¬ 
tion.  This  may  be  due  to  the  influence  of  a  slow  turn-off 
profile  for  the  H20  as  shown  in  Figs.  6  and  7.  A  possible 
effect  of  the  gradual  turn-off  of  H20  may  be  seen  in  the 
slope  of  the  sulfur  profiles  in  the  H2S  calibration  region 
in  Figs.  6  and  7.  Samples  grown  without  H20  present  or  with 
a  very  low  partial  pressure  of  H20  exhibit  flat  sulfur  pro¬ 
files  in  the  calibration  region.  There  are  two  possible 
reasons  for  the  dependence  of  the  sulfur  profile  on  the  H20: 

(1)  the  distribution  or  the  diffusion  of  sulfur  in  GaAs  is 
altered  through  some  interaction  with  the  H20  or  (2)  oxygen 
incorporated  in  the  layer  is  detected  in  the  form  of  02  which 
is  indistinguishable  from  S  in  normal  SIMS  analysis1®. 
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Figure  7.  Self-calibrating  layer  grown  with 
HnO  at  -40  °C. 
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Experiments  with  the  homogeneously  doped  samples  indicate 
that  the  first  explanation  is  more  plausible. 

Examination  of  the  behavior  of  the  sulfur  in  the 
sulfur  transition  region  shows  a  dependence  of  the  sulfur 
turn-on  profile  on  the  input  partial  pressure  of  the  H20. 

From  the  SIMS  analysis,  it  is  evident  that  the  turn-on  profile 
for  the  sulfur  becomes  more  gradual.  Associated  with  the  more 
gradual  turn-on  of  sulfur  is  the  elevation  of  the  sulfur  con¬ 
centration  in  the  H20  and  H2S  calibration  regions. 

Some  increase  in  the  sulfur  concentration  with  increased 
H20  partial  pressure  in  the  sulfur  calibration  region  was 
expected  from  the  work  done  with  the  homogeneously  doped 
layer  and  this  was  observed.  The  increase  of  sulfur  in  the 
H20  calibration  region,  however,  was  unexpected.  Fig.  5 
shows  the  sulfur  concentration  in  the  H20  calibration  region 
as  1x10 ^cm  ^  and  Fig.  6  shows  an  order  of  magnitude  increase 
of  sulfur  for  an  increase  in  H20  temperature  of  10°C. 

There  are  two  possible  explanations  for  this  behavior.  One 
explanation  is,  as  before,  that  02  may  be  mistaken  for  S  in 
the  SIMS  analysis.  The  other  possibilty  is  reactive  diffusion 
of  sulfur  with  the  impurity  introduced  or  defect  induced  by 
the  H20.  The  corresponding  softening  of  the  sulfur  turn-on 
profile  favors  the  latter  explanation. 

Finally,  notice  must  be  taken  of  the  large  concentration 
of  Si  throughout  the  layers  of  Figs.  6  and  7.  It  is  believed 
that  the  source  of  Si  is  the  substrate  and  that  the  profile 
is  the  result  of  autodoping  or  H20  enhanced  outdif fusion. 


The  C-V  profiles  confirm  that  Si  is  the  major  contributor 
to  the  carrier  concentration  when  incorporated  at  such  a 
high  concentration. 

11.3  DISCUSSION 

Experiments  performed  with  homogeneously  doped  layers 
demonstrate  that  H20  introduces  some  impurity  or  defect  into 
the  growing  layer  and  that  there  are  two  interactions  which 
should  be  considered.  These  are  (1)  the  K^O-residual  or 
intentionally-added  impurity  interaction  responsible  for  the 
highly  compensated  H20  doped  samples  and  deep  donor  behavior  and 
(2)  the  H2O-H2S  interaction  responsible  for  the  change  in 
the  sulfur  distribution  or  diffusion.  Analysis  of  self- 
calibrating  layers  yields  more  evidence  that  H20  increases 
the  distribution  or  diffusion  of  sulfur  in  GaAs  and,  at 
high  input  partial  pressures,  contributes  significantly  to 
the  free  carrier  concentration.  Experiments  are  currently 
being  performed  to  separate  growth  incorporation  from  reactive 
diffusion  effects. 
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